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 ;************************************************
11 ;************************************************
12 ; read in observed data
13 ;************************************************
14 diri = "/fis/cgd/cseg/people/jeff/clamp_data/biomass/ob/"
15 fili = "amazon_biomass.nc"
16 g = addfile (diri+fili,"r")
22 ;************************************************
24 ;************************************************
25 diri2 = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
26 fili2 = "b30.061n_1995-2004_ANN_climo_lnd.nc"
27 f = addfile (diri2+fili2,"r")
37 ind_lonL = min(ind(lon .ge. xi(0)))
38 ind_lonR = max(ind(lon .le. xi(dimsizes(xi)-1)))
45 print (xi(dimsizes(xi)-1))
48 ind_latS = min(ind(lat .ge. yi(0)))
49 ind_latN = max(ind(lat .le. yi(dimsizes(yi)-1)))
56 print (yi(dimsizes(yi)-1))
59 bo = new((/nlat,nlon/),float)
61 if (j.eq.0 .or. j.eq.nlat-1) then
64 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
67 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
71 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
72 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
75 ; CLAT = clat({LATS:LATN}) ; do once for *slight* efficiency
76 ; TEMP = bi({LATS:LATN},:) ; 2D [lat,lon]
79 if (i.eq.0 .or. i.eq.nlon-1) then
82 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
85 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
89 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
90 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
97 bo(j,i) = avg(bi({LATS:LATN},{LONL:LONR}))
98 ; bo(j,i) = wgt_areaave(TEMP(:,{LONL:LONR}), CLAT, 1.0, 0)
110 bo@long_name = bi@long_name
111 ; bo@_FillValue = bi@_FillValue
112 bo@_FillValue = 1.e+36