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