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")
15 diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
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"/)
244 mean = new((/ncase,ntime/),float)
245 lat = new((/ncase/),float)
246 lon = new((/ncase/),float)
247 elev = new((/ncase/),float)
249 do n = 0,dimsizes(b)-1
254 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
261 data2 = readAsciiTable(fili, ncol2, "float", 16)
263 mean(n,:) = data2(:,1)