forrest@0: ;*************************************************
forrest@0: ; ce_1.ncl
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: ;************************************************
forrest@0: begin
forrest@0:  
forrest@0: ;************************************************
forrest@0: ; read in model data
forrest@0: ;************************************************
forrest@0:   diri  = "/fis/cgd/cseg/people/jeff/clamp/class/"
forrest@0:   fili  = "T42.clamp.surface-data.060412.nc"
forrest@0:   f     = addfile(diri+fili,"r")
forrest@0: 
forrest@0:   y     = f->PCT_PFT
forrest@0:   printVarSummary(y)
forrest@0: 
forrest@0:   x = y(0,:,:)
forrest@0:   
forrest@0:   dimx = dimsizes(y)
forrest@0:   nlat = dimx(1)
forrest@0:   nlon = dimx(2)
forrest@0: 
forrest@0:   do j= 0,nlat-1
forrest@0:   do i= 0,nlon-1
forrest@0:      x(j,i) = sum(y(:,j,i))      
forrest@0:   end do
forrest@0:   end do 
forrest@0: 
forrest@0:   print (x)
forrest@0: 
forrest@0: end