Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
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 ;************************************************
10 ;************************************************
12 ;************************************************
18 diri = "/fis/cgd/cseg/people/jeff/clamp_data/surface_model/"
19 fili = "lnd_"+ model_grid +".nc"
20 f = addfile (diri+fili,"r")
27 ;************************************************
29 ;************************************************
30 diro = "/fis/cgd/cseg/people/jeff/clamp_data/observed/biomass/"
31 filo = "amazon_mask_"+model_grid + ".nc"
32 c = addfile(diro+filo,"c")
34 ;************************************************
35 ; read 1x1 amazon_mask data
36 ;************************************************
37 diri = "/fis/cgd/cseg/people/jeff/clamp_data/observed/biomass/"
38 fili = "amazon_mask_1x1.nc"
39 g = addfile(diri+fili,"r")
45 ;************************************************
46 ; change from -180-180, 90N-90S into 0-360E, 90S-90N
47 ; Observed NPP*scale_factor
48 ;************************************************
74 ;bo = linint2_Wrap(xi,yi,bi,True,xo,yo,0)
78 rad = 4.*atan(1.)/180.
82 clat@long_name = "cos(latitude)"
84 ; printVarSummary(clat)
86 bo = new ((/nlat,nlon/),"float")
89 if (j.eq.0 .or. j.eq.nlat-1) then
92 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
95 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
99 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
100 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
103 ; CLAT = clat({LATS:LATN}) ; do once for *slight* efficiency
104 ; TEMP = bi({LATS:LATN},:) ; 2D [lat,lon]
107 if (i.eq.0 .or. i.eq.nlon-1) then
110 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
112 if (i.eq.nlon-1) then
113 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
117 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
118 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
126 bo(j,i) = avg(bi({LATS:LATN},{LONL:LONR}))
127 ; bo(j,i) = wgt_areaave(TEMP(:,{LONL:LONR}), CLAT, 1.0, 0)
135 ; bo = where(bo.ge.0.5, 1.0,0.)
141 ; bo@units = bi@units
142 bo@long_name = "amazon mask"
143 ; bo@_FillValue = bi@_FillValue
144 bo@_FillValue = 1.e+36