1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/co2/04.test.read.ncl Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,270 @@
1.4 +;************************************************
1.5 +; Read ascii, Write nc
1.6 +;************************************************
1.7 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
1.8 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
1.9 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
1.10 +;************************************************
1.11 +begin
1.12 +; final output
1.13 +; diro = "/fis/cgd/cseg/people/jeff/clamp_data/co2/"
1.14 +; filo = "co2_globalView_222.nc"
1.15 +; c = addfile(diro+filo,"c")
1.16 +
1.17 +; input ascii files
1.18 + diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
1.19 + b = (/ \
1.20 +"aia005_02D2_seas.co2", \
1.21 +"aia015_02D2_seas.co2", \
1.22 +"aia025_02D2_seas.co2", \
1.23 +"aia035_02D2_seas.co2", \
1.24 +"aia045_02D2_seas.co2", \
1.25 +"aia055_02D2_seas.co2", \
1.26 +"aia065_02D2_seas.co2", \
1.27 +"alt_01D0_seas.co2", \
1.28 +"alt_02D0_seas.co2", \
1.29 +"alt_04D0_seas.co2", \
1.30 +"alt_06C0_seas.co2", \
1.31 +"alt_06D0_seas.co2", \
1.32 +"ams_01D0_seas.co2", \
1.33 +"ams_11C0_seas.co2", \
1.34 +"asc_01D0_seas.co2", \
1.35 +"ask_01D0_seas.co2", \
1.36 +"avi_01D0_seas.co2", \
1.37 +"azr_01D0_seas.co2", \
1.38 +"bal_01D1_seas.co2", \
1.39 +"bgu_11D0_seas.co2", \
1.40 +"bhd_15C0_seas.co2", \
1.41 +"bme_01D0_seas.co2", \
1.42 +"bmw_01D0_seas.co2", \
1.43 +"brw_01C0_seas.co2", \
1.44 +"brw_01D0_seas.co2", \
1.45 +"bsc_01D0_seas.co2", \
1.46 +"car030_01D2_seas.co2", \
1.47 +"car040_01D2_seas.co2", \
1.48 +"car050_01D2_seas.co2", \
1.49 +"car060_01D2_seas.co2", \
1.50 +"car070_01D2_seas.co2", \
1.51 +"car080_01D2_seas.co2", \
1.52 +"cba_01D0_seas.co2", \
1.53 +"cba_04D0_seas.co2", \
1.54 +"cfa_02D0_seas.co2", \
1.55 +"cgo_01D0_seas.co2", \
1.56 +"cgo_02D0_seas.co2", \
1.57 +"cgo_04D0_seas.co2", \
1.58 +"chr_01D0_seas.co2", \
1.59 +"cmn_17C0_seas.co2", \
1.60 +"cmo_01D0_seas.co2", \
1.61 +"coi_20C0_seas.co2", \
1.62 +"cpt_36C0_seas.co2", \
1.63 +"cri_02D0_seas.co2", \
1.64 +"crz_01D0_seas.co2", \
1.65 +"csj_06D0_seas.co2", \
1.66 +"daa_02D0_seas.co2", \
1.67 +"eic_01D0_seas.co2", \
1.68 +"esp005_01D2_seas.co2", \
1.69 +"esp015_01D2_seas.co2", \
1.70 +"esp025_01D2_seas.co2", \
1.71 +"esp035_01D2_seas.co2", \
1.72 +"esp045_01D2_seas.co2", \
1.73 +"esp055_01D2_seas.co2", \
1.74 +"esp_02D0_seas.co2", \
1.75 +"esp_06D0_seas.co2", \
1.76 +"frd040_06C3_seas.co2", \
1.77 +"gmi_01D0_seas.co2", \
1.78 +"goz_01D0_seas.co2", \
1.79 +"gsn_24D0_seas.co2", \
1.80 +"haa005_01D2_seas.co2", \
1.81 +"haa015_01D2_seas.co2", \
1.82 +"haa025_01D2_seas.co2", \
1.83 +"haa035_01D2_seas.co2", \
1.84 +"haa045_01D2_seas.co2", \
1.85 +"haa055_01D2_seas.co2", \
1.86 +"haa065_01D2_seas.co2", \
1.87 +"haa075_01D2_seas.co2", \
1.88 +"hat_20C0_seas.co2", \
1.89 +"hba_01D0_seas.co2", \
1.90 +"hfm005_01D2_seas.co2", \
1.91 +"hfm015_01D2_seas.co2", \
1.92 +"hfm025_01D2_seas.co2", \
1.93 +"hfm035_01D2_seas.co2", \
1.94 +"hfm045_01D2_seas.co2", \
1.95 +"hfm055_01D2_seas.co2", \
1.96 +"hfm065_01D2_seas.co2", \
1.97 +"hfm075_01D2_seas.co2", \
1.98 +"hun010_35C3_seas.co2", \
1.99 +"hun048_35C3_seas.co2", \
1.100 +"hun082_35C3_seas.co2", \
1.101 +"hun115_35C3_seas.co2", \
1.102 +"hun_01D0_seas.co2", \
1.103 +"ice_01D0_seas.co2", \
1.104 +"itn051_01C3_seas.co2", \
1.105 +"itn123_01C3_seas.co2", \
1.106 +"itn496_01C3_seas.co2", \
1.107 +"itn_01D0_seas.co2", \
1.108 +"izo_01D0_seas.co2", \
1.109 +"izo_27C0_seas.co2", \
1.110 +"jbn_29C0_seas.co2", \
1.111 +"key_01D0_seas.co2", \
1.112 +"kum_01D0_seas.co2", \
1.113 +"kum_04D0_seas.co2", \
1.114 +"kzd_01D0_seas.co2", \
1.115 +"kzm_01D0_seas.co2", \
1.116 +"lef011_01C3_seas.co2", \
1.117 +"lef030_01C3_seas.co2", \
1.118 +"lef076_01C3_seas.co2", \
1.119 +"lef122_01C3_seas.co2", \
1.120 +"lef244_01C3_seas.co2", \
1.121 +"lef396_01C3_seas.co2", \
1.122 +"lef_01D0_seas.co2", \
1.123 +"ljo_04D0_seas.co2", \
1.124 +"lmp_28D0_seas.co2", \
1.125 +"maa_02D0_seas.co2", \
1.126 +"mbc_01D0_seas.co2", \
1.127 +"mhd_01D0_seas.co2", \
1.128 +"mhdcbc_11C0_seas.co2", \
1.129 +"mhdrbc_11C0_seas.co2", \
1.130 +"mid_01D0_seas.co2", \
1.131 +"mlo_01C0_seas.co2", \
1.132 +"mlo_01D0_seas.co2", \
1.133 +"mlo_02D0_seas.co2", \
1.134 +"mlo_04D0_seas.co2", \
1.135 +"mnm_19C0_seas.co2", \
1.136 +"mqa_02D0_seas.co2", \
1.137 +"nwr_01D0_seas.co2", \
1.138 +"obs023_06C3_seas.co2", \
1.139 +"opw_01D0_seas.co2", \
1.140 +"orl005_11D2_seas.co2", \
1.141 +"orl015_11D2_seas.co2", \
1.142 +"orl025_11D2_seas.co2", \
1.143 +"orl035_11D2_seas.co2", \
1.144 +"palcbc_30C0_seas.co2", \
1.145 +"palmbc_30C0_seas.co2", \
1.146 +"pdm_11D0_seas.co2", \
1.147 +"pfa015_01D2_seas.co2", \
1.148 +"pfa025_01D2_seas.co2", \
1.149 +"pfa035_01D2_seas.co2", \
1.150 +"pfa045_01D2_seas.co2", \
1.151 +"pfa055_01D2_seas.co2", \
1.152 +"pfa065_01D2_seas.co2", \
1.153 +"pfa075_01D2_seas.co2", \
1.154 +"poc000_01D1_seas.co2", \
1.155 +"pocn05_01D1_seas.co2", \
1.156 +"pocn10_01D1_seas.co2", \
1.157 +"pocn15_01D1_seas.co2", \
1.158 +"pocn20_01D1_seas.co2", \
1.159 +"pocn25_01D1_seas.co2", \
1.160 +"pocn30_01D1_seas.co2", \
1.161 +"pocn35_01D1_seas.co2", \
1.162 +"pocn40_01D1_seas.co2", \
1.163 +"pocn45_01D1_seas.co2", \
1.164 +"pocs05_01D1_seas.co2", \
1.165 +"pocs10_01D1_seas.co2", \
1.166 +"pocs15_01D1_seas.co2", \
1.167 +"pocs20_01D1_seas.co2", \
1.168 +"pocs25_01D1_seas.co2", \
1.169 +"pocs30_01D1_seas.co2", \
1.170 +"pocs35_01D1_seas.co2", \
1.171 +"prs_21C0_seas.co2", \
1.172 +"prs_21D0_seas.co2", \
1.173 +"psa_01D0_seas.co2", \
1.174 +"psa_04D0_seas.co2", \
1.175 +"rpb_01D0_seas.co2", \
1.176 +"rta005_01D2_seas.co2", \
1.177 +"rta015_01D2_seas.co2", \
1.178 +"rta025_01D2_seas.co2", \
1.179 +"rta035_01D2_seas.co2", \
1.180 +"rta045_01D2_seas.co2", \
1.181 +"ryo_19C0_seas.co2", \
1.182 +"sbl_06C0_seas.co2", \
1.183 +"sbl_06D0_seas.co2", \
1.184 +"sch_23C0_seas.co2", \
1.185 +"scsn03_01D1_seas.co2", \
1.186 +"scsn06_01D1_seas.co2", \
1.187 +"scsn09_01D1_seas.co2", \
1.188 +"scsn12_01D1_seas.co2", \
1.189 +"scsn15_01D1_seas.co2", \
1.190 +"scsn18_01D1_seas.co2", \
1.191 +"scsn21_01D1_seas.co2", \
1.192 +"sey_01D0_seas.co2", \
1.193 +"shm_01D0_seas.co2", \
1.194 +"sis_02D0_seas.co2", \
1.195 +"smo_01C0_seas.co2", \
1.196 +"smo_01D0_seas.co2", \
1.197 +"smo_04D0_seas.co2", \
1.198 +"spo_01C0_seas.co2", \
1.199 +"spo_01D0_seas.co2", \
1.200 +"spo_02D0_seas.co2", \
1.201 +"spo_04D0_seas.co2", \
1.202 +"stm_01D0_seas.co2", \
1.203 +"stmebc_01D0_seas.co2", \
1.204 +"stp_12D0_seas.co2", \
1.205 +"sum_01D0_seas.co2", \
1.206 +"syo_01D0_seas.co2", \
1.207 +"syo_09C0_seas.co2", \
1.208 +"tap_01D0_seas.co2", \
1.209 +"tdf_01D0_seas.co2", \
1.210 +"trm_11D0_seas.co2", \
1.211 +"uta_01D0_seas.co2", \
1.212 +"uum_01D0_seas.co2", \
1.213 +"wes_23C0_seas.co2", \
1.214 +"wis_01D0_seas.co2", \
1.215 +"wkt009_01C3_seas.co2", \
1.216 +"wkt030_01C3_seas.co2", \
1.217 +"wkt061_01C3_seas.co2", \
1.218 +"wkt122_01C3_seas.co2", \
1.219 +"wkt244_01C3_seas.co2", \
1.220 +"wkt457_01C3_seas.co2", \
1.221 +"wlg_01D0_seas.co2", \
1.222 +"wlg_33C0_seas.co2", \
1.223 +"wpo000_10D2_seas.co2", \
1.224 +"wpon05_10D2_seas.co2", \
1.225 +"wpon10_10D2_seas.co2", \
1.226 +"wpon15_10D2_seas.co2", \
1.227 +"wpon20_10D2_seas.co2", \
1.228 +"wpon25_10D2_seas.co2", \
1.229 +"wpon30_10D2_seas.co2", \
1.230 +"wpos05_10D2_seas.co2", \
1.231 +"wpos10_10D2_seas.co2", \
1.232 +"wpos15_10D2_seas.co2", \
1.233 +"wpos20_10D2_seas.co2", \
1.234 +"wpos25_10D2_seas.co2", \
1.235 +"yon_19C0_seas.co2", \
1.236 +"zep_01D0_seas.co2", \
1.237 +"zep_31C0_seas.co2", \
1.238 +"zot005_45D2_seas.co2", \
1.239 +"zot015_45D2_seas.co2", \
1.240 +"zot025_45D2_seas.co2", \
1.241 +"zot035_45D2_seas.co2"/)
1.242 +
1.243 + ntime = 12
1.244 + ncase = dimsizes(b)
1.245 + print (ncase)
1.246 +
1.247 + mean = new((/ncase,ntime/),float)
1.248 + lat = new((/ncase/),float)
1.249 + lon = new((/ncase/),float)
1.250 + elev = new((/ncase/),float)
1.251 +
1.252 + do n = 0,dimsizes(b)-1
1.253 + fili = diri + b(n)
1.254 + print (fili)
1.255 +
1.256 + ncol1 = 4
1.257 + data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
1.258 +
1.259 + lat(n) = data1(0,0)
1.260 + lon(n) = data1(0,1)
1.261 + elev(n) = data1(0,2)
1.262 +
1.263 + ncol2 = 5
1.264 + data2 = readAsciiTable(fili, ncol2, "float", 16)
1.265 +
1.266 + mean(n,:) = data2(:,1)
1.267 + end do
1.268 +
1.269 + print (lat)
1.270 + print (mean)
1.271 +
1.272 +
1.273 +end