|
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_.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 "alt_01D0_seas.co2", \ |
|
20 "alt_06C0_seas.co2", \ |
|
21 "alt_06D0_seas.co2", \ |
|
22 "ams_01D0_seas.co2", \ |
|
23 "ams_11C0_seas.co2", \ |
|
24 "asc_01D0_seas.co2", \ |
|
25 "ask_01D0_seas.co2", \ |
|
26 "avi_01D0_seas.co2", \ |
|
27 "azr_01D0_seas.co2", \ |
|
28 "bal_01D1_seas.co2", \ |
|
29 "bgu_11D0_seas.co2", \ |
|
30 "bhd_15C0_seas.co2", \ |
|
31 "bme_01D0_seas.co2", \ |
|
32 "bmw_01D0_seas.co2", \ |
|
33 "brw_01C0_seas.co2", \ |
|
34 "brw_01D0_seas.co2", \ |
|
35 "bsc_01D0_seas.co2", \ |
|
36 "car030_01D2_seas.co2", \ |
|
37 "cba_01D0_seas.co2", \ |
|
38 "cba_04D0_seas.co2", \ |
|
39 "cfa_02D0_seas.co2", \ |
|
40 "cgo_01D0_seas.co2", \ |
|
41 "chr_01D0_seas.co2", \ |
|
42 "cmn_17C0_seas.co2", \ |
|
43 "cmo_01D0_seas.co2", \ |
|
44 "coi_20C0_seas.co2", \ |
|
45 "cpt_36C0_seas.co2", \ |
|
46 "cri_02D0_seas.co2", \ |
|
47 "crz_01D0_seas.co2", \ |
|
48 "csj_06D0_seas.co2", \ |
|
49 "daa_02D0_seas.co2", \ |
|
50 "eic_01D0_seas.co2", \ |
|
51 "esp005_01D2_seas.co2", \ |
|
52 "esp_02D0_seas.co2", \ |
|
53 "frd040_06C3_seas.co2", \ |
|
54 "gmi_01D0_seas.co2", \ |
|
55 "goz_01D0_seas.co2", \ |
|
56 "gsn_24D0_seas.co2", \ |
|
57 "haa005_01D2_seas.co2", \ |
|
58 "hat_20C0_seas.co2", \ |
|
59 "hba_01D0_seas.co2", \ |
|
60 "hfm005_01D2_seas.co2", \ |
|
61 "hun010_35C3_seas.co2", \ |
|
62 "hun_01D0_seas.co2", \ |
|
63 "ice_01D0_seas.co2", \ |
|
64 "itn051_01C3_seas.co2", \ |
|
65 "itn123_01C3_seas.co2", \ |
|
66 "itn496_01C3_seas.co2", \ |
|
67 "itn_01D0_seas.co2", \ |
|
68 "izo_01D0_seas.co2", \ |
|
69 "izo_27C0_seas.co2", \ |
|
70 "jbn_29C0_seas.co2", \ |
|
71 "key_01D0_seas.co2", \ |
|
72 "kum_01D0_seas.co2", \ |
|
73 "kzd_01D0_seas.co2", \ |
|
74 "kzm_01D0_seas.co2", \ |
|
75 "lef011_01C3_seas.co2", \ |
|
76 "lef_01D0_seas.co2", \ |
|
77 "ljo_04D0_seas.co2", \ |
|
78 "lmp_28D0_seas.co2", \ |
|
79 "maa_02D0_seas.co2", \ |
|
80 "mbc_01D0_seas.co2", \ |
|
81 "mhd_01D0_seas.co2", \ |
|
82 "mhdcbc_11C0_seas.co2", \ |
|
83 "mhdrbc_11C0_seas.co2", \ |
|
84 "mid_01D0_seas.co2", \ |
|
85 "mlo_01C0_seas.co2", \ |
|
86 "mlo_01D0_seas.co2", \ |
|
87 "mnm_19C0_seas.co2", \ |
|
88 "mqa_02D0_seas.co2", \ |
|
89 "nwr_01D0_seas.co2", \ |
|
90 "obs023_06C3_seas.co2", \ |
|
91 "opw_01D0_seas.co2", \ |
|
92 "orl005_11D2_seas.co2", \ |
|
93 "palcbc_30C0_seas.co2", \ |
|
94 "palmbc_30C0_seas.co2", \ |
|
95 "pdm_11D0_seas.co2", \ |
|
96 "pfa015_01D2_seas.co2", \ |
|
97 "poc000_01D1_seas.co2", \ |
|
98 "pocn05_01D1_seas.co2", \ |
|
99 "pocs05_01D1_seas.co2", \ |
|
100 "prs_21C0_seas.co2", \ |
|
101 "prs_21D0_seas.co2", \ |
|
102 "psa_01D0_seas.co2", \ |
|
103 "rpb_01D0_seas.co2", \ |
|
104 "rta005_01D2_seas.co2", \ |
|
105 "ryo_19C0_seas.co2", \ |
|
106 "sbl_06C0_seas.co2", \ |
|
107 "sbl_06D0_seas.co2", \ |
|
108 "sch_23C0_seas.co2", \ |
|
109 "scsn03_01D1_seas.co2", \ |
|
110 "sey_01D0_seas.co2", \ |
|
111 "shm_01D0_seas.co2", \ |
|
112 "sis_02D0_seas.co2", \ |
|
113 "smo_01C0_seas.co2", \ |
|
114 "smo_01D0_seas.co2", \ |
|
115 "spo_01C0_seas.co2", \ |
|
116 "spo_01D0_seas.co2", \ |
|
117 "stm_01D0_seas.co2", \ |
|
118 "stmebc_01D0_seas.co2", \ |
|
119 "stp_12D0_seas.co2", \ |
|
120 "sum_01D0_seas.co2", \ |
|
121 "syo_01D0_seas.co2", \ |
|
122 "syo_09C0_seas.co2", \ |
|
123 "tap_01D0_seas.co2", \ |
|
124 "tdf_01D0_seas.co2", \ |
|
125 "trm_11D0_seas.co2", \ |
|
126 "uta_01D0_seas.co2", \ |
|
127 "uum_01D0_seas.co2", \ |
|
128 "wes_23C0_seas.co2", \ |
|
129 "wis_01D0_seas.co2", \ |
|
130 "wkt009_01C3_seas.co2", \ |
|
131 "wlg_01D0_seas.co2", \ |
|
132 "wlg_33C0_seas.co2", \ |
|
133 "wpo000_10D2_seas.co2", \ |
|
134 "wpon05_10D2_seas.co2", \ |
|
135 "wpos05_10D2_seas.co2", \ |
|
136 "yon_19C0_seas.co2", \ |
|
137 "zep_01D0_seas.co2", \ |
|
138 "zep_31C0_seas.co2", \ |
|
139 "zot005_45D2_seas.co2"/) |
|
140 |
|
141 ntime = 12 |
|
142 ncase = dimsizes(b) |
|
143 print (ncase) |
|
144 |
|
145 mean = new((/ncase,ntime/),float) |
|
146 lat = new((/ncase/),float) |
|
147 lon = new((/ncase/),float) |
|
148 elev = new((/ncase/),float) |
|
149 |
|
150 do n = 0,dimsizes(b)-1 |
|
151 fili = diri + b(n) |
|
152 print (fili) |
|
153 |
|
154 ncol1 = 4 |
|
155 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/)) |
|
156 |
|
157 lat(n) = data1(0,0) |
|
158 lon(n) = data1(0,1) |
|
159 elev(n) = data1(0,2) |
|
160 |
|
161 ncol2 = 5 |
|
162 data2 = readAsciiTable(fili, ncol2, "float", 16) |
|
163 |
|
164 mean(n,:) = data2(:,1) |
|
165 end do |
|
166 |
|
167 ; mean!0 = "case" |
|
168 ; mean!1 = "time" |
|
169 ; mean&case= ispan(1,ncase,1) |
|
170 ; mean&time= ispan(1,ntime,1) |
|
171 mean@units = "u mol/mol" |
|
172 mean@long_name = "CO2 GlobalView seasonal" |
|
173 |
|
174 ; lat!0 = "case" |
|
175 ; lon!0 = "case" |
|
176 ; elev!0 = "case" |
|
177 ; lat&case = mean&case |
|
178 ; lon&case = mean&case |
|
179 ; elev&case = mean&case |
|
180 |
|
181 lat@units = "degrees_north" |
|
182 lat@long_name = "Latitude" |
|
183 lon@units = "degrees_ease" |
|
184 lon@long_name = "Longitude" |
|
185 elev@units = "m" |
|
186 elev@long_name = "Elevation above sea level" |
|
187 |
|
188 print (lat) |
|
189 print (mean) |
|
190 |
|
191 c->LAT = lat |
|
192 c->LON = lon |
|
193 c->ELEV = elev |
|
194 c->CO2_SEAS = mean |
|
195 end |