forrest@0: ; ***********************************************
forrest@0: ; xy_4.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:  g     = addfile ("Npp_T31_mean.nc","r")
forrest@0:  c     = g->NPP
forrest@0:  u     = zonalAve(c)  
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:  b     = f->NPP      
forrest@0: 
forrest@0:  sec_to_year = 86400.*365.
forrest@0: 
forrest@0:  b = b * sec_to_year
forrest@0:  v = zonalAve(b)
forrest@0: 
forrest@0:  b@long_name = "NPP (gC/m2/year)"
forrest@0:  
forrest@0:  ccr = esccr(u,v(0,:),0)
forrest@0: 
forrest@0:  print (ccr)
forrest@0: 
forrest@0: end