forrest@0: ; ***********************************************
forrest@0: ; zonal average plot
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 observed data
forrest@0: ;************************************************
forrest@0:  g     = addfile ("Npp_0.05deg_mean.nc","r")
forrest@0:  bi    = g->NPP   
forrest@0:  xi    = g->lon 
forrest@0:  yi    = g->lat
forrest@0:  yi    = (/ yi(::-1) /)
forrest@0:  bi    =  (/ bi(::-1,:) /)
forrest@0:  b2    = bi
forrest@0:  x2    = xi   
forrest@0:  
forrest@0:  nx = dimsizes(xi)
forrest@0:  do i= 0,nx-1
forrest@0:     if (i .lt. 3600) then
forrest@0:        p = i + 3600
forrest@0:        xi(p) = x2(i) + 360.      
forrest@0:     else
forrest@0:        p = i - 3600
forrest@0:        xi(p) = x2(i)
forrest@0:     end if
forrest@0:     bi(:,p)= b2(:,i) 
forrest@0:  end do
forrest@0: ;print (xi)
forrest@0: ;print (yi)
forrest@0: ;exit
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: 
forrest@0:  xo    = f->lon     
forrest@0:  yo    = f->lat      
forrest@0: 
forrest@0:  print (xi)
forrest@0:  print (yi)
forrest@0:  print (xo)
forrest@0:  print (yo)
forrest@0: 
forrest@0:  bo = linint2_Wrap(xi,yi,bi,True,xo,yo,0) 
forrest@0:  
forrest@0:  v = zonalAve(bo)
forrest@0: 
forrest@0:  v@long_name = "NPP (gC/m2/year)"
forrest@0:  
forrest@0: ;************************************************
forrest@0: ; plotting parameters
forrest@0: ;************************************************
forrest@0:  wks   = gsn_open_wks ("png","xy")                ; open workstation
forrest@0: 
forrest@0:  res                  = True                     ; plot mods desired
forrest@0:  res@tiMainString     = "Observed T31"          ; add title
forrest@0: 
forrest@0:  plot  = gsn_csm_xy (wks,v&lat,v,res) 
forrest@0: end