1 ; ***********************************************
2 ; interpolate into model grids (1.9x2.5)
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_1.9_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
55 ;************************************************
57 ;************************************************
58 diri2 = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
59 fili2 = "newcn05_ncep_1i_MONS_climo_lnd.nc"
60 f = addfile (diri2+fili2,"r")
72 bo = new((/nlat,nlon/),float)
75 if (j.eq.0 .or. j.eq.nlat-1) then
78 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
81 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
85 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
86 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
90 if (i.eq.0 .or. i.eq.nlon-1) then
93 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
96 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
100 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
101 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
109 bo(j,i) = avg(bi({LATS:LATN},{LONL:LONR}))
118 bo@long_name = bi@long_name
119 ; bo@_FillValue = bi@_FillValue
120 bo@_FillValue = 1.e+36