1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/co2/05.read_ascci_122.ncl Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,195 @@
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_.nc"
1.15 + c = addfile(diro+filo,"c")
1.16 + filedimdef(c,"time",-1,True)
1.17 +
1.18 +; input ascii files
1.19 + diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
1.20 + b = (/ \
1.21 +"aia005_02D2_seas.co2", \
1.22 +"alt_01D0_seas.co2", \
1.23 +"alt_06C0_seas.co2", \
1.24 +"alt_06D0_seas.co2", \
1.25 +"ams_01D0_seas.co2", \
1.26 +"ams_11C0_seas.co2", \
1.27 +"asc_01D0_seas.co2", \
1.28 +"ask_01D0_seas.co2", \
1.29 +"avi_01D0_seas.co2", \
1.30 +"azr_01D0_seas.co2", \
1.31 +"bal_01D1_seas.co2", \
1.32 +"bgu_11D0_seas.co2", \
1.33 +"bhd_15C0_seas.co2", \
1.34 +"bme_01D0_seas.co2", \
1.35 +"bmw_01D0_seas.co2", \
1.36 +"brw_01C0_seas.co2", \
1.37 +"brw_01D0_seas.co2", \
1.38 +"bsc_01D0_seas.co2", \
1.39 +"car030_01D2_seas.co2", \
1.40 +"cba_01D0_seas.co2", \
1.41 +"cba_04D0_seas.co2", \
1.42 +"cfa_02D0_seas.co2", \
1.43 +"cgo_01D0_seas.co2", \
1.44 +"chr_01D0_seas.co2", \
1.45 +"cmn_17C0_seas.co2", \
1.46 +"cmo_01D0_seas.co2", \
1.47 +"coi_20C0_seas.co2", \
1.48 +"cpt_36C0_seas.co2", \
1.49 +"cri_02D0_seas.co2", \
1.50 +"crz_01D0_seas.co2", \
1.51 +"csj_06D0_seas.co2", \
1.52 +"daa_02D0_seas.co2", \
1.53 +"eic_01D0_seas.co2", \
1.54 +"esp005_01D2_seas.co2", \
1.55 +"esp_02D0_seas.co2", \
1.56 +"frd040_06C3_seas.co2", \
1.57 +"gmi_01D0_seas.co2", \
1.58 +"goz_01D0_seas.co2", \
1.59 +"gsn_24D0_seas.co2", \
1.60 +"haa005_01D2_seas.co2", \
1.61 +"hat_20C0_seas.co2", \
1.62 +"hba_01D0_seas.co2", \
1.63 +"hfm005_01D2_seas.co2", \
1.64 +"hun010_35C3_seas.co2", \
1.65 +"hun_01D0_seas.co2", \
1.66 +"ice_01D0_seas.co2", \
1.67 +"itn051_01C3_seas.co2", \
1.68 +"itn123_01C3_seas.co2", \
1.69 +"itn496_01C3_seas.co2", \
1.70 +"itn_01D0_seas.co2", \
1.71 +"izo_01D0_seas.co2", \
1.72 +"izo_27C0_seas.co2", \
1.73 +"jbn_29C0_seas.co2", \
1.74 +"key_01D0_seas.co2", \
1.75 +"kum_01D0_seas.co2", \
1.76 +"kzd_01D0_seas.co2", \
1.77 +"kzm_01D0_seas.co2", \
1.78 +"lef011_01C3_seas.co2", \
1.79 +"lef_01D0_seas.co2", \
1.80 +"ljo_04D0_seas.co2", \
1.81 +"lmp_28D0_seas.co2", \
1.82 +"maa_02D0_seas.co2", \
1.83 +"mbc_01D0_seas.co2", \
1.84 +"mhd_01D0_seas.co2", \
1.85 +"mhdcbc_11C0_seas.co2", \
1.86 +"mhdrbc_11C0_seas.co2", \
1.87 +"mid_01D0_seas.co2", \
1.88 +"mlo_01C0_seas.co2", \
1.89 +"mlo_01D0_seas.co2", \
1.90 +"mnm_19C0_seas.co2", \
1.91 +"mqa_02D0_seas.co2", \
1.92 +"nwr_01D0_seas.co2", \
1.93 +"obs023_06C3_seas.co2", \
1.94 +"opw_01D0_seas.co2", \
1.95 +"orl005_11D2_seas.co2", \
1.96 +"palcbc_30C0_seas.co2", \
1.97 +"palmbc_30C0_seas.co2", \
1.98 +"pdm_11D0_seas.co2", \
1.99 +"pfa015_01D2_seas.co2", \
1.100 +"poc000_01D1_seas.co2", \
1.101 +"pocn05_01D1_seas.co2", \
1.102 +"pocs05_01D1_seas.co2", \
1.103 +"prs_21C0_seas.co2", \
1.104 +"prs_21D0_seas.co2", \
1.105 +"psa_01D0_seas.co2", \
1.106 +"rpb_01D0_seas.co2", \
1.107 +"rta005_01D2_seas.co2", \
1.108 +"ryo_19C0_seas.co2", \
1.109 +"sbl_06C0_seas.co2", \
1.110 +"sbl_06D0_seas.co2", \
1.111 +"sch_23C0_seas.co2", \
1.112 +"scsn03_01D1_seas.co2", \
1.113 +"sey_01D0_seas.co2", \
1.114 +"shm_01D0_seas.co2", \
1.115 +"sis_02D0_seas.co2", \
1.116 +"smo_01C0_seas.co2", \
1.117 +"smo_01D0_seas.co2", \
1.118 +"spo_01C0_seas.co2", \
1.119 +"spo_01D0_seas.co2", \
1.120 +"stm_01D0_seas.co2", \
1.121 +"stmebc_01D0_seas.co2", \
1.122 +"stp_12D0_seas.co2", \
1.123 +"sum_01D0_seas.co2", \
1.124 +"syo_01D0_seas.co2", \
1.125 +"syo_09C0_seas.co2", \
1.126 +"tap_01D0_seas.co2", \
1.127 +"tdf_01D0_seas.co2", \
1.128 +"trm_11D0_seas.co2", \
1.129 +"uta_01D0_seas.co2", \
1.130 +"uum_01D0_seas.co2", \
1.131 +"wes_23C0_seas.co2", \
1.132 +"wis_01D0_seas.co2", \
1.133 +"wkt009_01C3_seas.co2", \
1.134 +"wlg_01D0_seas.co2", \
1.135 +"wlg_33C0_seas.co2", \
1.136 +"wpo000_10D2_seas.co2", \
1.137 +"wpon05_10D2_seas.co2", \
1.138 +"wpos05_10D2_seas.co2", \
1.139 +"yon_19C0_seas.co2", \
1.140 +"zep_01D0_seas.co2", \
1.141 +"zep_31C0_seas.co2", \
1.142 +"zot005_45D2_seas.co2"/)
1.143 +
1.144 + ntime = 12
1.145 + ncase = dimsizes(b)
1.146 + print (ncase)
1.147 +
1.148 + mean = new((/ncase,ntime/),float)
1.149 + lat = new((/ncase/),float)
1.150 + lon = new((/ncase/),float)
1.151 + elev = new((/ncase/),float)
1.152 +
1.153 + do n = 0,dimsizes(b)-1
1.154 + fili = diri + b(n)
1.155 + print (fili)
1.156 +
1.157 + ncol1 = 4
1.158 + data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
1.159 +
1.160 + lat(n) = data1(0,0)
1.161 + lon(n) = data1(0,1)
1.162 + elev(n) = data1(0,2)
1.163 +
1.164 + ncol2 = 5
1.165 + data2 = readAsciiTable(fili, ncol2, "float", 16)
1.166 +
1.167 + mean(n,:) = data2(:,1)
1.168 + end do
1.169 +
1.170 +; mean!0 = "case"
1.171 +; mean!1 = "time"
1.172 +; mean&case= ispan(1,ncase,1)
1.173 +; mean&time= ispan(1,ntime,1)
1.174 + mean@units = "u mol/mol"
1.175 + mean@long_name = "CO2 GlobalView seasonal"
1.176 +
1.177 +; lat!0 = "case"
1.178 +; lon!0 = "case"
1.179 +; elev!0 = "case"
1.180 +; lat&case = mean&case
1.181 +; lon&case = mean&case
1.182 +; elev&case = mean&case
1.183 +
1.184 + lat@units = "degrees_north"
1.185 + lat@long_name = "Latitude"
1.186 + lon@units = "degrees_ease"
1.187 + lon@long_name = "Longitude"
1.188 + elev@units = "m"
1.189 + elev@long_name = "Elevation above sea level"
1.190 +
1.191 + print (lat)
1.192 + print (mean)
1.193 +
1.194 + c->LAT = lat
1.195 + c->LON = lon
1.196 + c->ELEV = elev
1.197 + c->CO2_SEAS = mean
1.198 +end