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