forrest@0: ; ***********************************************
forrest@0: ; xy_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: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
forrest@0: ;************************************************
forrest@0: begin
forrest@0: ;************************************************
forrest@0: ; read in data
forrest@0: ;************************************************
forrest@0:  f     = addfile ("Npp_T31_mean.nc","r")
forrest@0:  u     = f->NPP 
forrest@0:  v     = zonalAve(u)
forrest@0: 
forrest@0:  v@long_name = "NPP (gC/m2/year)"        
forrest@0: ;************************************************
forrest@0: ; plotting parameters
forrest@0: ;************************************************
forrest@0:  wks   = gsn_open_wks ("png","xy")                
forrest@0: 
forrest@0:  res                  = True                              
forrest@0:  res@tiMainString     = "Observed"         
forrest@0: 
forrest@0:  plot  = gsn_csm_xy (wks,v&lat,v,res) 
forrest@0: end