1 ; ***********************************************
2 ; interpolate into model grids (T31)
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 ;************************************************
9 ;************************************************
11 ;************************************************
12 diro = "/fis/cgd/cseg/people/jeff/clamp_data/npp/ob/"
13 filo = "Npp_T31_mean.nc"
14 c = addfile(diro+filo,"c")
16 ;************************************************
17 ; read in observed data
18 ;************************************************
19 diri = "/fis/cgd/cseg/people/jeff/clamp_data/npp/ob/"
20 fili = "Npp_0.05deg_mean.nc"
21 g = addfile (diri+fili,"r")
26 ;************************************************
27 ; change into 0-360E, 90S-90N
28 ; Observed NPP*scale_factor
29 ;************************************************
47 bi(:,p)= b2(:,i) * scale_factor
57 ;************************************************
59 ;************************************************
60 diri2 = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
61 fili2 = "b30.061n_1995-2004_ANN_climo_lnd.nc"
62 f = addfile (diri2+fili2,"r")
75 ;bo = linint2_Wrap(xi,yi,bi,True,xo,yo,0)
79 rad = 4.*atan(1.)/180.
82 clat@long_name = "cos(latitude)"
86 bo = new((/nlat,nlon/),float)
88 if (j.eq.0 .or. j.eq.nlat-1) then
91 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
94 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
98 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
99 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
102 ; CLAT = clat({LATS:LATN}) ; do once for *slight* efficiency
103 ; TEMP = bi({LATS:LATN},:) ; 2D [lat,lon]
106 if (i.eq.0 .or. i.eq.nlon-1) then
109 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
111 if (i.eq.nlon-1) then
112 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
116 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
117 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
124 bo(j,i) = avg(bi({LATS:LATN},{LONL:LONR}))
125 ; bo(j,i) = wgt_areaave(TEMP(:,{LONL:LONR}), CLAT, 1.0, 0)
137 bo@long_name = bi@long_name
138 ; bo@_FillValue = bi@_FillValue
139 bo@_FillValue = 1.e+36