;************************************************

 ;    Read ascii, Write nc                         

 ;************************************************

 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  

 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  

 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"   

 ;************************************************

 begin

   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"

   b = (/"zot005_45D2_seas.co2", "zot015_45D2_seas.co2"/)

   do n = 0,1

      print (b(n))                   

      ncol1 = 4

      data1 = readAsciiTable(diri+b(n), ncol1, "float", (/10,17/))

      print (data1)

      ncol2 = 5

      data2 = readAsciiTable(diri+b(n), ncol2, "float", 16)

      print (data2)

   end do

 end