|
forrest@0
|
1 |
;************************************************
|
|
forrest@0
|
2 |
; Read ascii, Write nc
|
|
forrest@0
|
3 |
;************************************************
|
|
forrest@0
|
4 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
|
|
forrest@0
|
5 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
|
|
forrest@0
|
6 |
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
|
|
forrest@0
|
7 |
;************************************************
|
|
forrest@0
|
8 |
begin
|
|
forrest@0
|
9 |
; final output
|
|
forrest@0
|
10 |
; diro = "/fis/cgd/cseg/people/jeff/clamp_data/co2/"
|
|
forrest@0
|
11 |
; filo = "co2_globalView_222.nc"
|
|
forrest@0
|
12 |
; c = addfile(diro+filo,"c")
|
|
forrest@0
|
13 |
|
|
forrest@0
|
14 |
; input ascii files
|
|
forrest@0
|
15 |
diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
|
|
forrest@0
|
16 |
b = (/ \
|
|
forrest@0
|
17 |
"aia005_02D2_seas.co2", \
|
|
forrest@0
|
18 |
"aia015_02D2_seas.co2", \
|
|
forrest@0
|
19 |
"aia025_02D2_seas.co2", \
|
|
forrest@0
|
20 |
"aia035_02D2_seas.co2", \
|
|
forrest@0
|
21 |
"aia045_02D2_seas.co2", \
|
|
forrest@0
|
22 |
"aia055_02D2_seas.co2", \
|
|
forrest@0
|
23 |
"aia065_02D2_seas.co2", \
|
|
forrest@0
|
24 |
"alt_01D0_seas.co2", \
|
|
forrest@0
|
25 |
"alt_02D0_seas.co2", \
|
|
forrest@0
|
26 |
"alt_04D0_seas.co2", \
|
|
forrest@0
|
27 |
"alt_06C0_seas.co2", \
|
|
forrest@0
|
28 |
"alt_06D0_seas.co2", \
|
|
forrest@0
|
29 |
"ams_01D0_seas.co2", \
|
|
forrest@0
|
30 |
"ams_11C0_seas.co2", \
|
|
forrest@0
|
31 |
"asc_01D0_seas.co2", \
|
|
forrest@0
|
32 |
"ask_01D0_seas.co2", \
|
|
forrest@0
|
33 |
"avi_01D0_seas.co2", \
|
|
forrest@0
|
34 |
"azr_01D0_seas.co2", \
|
|
forrest@0
|
35 |
"bal_01D1_seas.co2", \
|
|
forrest@0
|
36 |
"bgu_11D0_seas.co2", \
|
|
forrest@0
|
37 |
"bhd_15C0_seas.co2", \
|
|
forrest@0
|
38 |
"bme_01D0_seas.co2", \
|
|
forrest@0
|
39 |
"bmw_01D0_seas.co2", \
|
|
forrest@0
|
40 |
"brw_01C0_seas.co2", \
|
|
forrest@0
|
41 |
"brw_01D0_seas.co2", \
|
|
forrest@0
|
42 |
"bsc_01D0_seas.co2", \
|
|
forrest@0
|
43 |
"car030_01D2_seas.co2", \
|
|
forrest@0
|
44 |
"car040_01D2_seas.co2", \
|
|
forrest@0
|
45 |
"car050_01D2_seas.co2", \
|
|
forrest@0
|
46 |
"car060_01D2_seas.co2", \
|
|
forrest@0
|
47 |
"car070_01D2_seas.co2", \
|
|
forrest@0
|
48 |
"car080_01D2_seas.co2", \
|
|
forrest@0
|
49 |
"cba_01D0_seas.co2", \
|
|
forrest@0
|
50 |
"cba_04D0_seas.co2", \
|
|
forrest@0
|
51 |
"cfa_02D0_seas.co2", \
|
|
forrest@0
|
52 |
"cgo_01D0_seas.co2", \
|
|
forrest@0
|
53 |
"cgo_02D0_seas.co2", \
|
|
forrest@0
|
54 |
"cgo_04D0_seas.co2", \
|
|
forrest@0
|
55 |
"chr_01D0_seas.co2", \
|
|
forrest@0
|
56 |
"cmn_17C0_seas.co2", \
|
|
forrest@0
|
57 |
"cmo_01D0_seas.co2", \
|
|
forrest@0
|
58 |
"coi_20C0_seas.co2", \
|
|
forrest@0
|
59 |
"cpt_36C0_seas.co2", \
|
|
forrest@0
|
60 |
"cri_02D0_seas.co2", \
|
|
forrest@0
|
61 |
"crz_01D0_seas.co2", \
|
|
forrest@0
|
62 |
"csj_06D0_seas.co2", \
|
|
forrest@0
|
63 |
"daa_02D0_seas.co2", \
|
|
forrest@0
|
64 |
"eic_01D0_seas.co2", \
|
|
forrest@0
|
65 |
"esp005_01D2_seas.co2", \
|
|
forrest@0
|
66 |
"esp015_01D2_seas.co2", \
|
|
forrest@0
|
67 |
"esp025_01D2_seas.co2", \
|
|
forrest@0
|
68 |
"esp035_01D2_seas.co2", \
|
|
forrest@0
|
69 |
"esp045_01D2_seas.co2", \
|
|
forrest@0
|
70 |
"esp055_01D2_seas.co2", \
|
|
forrest@0
|
71 |
"esp_02D0_seas.co2", \
|
|
forrest@0
|
72 |
"esp_06D0_seas.co2", \
|
|
forrest@0
|
73 |
"frd040_06C3_seas.co2", \
|
|
forrest@0
|
74 |
"gmi_01D0_seas.co2", \
|
|
forrest@0
|
75 |
"goz_01D0_seas.co2", \
|
|
forrest@0
|
76 |
"gsn_24D0_seas.co2", \
|
|
forrest@0
|
77 |
"haa005_01D2_seas.co2", \
|
|
forrest@0
|
78 |
"haa015_01D2_seas.co2", \
|
|
forrest@0
|
79 |
"haa025_01D2_seas.co2", \
|
|
forrest@0
|
80 |
"haa035_01D2_seas.co2", \
|
|
forrest@0
|
81 |
"haa045_01D2_seas.co2", \
|
|
forrest@0
|
82 |
"haa055_01D2_seas.co2", \
|
|
forrest@0
|
83 |
"haa065_01D2_seas.co2", \
|
|
forrest@0
|
84 |
"haa075_01D2_seas.co2", \
|
|
forrest@0
|
85 |
"hat_20C0_seas.co2", \
|
|
forrest@0
|
86 |
"hba_01D0_seas.co2", \
|
|
forrest@0
|
87 |
"hfm005_01D2_seas.co2", \
|
|
forrest@0
|
88 |
"hfm015_01D2_seas.co2", \
|
|
forrest@0
|
89 |
"hfm025_01D2_seas.co2", \
|
|
forrest@0
|
90 |
"hfm035_01D2_seas.co2", \
|
|
forrest@0
|
91 |
"hfm045_01D2_seas.co2", \
|
|
forrest@0
|
92 |
"hfm055_01D2_seas.co2", \
|
|
forrest@0
|
93 |
"hfm065_01D2_seas.co2", \
|
|
forrest@0
|
94 |
"hfm075_01D2_seas.co2", \
|
|
forrest@0
|
95 |
"hun010_35C3_seas.co2", \
|
|
forrest@0
|
96 |
"hun048_35C3_seas.co2", \
|
|
forrest@0
|
97 |
"hun082_35C3_seas.co2", \
|
|
forrest@0
|
98 |
"hun115_35C3_seas.co2", \
|
|
forrest@0
|
99 |
"hun_01D0_seas.co2", \
|
|
forrest@0
|
100 |
"ice_01D0_seas.co2", \
|
|
forrest@0
|
101 |
"itn051_01C3_seas.co2", \
|
|
forrest@0
|
102 |
"itn123_01C3_seas.co2", \
|
|
forrest@0
|
103 |
"itn496_01C3_seas.co2", \
|
|
forrest@0
|
104 |
"itn_01D0_seas.co2", \
|
|
forrest@0
|
105 |
"izo_01D0_seas.co2", \
|
|
forrest@0
|
106 |
"izo_27C0_seas.co2", \
|
|
forrest@0
|
107 |
"jbn_29C0_seas.co2", \
|
|
forrest@0
|
108 |
"key_01D0_seas.co2", \
|
|
forrest@0
|
109 |
"kum_01D0_seas.co2", \
|
|
forrest@0
|
110 |
"kum_04D0_seas.co2", \
|
|
forrest@0
|
111 |
"kzd_01D0_seas.co2", \
|
|
forrest@0
|
112 |
"kzm_01D0_seas.co2", \
|
|
forrest@0
|
113 |
"lef011_01C3_seas.co2", \
|
|
forrest@0
|
114 |
"lef030_01C3_seas.co2", \
|
|
forrest@0
|
115 |
"lef076_01C3_seas.co2", \
|
|
forrest@0
|
116 |
"lef122_01C3_seas.co2", \
|
|
forrest@0
|
117 |
"lef244_01C3_seas.co2", \
|
|
forrest@0
|
118 |
"lef396_01C3_seas.co2", \
|
|
forrest@0
|
119 |
"lef_01D0_seas.co2", \
|
|
forrest@0
|
120 |
"ljo_04D0_seas.co2", \
|
|
forrest@0
|
121 |
"lmp_28D0_seas.co2", \
|
|
forrest@0
|
122 |
"maa_02D0_seas.co2", \
|
|
forrest@0
|
123 |
"mbc_01D0_seas.co2", \
|
|
forrest@0
|
124 |
"mhd_01D0_seas.co2", \
|
|
forrest@0
|
125 |
"mhdcbc_11C0_seas.co2", \
|
|
forrest@0
|
126 |
"mhdrbc_11C0_seas.co2", \
|
|
forrest@0
|
127 |
"mid_01D0_seas.co2", \
|
|
forrest@0
|
128 |
"mlo_01C0_seas.co2", \
|
|
forrest@0
|
129 |
"mlo_01D0_seas.co2", \
|
|
forrest@0
|
130 |
"mlo_02D0_seas.co2", \
|
|
forrest@0
|
131 |
"mlo_04D0_seas.co2", \
|
|
forrest@0
|
132 |
"mnm_19C0_seas.co2", \
|
|
forrest@0
|
133 |
"mqa_02D0_seas.co2", \
|
|
forrest@0
|
134 |
"nwr_01D0_seas.co2", \
|
|
forrest@0
|
135 |
"obs023_06C3_seas.co2", \
|
|
forrest@0
|
136 |
"opw_01D0_seas.co2", \
|
|
forrest@0
|
137 |
"orl005_11D2_seas.co2", \
|
|
forrest@0
|
138 |
"orl015_11D2_seas.co2", \
|
|
forrest@0
|
139 |
"orl025_11D2_seas.co2", \
|
|
forrest@0
|
140 |
"orl035_11D2_seas.co2", \
|
|
forrest@0
|
141 |
"palcbc_30C0_seas.co2", \
|
|
forrest@0
|
142 |
"palmbc_30C0_seas.co2", \
|
|
forrest@0
|
143 |
"pdm_11D0_seas.co2", \
|
|
forrest@0
|
144 |
"pfa015_01D2_seas.co2", \
|
|
forrest@0
|
145 |
"pfa025_01D2_seas.co2", \
|
|
forrest@0
|
146 |
"pfa035_01D2_seas.co2", \
|
|
forrest@0
|
147 |
"pfa045_01D2_seas.co2", \
|
|
forrest@0
|
148 |
"pfa055_01D2_seas.co2", \
|
|
forrest@0
|
149 |
"pfa065_01D2_seas.co2", \
|
|
forrest@0
|
150 |
"pfa075_01D2_seas.co2", \
|
|
forrest@0
|
151 |
"poc000_01D1_seas.co2", \
|
|
forrest@0
|
152 |
"pocn05_01D1_seas.co2", \
|
|
forrest@0
|
153 |
"pocn10_01D1_seas.co2", \
|
|
forrest@0
|
154 |
"pocn15_01D1_seas.co2", \
|
|
forrest@0
|
155 |
"pocn20_01D1_seas.co2", \
|
|
forrest@0
|
156 |
"pocn25_01D1_seas.co2", \
|
|
forrest@0
|
157 |
"pocn30_01D1_seas.co2", \
|
|
forrest@0
|
158 |
"pocn35_01D1_seas.co2", \
|
|
forrest@0
|
159 |
"pocn40_01D1_seas.co2", \
|
|
forrest@0
|
160 |
"pocn45_01D1_seas.co2", \
|
|
forrest@0
|
161 |
"pocs05_01D1_seas.co2", \
|
|
forrest@0
|
162 |
"pocs10_01D1_seas.co2", \
|
|
forrest@0
|
163 |
"pocs15_01D1_seas.co2", \
|
|
forrest@0
|
164 |
"pocs20_01D1_seas.co2", \
|
|
forrest@0
|
165 |
"pocs25_01D1_seas.co2", \
|
|
forrest@0
|
166 |
"pocs30_01D1_seas.co2", \
|
|
forrest@0
|
167 |
"pocs35_01D1_seas.co2", \
|
|
forrest@0
|
168 |
"prs_21C0_seas.co2", \
|
|
forrest@0
|
169 |
"prs_21D0_seas.co2", \
|
|
forrest@0
|
170 |
"psa_01D0_seas.co2", \
|
|
forrest@0
|
171 |
"psa_04D0_seas.co2", \
|
|
forrest@0
|
172 |
"rpb_01D0_seas.co2", \
|
|
forrest@0
|
173 |
"rta005_01D2_seas.co2", \
|
|
forrest@0
|
174 |
"rta015_01D2_seas.co2", \
|
|
forrest@0
|
175 |
"rta025_01D2_seas.co2", \
|
|
forrest@0
|
176 |
"rta035_01D2_seas.co2", \
|
|
forrest@0
|
177 |
"rta045_01D2_seas.co2", \
|
|
forrest@0
|
178 |
"ryo_19C0_seas.co2", \
|
|
forrest@0
|
179 |
"sbl_06C0_seas.co2", \
|
|
forrest@0
|
180 |
"sbl_06D0_seas.co2", \
|
|
forrest@0
|
181 |
"sch_23C0_seas.co2", \
|
|
forrest@0
|
182 |
"scsn03_01D1_seas.co2", \
|
|
forrest@0
|
183 |
"scsn06_01D1_seas.co2", \
|
|
forrest@0
|
184 |
"scsn09_01D1_seas.co2", \
|
|
forrest@0
|
185 |
"scsn12_01D1_seas.co2", \
|
|
forrest@0
|
186 |
"scsn15_01D1_seas.co2", \
|
|
forrest@0
|
187 |
"scsn18_01D1_seas.co2", \
|
|
forrest@0
|
188 |
"scsn21_01D1_seas.co2", \
|
|
forrest@0
|
189 |
"sey_01D0_seas.co2", \
|
|
forrest@0
|
190 |
"shm_01D0_seas.co2", \
|
|
forrest@0
|
191 |
"sis_02D0_seas.co2", \
|
|
forrest@0
|
192 |
"smo_01C0_seas.co2", \
|
|
forrest@0
|
193 |
"smo_01D0_seas.co2", \
|
|
forrest@0
|
194 |
"smo_04D0_seas.co2", \
|
|
forrest@0
|
195 |
"spo_01C0_seas.co2", \
|
|
forrest@0
|
196 |
"spo_01D0_seas.co2", \
|
|
forrest@0
|
197 |
"spo_02D0_seas.co2", \
|
|
forrest@0
|
198 |
"spo_04D0_seas.co2", \
|
|
forrest@0
|
199 |
"stm_01D0_seas.co2", \
|
|
forrest@0
|
200 |
"stmebc_01D0_seas.co2", \
|
|
forrest@0
|
201 |
"stp_12D0_seas.co2", \
|
|
forrest@0
|
202 |
"sum_01D0_seas.co2", \
|
|
forrest@0
|
203 |
"syo_01D0_seas.co2", \
|
|
forrest@0
|
204 |
"syo_09C0_seas.co2", \
|
|
forrest@0
|
205 |
"tap_01D0_seas.co2", \
|
|
forrest@0
|
206 |
"tdf_01D0_seas.co2", \
|
|
forrest@0
|
207 |
"trm_11D0_seas.co2", \
|
|
forrest@0
|
208 |
"uta_01D0_seas.co2", \
|
|
forrest@0
|
209 |
"uum_01D0_seas.co2", \
|
|
forrest@0
|
210 |
"wes_23C0_seas.co2", \
|
|
forrest@0
|
211 |
"wis_01D0_seas.co2", \
|
|
forrest@0
|
212 |
"wkt009_01C3_seas.co2", \
|
|
forrest@0
|
213 |
"wkt030_01C3_seas.co2", \
|
|
forrest@0
|
214 |
"wkt061_01C3_seas.co2", \
|
|
forrest@0
|
215 |
"wkt122_01C3_seas.co2", \
|
|
forrest@0
|
216 |
"wkt244_01C3_seas.co2", \
|
|
forrest@0
|
217 |
"wkt457_01C3_seas.co2", \
|
|
forrest@0
|
218 |
"wlg_01D0_seas.co2", \
|
|
forrest@0
|
219 |
"wlg_33C0_seas.co2", \
|
|
forrest@0
|
220 |
"wpo000_10D2_seas.co2", \
|
|
forrest@0
|
221 |
"wpon05_10D2_seas.co2", \
|
|
forrest@0
|
222 |
"wpon10_10D2_seas.co2", \
|
|
forrest@0
|
223 |
"wpon15_10D2_seas.co2", \
|
|
forrest@0
|
224 |
"wpon20_10D2_seas.co2", \
|
|
forrest@0
|
225 |
"wpon25_10D2_seas.co2", \
|
|
forrest@0
|
226 |
"wpon30_10D2_seas.co2", \
|
|
forrest@0
|
227 |
"wpos05_10D2_seas.co2", \
|
|
forrest@0
|
228 |
"wpos10_10D2_seas.co2", \
|
|
forrest@0
|
229 |
"wpos15_10D2_seas.co2", \
|
|
forrest@0
|
230 |
"wpos20_10D2_seas.co2", \
|
|
forrest@0
|
231 |
"wpos25_10D2_seas.co2", \
|
|
forrest@0
|
232 |
"yon_19C0_seas.co2", \
|
|
forrest@0
|
233 |
"zep_01D0_seas.co2", \
|
|
forrest@0
|
234 |
"zep_31C0_seas.co2", \
|
|
forrest@0
|
235 |
"zot005_45D2_seas.co2", \
|
|
forrest@0
|
236 |
"zot015_45D2_seas.co2", \
|
|
forrest@0
|
237 |
"zot025_45D2_seas.co2", \
|
|
forrest@0
|
238 |
"zot035_45D2_seas.co2"/)
|
|
forrest@0
|
239 |
|
|
forrest@0
|
240 |
ntime = 12
|
|
forrest@0
|
241 |
ncase = dimsizes(b)
|
|
forrest@0
|
242 |
print (ncase)
|
|
forrest@0
|
243 |
|
|
forrest@0
|
244 |
mean = new((/ncase,ntime/),float)
|
|
forrest@0
|
245 |
lat = new((/ncase/),float)
|
|
forrest@0
|
246 |
lon = new((/ncase/),float)
|
|
forrest@0
|
247 |
elev = new((/ncase/),float)
|
|
forrest@0
|
248 |
|
|
forrest@0
|
249 |
do n = 0,dimsizes(b)-1
|
|
forrest@0
|
250 |
fili = diri + b(n)
|
|
forrest@0
|
251 |
print (fili)
|
|
forrest@0
|
252 |
|
|
forrest@0
|
253 |
ncol1 = 4
|
|
forrest@0
|
254 |
data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
|
|
forrest@0
|
255 |
|
|
forrest@0
|
256 |
lat(n) = data1(0,0)
|
|
forrest@0
|
257 |
lon(n) = data1(0,1)
|
|
forrest@0
|
258 |
elev(n) = data1(0,2)
|
|
forrest@0
|
259 |
|
|
forrest@0
|
260 |
ncol2 = 5
|
|
forrest@0
|
261 |
data2 = readAsciiTable(fili, ncol2, "float", 16)
|
|
forrest@0
|
262 |
|
|
forrest@0
|
263 |
mean(n,:) = data2(:,1)
|
|
forrest@0
|
264 |
end do
|
|
forrest@0
|
265 |
|
|
forrest@0
|
266 |
print (lat)
|
|
forrest@0
|
267 |
print (mean)
|
|
forrest@0
|
268 |
|
|
forrest@0
|
269 |
|
|
forrest@0
|
270 |
end
|