1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/co2/02.test.read.ncl	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,22 @@
     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 +  diri  = "/fis/cgd/cseg/people/jeff/clamp_data/co2/globalView-C4Mip/seas/"
    1.13 +  b = (/"zot005_45D2_seas.co2", "zot015_45D2_seas.co2"/)
    1.14 +  do n = 0,1
    1.15 +     print (b(n))                   
    1.16 +     ncol1 = 4
    1.17 +     data1 = readAsciiTable(diri+b(n), ncol1, "float", (/10,17/))
    1.18 +     print (data1)
    1.19 +
    1.20 +     ncol2 = 5
    1.21 +     data2 = readAsciiTable(diri+b(n), ncol2, "float", 16)
    1.22 +     print (data2)
    1.23 +  end do
    1.24 +  
    1.25 +end