1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/co2/06.read_ascci_98.ncl Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,175 @@
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_98.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", \
1.22 +"alt_01D0", \
1.23 +"ams_01D0", \
1.24 +"asc_01D0", \
1.25 +"ask_01D0", \
1.26 +"avi_01D0", \
1.27 +"azr_01D0", \
1.28 +"bal_01D1", \
1.29 +"bgu_11D0", \
1.30 +"bhd_15C0", \
1.31 +"bme_01D0", \
1.32 +"bmw_01D0", \
1.33 +"brw_01C0", \
1.34 +"bsc_01D0", \
1.35 +"car030_01D2", \
1.36 +"cba_01D0", \
1.37 +"cfa_02D0", \
1.38 +"cgo_01D0", \
1.39 +"chr_01D0", \
1.40 +"cmn_17C0", \
1.41 +"cmo_01D0", \
1.42 +"coi_20C0", \
1.43 +"cpt_36C0", \
1.44 +"cri_02D0", \
1.45 +"crz_01D0", \
1.46 +"csj_06D0", \
1.47 +"daa_02D0", \
1.48 +"eic_01D0", \
1.49 +"esp_02D0", \
1.50 +"frd040_06C3", \
1.51 +"gmi_01D0", \
1.52 +"goz_01D0", \
1.53 +"gsn_24D0", \
1.54 +"haa005_01D2", \
1.55 +"hat_20C0", \
1.56 +"hba_01D0", \
1.57 +"hfm005_01D2", \
1.58 +"hun010_35C3", \
1.59 +"ice_01D0", \
1.60 +"itn051_01C3", \
1.61 +"izo_01D0", \
1.62 +"jbn_29C0", \
1.63 +"key_01D0", \
1.64 +"kum_01D0", \
1.65 +"kzd_01D0", \
1.66 +"kzm_01D0", \
1.67 +"lef011_01C3", \
1.68 +"ljo_04D0", \
1.69 +"lmp_28D0", \
1.70 +"maa_02D0", \
1.71 +"mbc_01D0", \
1.72 +"mhd_01D0", \
1.73 +"mid_01D0", \
1.74 +"mlo_01C0", \
1.75 +"mnm_19C0", \
1.76 +"mqa_02D0", \
1.77 +"nwr_01D0", \
1.78 +"obs023_06C3", \
1.79 +"opw_01D0", \
1.80 +"orl005_11D2", \
1.81 +"palcbc_30C0", \
1.82 +"pdm_11D0", \
1.83 +"pfa015_01D2", \
1.84 +"poc000_01D1", \
1.85 +"pocn05_01D1", \
1.86 +"pocs05_01D1", \
1.87 +"prs_21C0", \
1.88 +"psa_01D0", \
1.89 +"rpb_01D0", \
1.90 +"rta005_01D2", \
1.91 +"ryo_19C0", \
1.92 +"sbl_06C0", \
1.93 +"sch_23C0", \
1.94 +"scsn03_01D1", \
1.95 +"sey_01D0", \
1.96 +"shm_01D0", \
1.97 +"sis_02D0", \
1.98 +"smo_01C0", \
1.99 +"spo_01C0", \
1.100 +"stm_01D0", \
1.101 +"stp_12D0", \
1.102 +"sum_01D0", \
1.103 +"syo_01D0", \
1.104 +"tap_01D0", \
1.105 +"tdf_01D0", \
1.106 +"trm_11D0", \
1.107 +"uta_01D0", \
1.108 +"uum_01D0", \
1.109 +"wes_23C0", \
1.110 +"wis_01D0", \
1.111 +"wkt009_01C3", \
1.112 +"wlg_01D0", \
1.113 +"wpo000_10D2", \
1.114 +"wpon05_10D2", \
1.115 +"wpos05_10D2", \
1.116 +"yon_19C0", \
1.117 +"zep_01D0", \
1.118 +"zot005_45D2"/)
1.119 +
1.120 + ntime = 12
1.121 + ncase = dimsizes(b)
1.122 + print (ncase)
1.123 +
1.124 + mean = new((/ncase,ntime/),float)
1.125 + lat = new((/ncase/),float)
1.126 + lon = new((/ncase/),float)
1.127 + elev = new((/ncase/),float)
1.128 +
1.129 + do n = 0,dimsizes(b)-1
1.130 + fili = diri + b(n) + "_seas.co2"
1.131 + print (fili)
1.132 +
1.133 + ncol1 = 4
1.134 + data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
1.135 +
1.136 + lat(n) = data1(0,0)
1.137 + lon(n) = data1(0,1)
1.138 + elev(n) = data1(0,2)
1.139 +
1.140 + ncol2 = 5
1.141 + data2 = readAsciiTable(fili, ncol2, "float", 16)
1.142 +
1.143 + mean(n,:) = data2(:,1)
1.144 + end do
1.145 +
1.146 +; mean!0 = "case"
1.147 +; mean!1 = "time"
1.148 +; mean&case= ispan(1,ncase,1)
1.149 +; mean&time= ispan(1,ntime,1)
1.150 + mean@units = "u mol/mol"
1.151 + mean@long_name = "CO2 GlobalView seasonal"
1.152 +
1.153 +; lat!0 = "case"
1.154 +; lon!0 = "case"
1.155 +; elev!0 = "case"
1.156 +; lat&case = mean&case
1.157 +; lon&case = mean&case
1.158 +; elev&case = mean&case
1.159 +
1.160 + lat@units = "degrees_north"
1.161 + lat@long_name = "Latitude"
1.162 + lon@units = "degrees_ease"
1.163 + lon@long_name = "Longitude"
1.164 + elev@units = "m"
1.165 + elev@long_name = "Elevation above sea level"
1.166 + b@units = ""
1.167 + b@long_name = "Station"
1.168 +
1.169 + print (lat)
1.170 + print (mean)
1.171 + print (b)
1.172 +
1.173 + c->LAT = lat
1.174 + c->LON = lon
1.175 + c->ELEV = elev
1.176 + c->CO2_SEAS = mean
1.177 + c->STATION = stringtochar(b)
1.178 +end