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 ("i01.03cn_1545-1569_ANN_climo.nc","r")
forrest@0:  f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
forrest@0:  u     = f->NPP 
forrest@0:  v     = zonalAve(u)
forrest@0: 
forrest@0:  nsec_per_year = 60*60*24*365
forrest@0: 
forrest@0:  v = v * nsec_per_year
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     = "Model i01.03cn"         
forrest@0:  res@tiMainString     = "Model i01.04casa"         
forrest@0: 
forrest@0:  plot  = gsn_csm_xy (wks,v&lat,v(0,:),res) 
forrest@0: end