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 ;************************************************
12 ;************************************************
14 ;************************************************
15 diro = "/fis/cgd/cseg/people/jeff/clamp_data/"
16 ; filo = "LAI_" + year + "_monthly_T42.nc"
17 ; filo = "LAI_2000-2005_mean_T42.nc"
18 filo = "LAI_2000-2005_ensemble_T42.nc"
19 c = addfile(diro+filo,"c")
20 filedimdef(c,"time",-1,True)
22 ;************************************************
23 ; read in observed data
24 ;************************************************
25 diri = "/fis/cgd/cseg/people/jeff/clamp_data/"
26 ; fili = "LAI_" + year + "_monthly.nc"
27 ; fili = "LAI_2000-2005_mean_0.05deg.nc"
28 fili = "LAI_2000-2005_ensemble_0.05deg.nc"
29 g = addfile (diri+fili,"r")
34 ;************************************************
35 ; change into 0-360E, 90S-90N
36 ; Observed NPP*scale_factor
37 ;************************************************
40 bi = (/ bi(:,::-1,:) /)
65 ;************************************************
67 ;************************************************
68 ;fili2 = "i01.04casa_1605-1629_ANN_climo.nc"
69 fili2 = "i01.03cn_1545-1569_ANN_climo.nc"
70 f = addfile (diri+fili2,"r")
83 ;bo = linint2_Wrap(xi,yi,bi,True,xo,yo,0)
87 ; rad = 4.*atan(1.)/180.
90 ; clat@long_name = "cos(latitude)"
92 ; printVarSummary(clat)
94 bo = new((/12,nlat,nlon/),float)
95 ;bo = new((/1,nlat,nlon/),float)
98 if (j.eq.0 .or. j.eq.nlat-1) then
101 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
103 if (j.eq.nlat-1) then
104 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
108 LATS = lat(j)-0.5*(lat(j)-lat(j-1))
109 LATN = lat(j)+0.5*(lat(j+1)-lat(j))
112 ; CLAT = clat({LATS:LATN}) ; do once for *slight* efficiency
113 ; TEMP = bi(:,{LATS:LATN},:) ; 2D [lat,lon]
116 if (i.eq.0 .or. i.eq.nlon-1) then
119 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
121 if (i.eq.nlon-1) then
122 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
126 LONL = lon(i)-0.5*(lon(i)-lon(i-1))
127 LONR = lon(i)+0.5*(lon(i+1)-lon(i))
134 ; bo(:,j,i) = avg(bi(:,{LATS:LATN},{LONL:LONR}))
135 bo(m,j,i) = avg(bi(m,{LATS:LATN},{LONL:LONR}))
136 ; bo(:,j,i) = wgt_areaave(TEMP(:,{LONL:LONR}), CLAT, 1.0, 0)
151 ; bo@units = bi@units
152 ; bo@long_name = bi@long_name
154 bo@long_name = "Leaf Area Index"
155 bo@_FillValue = bi@_FillValue