1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/npp/51.zonavg.ncl	Mon Jan 26 22:08:20 2009 -0500
     1.3 @@ -0,0 +1,32 @@
     1.4 +; ***********************************************
     1.5 +; xy_1.ncl
     1.6 +; ***********************************************
     1.7 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
     1.8 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
     1.9 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
    1.10 +;************************************************
    1.11 +begin
    1.12 +;************************************************
    1.13 +; read in data
    1.14 +;************************************************
    1.15 +;f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
    1.16 + f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
    1.17 + u     = f->NPP 
    1.18 + v     = zonalAve(u)
    1.19 +
    1.20 + nsec_per_year = 60*60*24*365
    1.21 +
    1.22 + v = v * nsec_per_year
    1.23 +
    1.24 + v@long_name = "NPP (gC/m2/year)"        
    1.25 +;************************************************
    1.26 +; plotting parameters
    1.27 +;************************************************
    1.28 + wks   = gsn_open_wks ("png","xy")                
    1.29 +
    1.30 + res                  = True                     
    1.31 +;res@tiMainString     = "Model i01.03cn"         
    1.32 + res@tiMainString     = "Model i01.04casa"         
    1.33 +
    1.34 + plot  = gsn_csm_xy (wks,v&lat,v(0,:),res) 
    1.35 +end
    1.36 \ No newline at end of file