Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;************************************************
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 ;************************************************
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)
16 diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
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"/)
245 mean = new((/ncase,ntime/),float)
246 lat = new((/ncase/),float)
247 lon = new((/ncase/),float)
248 elev = new((/ncase/),float)
250 do n = 0,dimsizes(b)-1
255 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
262 data2 = readAsciiTable(fili, ncol2, "float", 16)
264 mean(n,:) = data2(:,1)
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"
277 ; lat&case = mean&case
278 ; lon&case = mean&case
279 ; elev&case = mean&case
281 lat@units = "degrees_north"
282 lat@long_name = "Latitude"
283 lon@units = "degrees_ease"
284 lon@long_name = "Longitude"
286 elev@long_name = "Elevation above sea level"