npp/51.zonavg_ob.ncl
changeset 1 4be95183fbcd
equal deleted inserted replaced
-1:000000000000 0:6674eb456744
       
     1 ; ***********************************************
       
     2 ; xy_1.ncl
       
     3 ; ***********************************************
       
     4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
       
     5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
       
     6 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
       
     7 ;************************************************
       
     8 begin
       
     9 ;************************************************
       
    10 ; read in data
       
    11 ;************************************************
       
    12  f     = addfile ("Npp_T31_mean.nc","r")
       
    13  u     = f->NPP 
       
    14  v     = zonalAve(u)
       
    15 
       
    16  v@long_name = "NPP (gC/m2/year)"        
       
    17 ;************************************************
       
    18 ; plotting parameters
       
    19 ;************************************************
       
    20  wks   = gsn_open_wks ("png","xy")                
       
    21 
       
    22  res                  = True                              
       
    23  res@tiMainString     = "Observed"         
       
    24 
       
    25  plot  = gsn_csm_xy (wks,v&lat,v,res) 
       
    26 end