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_222.nc"
forrest@0: ; c = addfile(diro+filo,"c")
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_seas.co2", \
forrest@0: "aia015_02D2_seas.co2", \
forrest@0: "aia025_02D2_seas.co2", \
forrest@0: "aia035_02D2_seas.co2", \
forrest@0: "aia045_02D2_seas.co2", \
forrest@0: "aia055_02D2_seas.co2", \
forrest@0: "aia065_02D2_seas.co2", \
forrest@0: "alt_01D0_seas.co2", \
forrest@0: "alt_02D0_seas.co2", \
forrest@0: "alt_04D0_seas.co2", \
forrest@0: "alt_06C0_seas.co2", \
forrest@0: "alt_06D0_seas.co2", \
forrest@0: "ams_01D0_seas.co2", \
forrest@0: "ams_11C0_seas.co2", \
forrest@0: "asc_01D0_seas.co2", \
forrest@0: "ask_01D0_seas.co2", \
forrest@0: "avi_01D0_seas.co2", \
forrest@0: "azr_01D0_seas.co2", \
forrest@0: "bal_01D1_seas.co2", \
forrest@0: "bgu_11D0_seas.co2", \
forrest@0: "bhd_15C0_seas.co2", \
forrest@0: "bme_01D0_seas.co2", \
forrest@0: "bmw_01D0_seas.co2", \
forrest@0: "brw_01C0_seas.co2", \
forrest@0: "brw_01D0_seas.co2", \
forrest@0: "bsc_01D0_seas.co2", \
forrest@0: "car030_01D2_seas.co2", \
forrest@0: "car040_01D2_seas.co2", \
forrest@0: "car050_01D2_seas.co2", \
forrest@0: "car060_01D2_seas.co2", \
forrest@0: "car070_01D2_seas.co2", \
forrest@0: "car080_01D2_seas.co2", \
forrest@0: "cba_01D0_seas.co2", \
forrest@0: "cba_04D0_seas.co2", \
forrest@0: "cfa_02D0_seas.co2", \
forrest@0: "cgo_01D0_seas.co2", \
forrest@0: "cgo_02D0_seas.co2", \
forrest@0: "cgo_04D0_seas.co2", \
forrest@0: "chr_01D0_seas.co2", \
forrest@0: "cmn_17C0_seas.co2", \
forrest@0: "cmo_01D0_seas.co2", \
forrest@0: "coi_20C0_seas.co2", \
forrest@0: "cpt_36C0_seas.co2", \
forrest@0: "cri_02D0_seas.co2", \
forrest@0: "crz_01D0_seas.co2", \
forrest@0: "csj_06D0_seas.co2", \
forrest@0: "daa_02D0_seas.co2", \
forrest@0: "eic_01D0_seas.co2", \
forrest@0: "esp005_01D2_seas.co2", \
forrest@0: "esp015_01D2_seas.co2", \
forrest@0: "esp025_01D2_seas.co2", \
forrest@0: "esp035_01D2_seas.co2", \
forrest@0: "esp045_01D2_seas.co2", \
forrest@0: "esp055_01D2_seas.co2", \
forrest@0: "esp_02D0_seas.co2", \
forrest@0: "esp_06D0_seas.co2", \
forrest@0: "frd040_06C3_seas.co2", \
forrest@0: "gmi_01D0_seas.co2", \
forrest@0: "goz_01D0_seas.co2", \
forrest@0: "gsn_24D0_seas.co2", \
forrest@0: "haa005_01D2_seas.co2", \
forrest@0: "haa015_01D2_seas.co2", \
forrest@0: "haa025_01D2_seas.co2", \
forrest@0: "haa035_01D2_seas.co2", \
forrest@0: "haa045_01D2_seas.co2", \
forrest@0: "haa055_01D2_seas.co2", \
forrest@0: "haa065_01D2_seas.co2", \
forrest@0: "haa075_01D2_seas.co2", \
forrest@0: "hat_20C0_seas.co2", \
forrest@0: "hba_01D0_seas.co2", \
forrest@0: "hfm005_01D2_seas.co2", \
forrest@0: "hfm015_01D2_seas.co2", \
forrest@0: "hfm025_01D2_seas.co2", \
forrest@0: "hfm035_01D2_seas.co2", \
forrest@0: "hfm045_01D2_seas.co2", \
forrest@0: "hfm055_01D2_seas.co2", \
forrest@0: "hfm065_01D2_seas.co2", \
forrest@0: "hfm075_01D2_seas.co2", \
forrest@0: "hun010_35C3_seas.co2", \
forrest@0: "hun048_35C3_seas.co2", \
forrest@0: "hun082_35C3_seas.co2", \
forrest@0: "hun115_35C3_seas.co2", \
forrest@0: "hun_01D0_seas.co2", \
forrest@0: "ice_01D0_seas.co2", \
forrest@0: "itn051_01C3_seas.co2", \
forrest@0: "itn123_01C3_seas.co2", \
forrest@0: "itn496_01C3_seas.co2", \
forrest@0: "itn_01D0_seas.co2", \
forrest@0: "izo_01D0_seas.co2", \
forrest@0: "izo_27C0_seas.co2", \
forrest@0: "jbn_29C0_seas.co2", \
forrest@0: "key_01D0_seas.co2", \
forrest@0: "kum_01D0_seas.co2", \
forrest@0: "kum_04D0_seas.co2", \
forrest@0: "kzd_01D0_seas.co2", \
forrest@0: "kzm_01D0_seas.co2", \
forrest@0: "lef011_01C3_seas.co2", \
forrest@0: "lef030_01C3_seas.co2", \
forrest@0: "lef076_01C3_seas.co2", \
forrest@0: "lef122_01C3_seas.co2", \
forrest@0: "lef244_01C3_seas.co2", \
forrest@0: "lef396_01C3_seas.co2", \
forrest@0: "lef_01D0_seas.co2", \
forrest@0: "ljo_04D0_seas.co2", \
forrest@0: "lmp_28D0_seas.co2", \
forrest@0: "maa_02D0_seas.co2", \
forrest@0: "mbc_01D0_seas.co2", \
forrest@0: "mhd_01D0_seas.co2", \
forrest@0: "mhdcbc_11C0_seas.co2", \
forrest@0: "mhdrbc_11C0_seas.co2", \
forrest@0: "mid_01D0_seas.co2", \
forrest@0: "mlo_01C0_seas.co2", \
forrest@0: "mlo_01D0_seas.co2", \
forrest@0: "mlo_02D0_seas.co2", \
forrest@0: "mlo_04D0_seas.co2", \
forrest@0: "mnm_19C0_seas.co2", \
forrest@0: "mqa_02D0_seas.co2", \
forrest@0: "nwr_01D0_seas.co2", \
forrest@0: "obs023_06C3_seas.co2", \
forrest@0: "opw_01D0_seas.co2", \
forrest@0: "orl005_11D2_seas.co2", \
forrest@0: "orl015_11D2_seas.co2", \
forrest@0: "orl025_11D2_seas.co2", \
forrest@0: "orl035_11D2_seas.co2", \
forrest@0: "palcbc_30C0_seas.co2", \
forrest@0: "palmbc_30C0_seas.co2", \
forrest@0: "pdm_11D0_seas.co2", \
forrest@0: "pfa015_01D2_seas.co2", \
forrest@0: "pfa025_01D2_seas.co2", \
forrest@0: "pfa035_01D2_seas.co2", \
forrest@0: "pfa045_01D2_seas.co2", \
forrest@0: "pfa055_01D2_seas.co2", \
forrest@0: "pfa065_01D2_seas.co2", \
forrest@0: "pfa075_01D2_seas.co2", \
forrest@0: "poc000_01D1_seas.co2", \
forrest@0: "pocn05_01D1_seas.co2", \
forrest@0: "pocn10_01D1_seas.co2", \
forrest@0: "pocn15_01D1_seas.co2", \
forrest@0: "pocn20_01D1_seas.co2", \
forrest@0: "pocn25_01D1_seas.co2", \
forrest@0: "pocn30_01D1_seas.co2", \
forrest@0: "pocn35_01D1_seas.co2", \
forrest@0: "pocn40_01D1_seas.co2", \
forrest@0: "pocn45_01D1_seas.co2", \
forrest@0: "pocs05_01D1_seas.co2", \
forrest@0: "pocs10_01D1_seas.co2", \
forrest@0: "pocs15_01D1_seas.co2", \
forrest@0: "pocs20_01D1_seas.co2", \
forrest@0: "pocs25_01D1_seas.co2", \
forrest@0: "pocs30_01D1_seas.co2", \
forrest@0: "pocs35_01D1_seas.co2", \
forrest@0: "prs_21C0_seas.co2", \
forrest@0: "prs_21D0_seas.co2", \
forrest@0: "psa_01D0_seas.co2", \
forrest@0: "psa_04D0_seas.co2", \
forrest@0: "rpb_01D0_seas.co2", \
forrest@0: "rta005_01D2_seas.co2", \
forrest@0: "rta015_01D2_seas.co2", \
forrest@0: "rta025_01D2_seas.co2", \
forrest@0: "rta035_01D2_seas.co2", \
forrest@0: "rta045_01D2_seas.co2", \
forrest@0: "ryo_19C0_seas.co2", \
forrest@0: "sbl_06C0_seas.co2", \
forrest@0: "sbl_06D0_seas.co2", \
forrest@0: "sch_23C0_seas.co2", \
forrest@0: "scsn03_01D1_seas.co2", \
forrest@0: "scsn06_01D1_seas.co2", \
forrest@0: "scsn09_01D1_seas.co2", \
forrest@0: "scsn12_01D1_seas.co2", \
forrest@0: "scsn15_01D1_seas.co2", \
forrest@0: "scsn18_01D1_seas.co2", \
forrest@0: "scsn21_01D1_seas.co2", \
forrest@0: "sey_01D0_seas.co2", \
forrest@0: "shm_01D0_seas.co2", \
forrest@0: "sis_02D0_seas.co2", \
forrest@0: "smo_01C0_seas.co2", \
forrest@0: "smo_01D0_seas.co2", \
forrest@0: "smo_04D0_seas.co2", \
forrest@0: "spo_01C0_seas.co2", \
forrest@0: "spo_01D0_seas.co2", \
forrest@0: "spo_02D0_seas.co2", \
forrest@0: "spo_04D0_seas.co2", \
forrest@0: "stm_01D0_seas.co2", \
forrest@0: "stmebc_01D0_seas.co2", \
forrest@0: "stp_12D0_seas.co2", \
forrest@0: "sum_01D0_seas.co2", \
forrest@0: "syo_01D0_seas.co2", \
forrest@0: "syo_09C0_seas.co2", \
forrest@0: "tap_01D0_seas.co2", \
forrest@0: "tdf_01D0_seas.co2", \
forrest@0: "trm_11D0_seas.co2", \
forrest@0: "uta_01D0_seas.co2", \
forrest@0: "uum_01D0_seas.co2", \
forrest@0: "wes_23C0_seas.co2", \
forrest@0: "wis_01D0_seas.co2", \
forrest@0: "wkt009_01C3_seas.co2", \
forrest@0: "wkt030_01C3_seas.co2", \
forrest@0: "wkt061_01C3_seas.co2", \
forrest@0: "wkt122_01C3_seas.co2", \
forrest@0: "wkt244_01C3_seas.co2", \
forrest@0: "wkt457_01C3_seas.co2", \
forrest@0: "wlg_01D0_seas.co2", \
forrest@0: "wlg_33C0_seas.co2", \
forrest@0: "wpo000_10D2_seas.co2", \
forrest@0: "wpon05_10D2_seas.co2", \
forrest@0: "wpon10_10D2_seas.co2", \
forrest@0: "wpon15_10D2_seas.co2", \
forrest@0: "wpon20_10D2_seas.co2", \
forrest@0: "wpon25_10D2_seas.co2", \
forrest@0: "wpon30_10D2_seas.co2", \
forrest@0: "wpos05_10D2_seas.co2", \
forrest@0: "wpos10_10D2_seas.co2", \
forrest@0: "wpos15_10D2_seas.co2", \
forrest@0: "wpos20_10D2_seas.co2", \
forrest@0: "wpos25_10D2_seas.co2", \
forrest@0: "yon_19C0_seas.co2", \
forrest@0: "zep_01D0_seas.co2", \
forrest@0: "zep_31C0_seas.co2", \
forrest@0: "zot005_45D2_seas.co2", \
forrest@0: "zot015_45D2_seas.co2", \
forrest@0: "zot025_45D2_seas.co2", \
forrest@0: "zot035_45D2_seas.co2"/)
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)
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:      print (lat)
forrest@0:      print (mean)
forrest@0: 
forrest@0:   
forrest@0: end