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.39cn"
16 diro = "/lustre/wolf-ddn/scratch/hof/cases1/i01.39cn_out/"
17 filo = model_name + "_ameriflux_2101-2125_monthly.nc"
18 c = addfile(diro+filo,"c")
19 filedimdef(c,"year",-1,True)
29 nyear = year_end - year_start + 1
32 data1 = new((/nyear,nmonth,nlat,nlon/),float)
33 data2 = new((/nyear,nmonth,nlat,nlon/),float)
34 data3 = new((/nyear,nmonth,nlat,nlon/),float)
35 ;data4 = new((/nyear,nmonth,nlat,nlon/),float)
36 ;data5 = new((/nyear,nmonth,nlat,nlon/),float)
37 data6 = new((/nyear,nmonth,nlat,nlon/),float)
38 data7 = new((/nyear,nmonth,nlat,nlon/),float)
39 data8 = new((/nyear,nmonth,nlat,nlon/),float)
40 data9 = new((/nyear,nmonth,nlat,nlon/),float)
41 data10 = new((/nyear,nmonth,nlat,nlon/),float)
42 date = new((/nyear,nmonth/),integer)
44 ; name dimensions and assign coordinate variables
98 date@long_name = "date: yyyymm"
101 diri = "/lustre/wolf-ddn/scratch/hof/cases1/"+model_name+"/"
104 year = year_start + m
107 date(m,n) = year*100 +nn
111 fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
113 fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
118 g = addfile (diri+fili,"r")
121 data1(m,n,:,:) = (/b(0,:,:)/)
122 data1@long_name = b@long_name
123 data1@units = b@units
124 data1@_FillValue = b@_FillValue
128 data2(m,n,:,:) = (/b(0,:,:)/)
129 data2@long_name = b@long_name
130 data2@units = b@units
131 data2@_FillValue = b@_FillValue
135 data3(m,n,:,:) = (/b(0,:,:)/)
136 data3@long_name = b@long_name
137 data3@units = b@units
138 data3@_FillValue = b@_FillValue
142 ;data4(m,n,:,:) = (/b(0,:,:)/)
143 ;data4@long_name = b@long_name
144 ;data4@units = b@units
145 ;data4@_FillValue = b@_FillValue
149 ;data5(m,n,:,:) = (/b(0,:,:)/)
150 ;data5@long_name = b@long_name
151 ;data5@units = b@units
152 ;data5@_FillValue = b@_FillValue
156 data6(m,n,:,:) = (/b(0,:,:)/)
157 data6@long_name = b@long_name
158 data6@units = b@units
159 data6@_FillValue = b@_FillValue
163 data7(m,n,:,:) = (/b(0,:,:)/)
164 data7@long_name = b@long_name
165 data7@units = b@units
166 data7@_FillValue = b@_FillValue
170 data8(m,n,:,:) = (/b(0,:,:)/)
171 data8@long_name = b@long_name
172 data8@units = b@units
173 data8@_FillValue = b@_FillValue
177 data9(m,n,:,:) = (/b(0,:,:)/)
178 data9@long_name = b@long_name
179 data9@units = b@units
180 data9@_FillValue = b@_FillValue
184 data10(m,n,:,:) = (/b(0,:,:)/)
185 data10@long_name = b@long_name
186 data10@units = b@units
187 data10@_FillValue = b@_FillValue