; ***********************************************

 ; interpolate into model grids (T31)

 ; ***********************************************

 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

 ;************************************************

 ; output data

 ;************************************************

   diro  = "/fis/cgd/cseg/people/jeff/clamp/"

   filo  = "Npp_T31_mean.nc"

   c = addfile(diro+filo,"c")

 ;************************************************

 ; read in observed data

 ;************************************************

  g     = addfile ("Npp_0.05deg_mean.nc","r")

  bi    = g->NPP   

  xi    = g->lon 

  yi    = g->lat

 ;************************************************

 ; change into 0-360E, 90S-90N

 ; Observed NPP*scale_factor

 ;************************************************

  scale_factor = 0.1

  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) * scale_factor

  end do

 ;print (xi)

 ;print (yi)

 ;exit

 ;************************************************

 ; read in model data

 ;************************************************

  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) 

   bo!0   = "lat"

   bo!1   = "lon"

   bo&lat =  yo

   bo&lon =  xo

   bo@units      = bi@units

   bo@long_name  = bi@long_name

   bo@_FillValue = bi@_FillValue

   c->NPP  = bo

 end