forrest@0: ;************************************************
forrest@0: ;    Read ascii, Write nc                         
forrest@0: ;************************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"   
forrest@0: ;************************************************
forrest@0: begin
forrest@0: ; final output
forrest@0:   diro  = "/fis/cgd/cseg/people/jeff/clamp_data/co2/"
forrest@0:   filo =  "co2_globalView_98.nc"
forrest@0:   c = addfile(diro+filo,"c")
forrest@0:   filedimdef(c,"time",-1,True)
forrest@0: 
forrest@0: ; input ascii files
forrest@0:   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
forrest@0:   b = (/ \
forrest@0: "aia005_02D2", \
forrest@0: "alt_01D0", \
forrest@0: "ams_01D0", \
forrest@0: "asc_01D0", \
forrest@0: "ask_01D0", \
forrest@0: "avi_01D0", \
forrest@0: "azr_01D0", \
forrest@0: "bal_01D1", \
forrest@0: "bgu_11D0", \
forrest@0: "bhd_15C0", \
forrest@0: "bme_01D0", \
forrest@0: "bmw_01D0", \
forrest@0: "brw_01C0", \
forrest@0: "bsc_01D0", \
forrest@0: "car030_01D2", \
forrest@0: "cba_01D0", \
forrest@0: "cfa_02D0", \
forrest@0: "cgo_01D0", \
forrest@0: "chr_01D0", \
forrest@0: "cmn_17C0", \
forrest@0: "cmo_01D0", \
forrest@0: "coi_20C0", \
forrest@0: "cpt_36C0", \
forrest@0: "cri_02D0", \
forrest@0: "crz_01D0", \
forrest@0: "csj_06D0", \
forrest@0: "daa_02D0", \
forrest@0: "eic_01D0", \
forrest@0: "esp_02D0", \
forrest@0: "frd040_06C3", \
forrest@0: "gmi_01D0", \
forrest@0: "goz_01D0", \
forrest@0: "gsn_24D0", \
forrest@0: "haa005_01D2", \
forrest@0: "hat_20C0", \
forrest@0: "hba_01D0", \
forrest@0: "hfm005_01D2", \
forrest@0: "hun010_35C3", \
forrest@0: "ice_01D0", \
forrest@0: "itn051_01C3", \
forrest@0: "izo_01D0", \
forrest@0: "jbn_29C0", \
forrest@0: "key_01D0", \
forrest@0: "kum_01D0", \
forrest@0: "kzd_01D0", \
forrest@0: "kzm_01D0", \
forrest@0: "lef011_01C3", \
forrest@0: "ljo_04D0", \
forrest@0: "lmp_28D0", \
forrest@0: "maa_02D0", \
forrest@0: "mbc_01D0", \
forrest@0: "mhd_01D0", \
forrest@0: "mid_01D0", \
forrest@0: "mlo_01C0", \
forrest@0: "mnm_19C0", \
forrest@0: "mqa_02D0", \
forrest@0: "nwr_01D0", \
forrest@0: "obs023_06C3", \
forrest@0: "opw_01D0", \
forrest@0: "orl005_11D2", \
forrest@0: "palcbc_30C0", \
forrest@0: "pdm_11D0", \
forrest@0: "pfa015_01D2", \
forrest@0: "poc000_01D1", \
forrest@0: "pocn05_01D1", \
forrest@0: "pocs05_01D1", \
forrest@0: "prs_21C0", \
forrest@0: "psa_01D0", \
forrest@0: "rpb_01D0", \
forrest@0: "rta005_01D2", \
forrest@0: "ryo_19C0", \
forrest@0: "sbl_06C0", \
forrest@0: "sch_23C0", \
forrest@0: "scsn03_01D1", \
forrest@0: "sey_01D0", \
forrest@0: "shm_01D0", \
forrest@0: "sis_02D0", \
forrest@0: "smo_01C0", \
forrest@0: "spo_01C0", \
forrest@0: "stm_01D0", \
forrest@0: "stp_12D0", \
forrest@0: "sum_01D0", \
forrest@0: "syo_01D0", \
forrest@0: "tap_01D0", \
forrest@0: "tdf_01D0", \
forrest@0: "trm_11D0", \
forrest@0: "uta_01D0", \
forrest@0: "uum_01D0", \
forrest@0: "wes_23C0", \
forrest@0: "wis_01D0", \
forrest@0: "wkt009_01C3", \
forrest@0: "wlg_01D0", \
forrest@0: "wpo000_10D2", \
forrest@0: "wpon05_10D2", \
forrest@0: "wpos05_10D2", \
forrest@0: "yon_19C0", \
forrest@0: "zep_01D0", \
forrest@0: "zot005_45D2"/) 
forrest@0: 
forrest@0:   ntime = 12
forrest@0:   ncase = dimsizes(b)
forrest@0:   print (ncase)
forrest@0: 
forrest@0:   mean = new((/ncase,ntime/),float)
forrest@0:   lat  = new((/ncase/),float)
forrest@0:   lon  = new((/ncase/),float)
forrest@0:   elev = new((/ncase/),float) 
forrest@0: 
forrest@0:   do n = 0,dimsizes(b)-1
forrest@0:      fili = diri + b(n) + "_seas.co2"
forrest@0:      print (fili)
forrest@0:                     
forrest@0:      ncol1 = 4
forrest@0:      data1 = readAsciiTable(fili, ncol1, "float", (/10,17/))
forrest@0: 
forrest@0:      lat(n)  = data1(0,0)
forrest@0:      lon(n)  = data1(0,1)
forrest@0:      elev(n) = data1(0,2)
forrest@0: 
forrest@0:      ncol2 = 5
forrest@0:      data2 = readAsciiTable(fili, ncol2, "float", 16)
forrest@0: 
forrest@0:      mean(n,:) = data2(:,1)
forrest@0:   end do
forrest@0: 
forrest@0: ; mean!0  = "case"
forrest@0: ; mean!1  = "time"
forrest@0: ; mean&case= ispan(1,ncase,1)
forrest@0: ; mean&time= ispan(1,ntime,1)
forrest@0:   mean@units      = "u mol/mol"
forrest@0:   mean@long_name  = "CO2 GlobalView seasonal"
forrest@0: 
forrest@0: ; lat!0     = "case"
forrest@0: ; lon!0     = "case"
forrest@0: ; elev!0    = "case"
forrest@0: ; lat&case  = mean&case
forrest@0: ; lon&case  = mean&case
forrest@0: ; elev&case = mean&case
forrest@0: 
forrest@0:   lat@units      = "degrees_north"
forrest@0:   lat@long_name  = "Latitude"
forrest@0:   lon@units      = "degrees_ease"
forrest@0:   lon@long_name  = "Longitude"        
forrest@0:   elev@units      = "m"
forrest@0:   elev@long_name  = "Elevation above sea level"
forrest@0:   b@units      = ""
forrest@0:   b@long_name  = "Station"
forrest@0: 
forrest@0:   print (lat)
forrest@0:   print (mean)
forrest@0:   print (b)
forrest@0:   
forrest@0:   c->LAT       = lat
forrest@0:   c->LON       = lon
forrest@0:   c->ELEV      = elev
forrest@0:   c->CO2_SEAS  = mean
forrest@0:   c->STATION   = stringtochar(b)  
forrest@0: end