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_.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 "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"/)
145 mean = new((/ncase,ntime/),float)
146 lat = new((/ncase/),float)
147 lon = new((/ncase/),float)
148 elev = new((/ncase/),float)
150 do n = 0,dimsizes(b)-1
155 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
162 data2 = readAsciiTable(fili, ncol2, "float", 16)
164 mean(n,:) = data2(:,1)
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"
177 ; lat&case = mean&case
178 ; lon&case = mean&case
179 ; elev&case = mean&case
181 lat@units = "degrees_north"
182 lat@long_name = "Latitude"
183 lon@units = "degrees_ease"
184 lon@long_name = "Longitude"
186 elev@long_name = "Elevation above sea level"