|
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_98.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", \ |
|
19 "alt_01D0", \ |
|
20 "ams_01D0", \ |
|
21 "asc_01D0", \ |
|
22 "ask_01D0", \ |
|
23 "avi_01D0", \ |
|
24 "azr_01D0", \ |
|
25 "bal_01D1", \ |
|
26 "bgu_11D0", \ |
|
27 "bhd_15C0", \ |
|
28 "bme_01D0", \ |
|
29 "bmw_01D0", \ |
|
30 "brw_01C0", \ |
|
31 "bsc_01D0", \ |
|
32 "car030_01D2", \ |
|
33 "cba_01D0", \ |
|
34 "cfa_02D0", \ |
|
35 "cgo_01D0", \ |
|
36 "chr_01D0", \ |
|
37 "cmn_17C0", \ |
|
38 "cmo_01D0", \ |
|
39 "coi_20C0", \ |
|
40 "cpt_36C0", \ |
|
41 "cri_02D0", \ |
|
42 "crz_01D0", \ |
|
43 "csj_06D0", \ |
|
44 "daa_02D0", \ |
|
45 "eic_01D0", \ |
|
46 "esp_02D0", \ |
|
47 "frd040_06C3", \ |
|
48 "gmi_01D0", \ |
|
49 "goz_01D0", \ |
|
50 "gsn_24D0", \ |
|
51 "haa005_01D2", \ |
|
52 "hat_20C0", \ |
|
53 "hba_01D0", \ |
|
54 "hfm005_01D2", \ |
|
55 "hun010_35C3", \ |
|
56 "ice_01D0", \ |
|
57 "itn051_01C3", \ |
|
58 "izo_01D0", \ |
|
59 "jbn_29C0", \ |
|
60 "key_01D0", \ |
|
61 "kum_01D0", \ |
|
62 "kzd_01D0", \ |
|
63 "kzm_01D0", \ |
|
64 "lef011_01C3", \ |
|
65 "ljo_04D0", \ |
|
66 "lmp_28D0", \ |
|
67 "maa_02D0", \ |
|
68 "mbc_01D0", \ |
|
69 "mhd_01D0", \ |
|
70 "mid_01D0", \ |
|
71 "mlo_01C0", \ |
|
72 "mnm_19C0", \ |
|
73 "mqa_02D0", \ |
|
74 "nwr_01D0", \ |
|
75 "obs023_06C3", \ |
|
76 "opw_01D0", \ |
|
77 "orl005_11D2", \ |
|
78 "palcbc_30C0", \ |
|
79 "pdm_11D0", \ |
|
80 "pfa015_01D2", \ |
|
81 "poc000_01D1", \ |
|
82 "pocn05_01D1", \ |
|
83 "pocs05_01D1", \ |
|
84 "prs_21C0", \ |
|
85 "psa_01D0", \ |
|
86 "rpb_01D0", \ |
|
87 "rta005_01D2", \ |
|
88 "ryo_19C0", \ |
|
89 "sbl_06C0", \ |
|
90 "sch_23C0", \ |
|
91 "scsn03_01D1", \ |
|
92 "sey_01D0", \ |
|
93 "shm_01D0", \ |
|
94 "sis_02D0", \ |
|
95 "smo_01C0", \ |
|
96 "spo_01C0", \ |
|
97 "stm_01D0", \ |
|
98 "stp_12D0", \ |
|
99 "sum_01D0", \ |
|
100 "syo_01D0", \ |
|
101 "tap_01D0", \ |
|
102 "tdf_01D0", \ |
|
103 "trm_11D0", \ |
|
104 "uta_01D0", \ |
|
105 "uum_01D0", \ |
|
106 "wes_23C0", \ |
|
107 "wis_01D0", \ |
|
108 "wkt009_01C3", \ |
|
109 "wlg_01D0", \ |
|
110 "wpo000_10D2", \ |
|
111 "wpon05_10D2", \ |
|
112 "wpos05_10D2", \ |
|
113 "yon_19C0", \ |
|
114 "zep_01D0", \ |
|
115 "zot005_45D2"/) |
|
116 |
|
117 ntime = 12 |
|
118 ncase = dimsizes(b) |
|
119 print (ncase) |
|
120 |
|
121 mean = new((/ncase,ntime/),float) |
|
122 lat = new((/ncase/),float) |
|
123 lon = new((/ncase/),float) |
|
124 elev = new((/ncase/),float) |
|
125 |
|
126 do n = 0,dimsizes(b)-1 |
|
127 fili = diri + b(n) + "_seas.co2" |
|
128 print (fili) |
|
129 |
|
130 ncol1 = 4 |
|
131 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/)) |
|
132 |
|
133 lat(n) = data1(0,0) |
|
134 lon(n) = data1(0,1) |
|
135 elev(n) = data1(0,2) |
|
136 |
|
137 ncol2 = 5 |
|
138 data2 = readAsciiTable(fili, ncol2, "float", 16) |
|
139 |
|
140 mean(n,:) = data2(:,1) |
|
141 end do |
|
142 |
|
143 ; mean!0 = "case" |
|
144 ; mean!1 = "time" |
|
145 ; mean&case= ispan(1,ncase,1) |
|
146 ; mean&time= ispan(1,ntime,1) |
|
147 mean@units = "u mol/mol" |
|
148 mean@long_name = "CO2 GlobalView seasonal" |
|
149 |
|
150 ; lat!0 = "case" |
|
151 ; lon!0 = "case" |
|
152 ; elev!0 = "case" |
|
153 ; lat&case = mean&case |
|
154 ; lon&case = mean&case |
|
155 ; elev&case = mean&case |
|
156 |
|
157 lat@units = "degrees_north" |
|
158 lat@long_name = "Latitude" |
|
159 lon@units = "degrees_ease" |
|
160 lon@long_name = "Longitude" |
|
161 elev@units = "m" |
|
162 elev@long_name = "Elevation above sea level" |
|
163 b@units = "" |
|
164 b@long_name = "Station" |
|
165 |
|
166 print (lat) |
|
167 print (mean) |
|
168 print (b) |
|
169 |
|
170 c->LAT = lat |
|
171 c->LON = lon |
|
172 c->ELEV = elev |
|
173 c->CO2_SEAS = mean |
|
174 c->STATION = stringtochar(b) |
|
175 end |