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: ;----------------------------------------------------------
forrest@0:   year_start = 2000
forrest@0:   year_end   = 2005
forrest@0: ;----------------------------------------------------------
forrest@0:   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/"
forrest@0:   diro  = "/fis/cgd/cseg/people/jeff/clamp_data/"
forrest@0:   filo1 = "LAI_"+year_start+"-"+year_end+"_ensemble_T42.nc"
forrest@0:   filo2 = "LAI_"+year_start+"-"+year_end+"_mean_T42.nc"
forrest@0:   c = addfile(diro+filo1,"c")
forrest@0:   d = addfile(diro+filo2,"c")
forrest@0: 
forrest@0:   nyear = year_end - year_start + 1
forrest@0: 
forrest@0:   do n = year_start,year_end
forrest@0:    
forrest@0: ;    fili = "LAI_" + n + "_monthly.nc"
forrest@0:      fili = "LAI_" + n + "_monthly_T42.nc"
forrest@0:      print (fili)
forrest@0:   
forrest@0:      a = addfile(diri+fili,"r")
forrest@0:      x = a->LAI
forrest@0:     
forrest@0:      if (n .eq. year_start) then
forrest@0:         y = x
forrest@0:         printVarSummary (x) 
forrest@0:      else
forrest@0:         y = y + x
forrest@0:      end if
forrest@0: 
forrest@0:      delete (a)
forrest@0:      delete (x)
forrest@0:      delete (fili)
forrest@0:   end do
forrest@0:       
forrest@0:   y = y/nyear
forrest@0:   printVarSummary (y)
forrest@0:  
forrest@0:   c->LAI  = y
forrest@0: 
forrest@0:   z = dim_avg_Wrap(y(lat|:,lon|:,time|:))
forrest@0:   printVarSummary (z)
forrest@0:   
forrest@0:   d->LAI  = z     
forrest@0: end
forrest@0: