Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;************************************************
3 ; output: lat: N->S lon: -180W->180E
4 ;************************************************
5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
6 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
7 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
8 ;************************************************
13 model_name = "i01.10casa"
16 filo = model_name + "_ameriflux_1990-2004_monthly.nc"
17 c = addfile(diro+filo,"c")
18 filedimdef(c,"year",-1,True)
28 nyear = year_end - year_start + 1
31 data1 = new((/nyear,nmonth,nlat,nlon/),float)
32 data2 = new((/nyear,nmonth,nlat,nlon/),float)
33 data3 = new((/nyear,nmonth,nlat,nlon/),float)
34 data4 = new((/nyear,nmonth,nlat,nlon/),float)
35 data5 = new((/nyear,nmonth,nlat,nlon/),float)
36 data6 = new((/nyear,nmonth,nlat,nlon/),float)
37 data7 = new((/nyear,nmonth,nlat,nlon/),float)
38 date = new((/nyear,nmonth/),integer)
40 ; name dimensions and assign coordinate variables
79 date@long_name = "date: yyyymm"
82 diri = "/ptmp/jeff/"+model_name+"/"
88 date(m,n) = year*100 +nn
92 fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
94 fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
99 g = addfile (diri+fili,"r")
102 data1(m,n,:,:) = (/b(0,:,:)/)
103 data1@long_name = b@long_name
104 data1@units = b@units
105 data1@_FillValue = b@_FillValue
109 data2(m,n,:,:) = (/b(0,:,:)/)
110 data2@long_name = b@long_name
111 data2@units = b@units
112 data2@_FillValue = b@_FillValue
116 data3(m,n,:,:) = (/b(0,:,:)/)
118 data3(m,n,:,:) = data3(m,n,:,:)+(/b(0,:,:)/)
119 data3@long_name = "total ecosystem respiration, autotrophic + heterotrophic"
120 data3@units = b@units
121 data3@_FillValue = b@_FillValue
125 data4(m,n,:,:) = (/b(0,:,:)/)
126 data4@long_name = b@long_name
127 data4@units = b@units
128 data4@_FillValue = b@_FillValue
132 data5(m,n,:,:) = (/b(0,:,:)/)
133 data5@long_name = b@long_name
134 data5@units = b@units
135 data5@_FillValue = b@_FillValue
139 data6(m,n,:,:) = (/b(0,:,:)/)
140 data6@long_name = b@long_name
141 data6@units = b@units
142 data6@_FillValue = b@_FillValue
146 data7(m,n,:,:) = (/b(0,:,:)/)
147 data7@long_name = b@long_name
148 data7@units = b@units
149 data7@_FillValue = b@_FillValue