author | Forrest Hoffman <forrest@climatemodeling.org> |
Mon, 26 Jan 2009 22:08:20 -0500 | |
changeset 0 | 0c6405ab2ff4 |
permissions | -rw-r--r-- |
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)
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"/)
117 ntime = 12
118 ncase = dimsizes(b)
119 print (ncase)
121 mean = new((/ncase,ntime/),float)
122 lat = new((/ncase/),float)
123 lon = new((/ncase/),float)
124 elev = new((/ncase/),float)
126 do n = 0,dimsizes(b)-1
127 fili = diri + b(n) + "_seas.co2"
128 print (fili)
130 ncol1 = 4
131 data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
133 lat(n) = data1(0,0)
134 lon(n) = data1(0,1)
135 elev(n) = data1(0,2)
137 ncol2 = 5
138 data2 = readAsciiTable(fili, ncol2, "float", 16)
140 mean(n,:) = data2(:,1)
141 end do
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"
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
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"
166 print (lat)
167 print (mean)
168 print (b)
170 c->LAT = lat
171 c->LON = lon
172 c->ELEV = elev
173 c->CO2_SEAS = mean
174 c->STATION = stringtochar(b)
175 end