diff -r 000000000000 -r 0c6405ab2ff4 npp/52.zonavg_to_T31.ncl
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/npp/52.zonavg_to_T31.ncl	Mon Jan 26 22:08:20 2009 -0500
@@ -0,0 +1,62 @@
+; ***********************************************
+; zonal average plot
+; ***********************************************
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
+;************************************************
+begin
+;************************************************
+; read in observed data
+;************************************************
+ g     = addfile ("Npp_0.05deg_mean.nc","r")
+ bi    = g->NPP   
+ xi    = g->lon 
+ yi    = g->lat
+ yi    = (/ yi(::-1) /)
+ bi    =  (/ bi(::-1,:) /)
+ b2    = bi
+ x2    = xi   
+ 
+ nx = dimsizes(xi)
+ do i= 0,nx-1
+    if (i .lt. 3600) then
+       p = i + 3600
+       xi(p) = x2(i) + 360.      
+    else
+       p = i - 3600
+       xi(p) = x2(i)
+    end if
+    bi(:,p)= b2(:,i) 
+ end do
+;print (xi)
+;print (yi)
+;exit
+
+ f     = addfile ("i01.03cn_1545-1569_ANN_climo.nc","r")
+;f     = addfile ("i01.04casa_1605-1629_ANN_climo.nc","r")
+
+ xo    = f->lon     
+ yo    = f->lat      
+
+ print (xi)
+ print (yi)
+ print (xo)
+ print (yo)
+
+ bo = linint2_Wrap(xi,yi,bi,True,xo,yo,0) 
+ 
+ v = zonalAve(bo)
+
+ v@long_name = "NPP (gC/m2/year)"
+ 
+;************************************************
+; plotting parameters
+;************************************************
+ wks   = gsn_open_wks ("png","xy")                ; open workstation
+
+ res                  = True                     ; plot mods desired
+ res@tiMainString     = "Observed T31"          ; add title
+
+ plot  = gsn_csm_xy (wks,v&lat,v,res) 
+end