Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;************************************************
2 ; output: lat: N->S lon: -180W->180E
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 model_name = "i01.10cn"
15 filo = model_name + "_Fire_C_1979-2004_monthly.nc"
16 c = addfile(diro+filo,"c")
17 filedimdef(c,"year",-1,True)
27 nyear = year_end - year_start + 1
30 data1 = new((/nyear,nmonth,nlat,nlon/),float)
31 data2 = new((/nyear,nmonth,nlat,nlon/),float)
32 data3 = new((/nyear,nmonth,nlat,nlon/),float)
33 date = new((/nyear,nmonth/),integer)
36 diri = "/ptmp/jeff/"+model_name+"/"
42 date(m,n) = year*100 +nn
46 fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
48 fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
53 g = addfile (diri+fili,"r")
55 b1 = g->COL_FIRE_CLOSS
56 data1(m,n,:,:) = (/b1(0,:,:)/)
59 data2(m,n,:,:) = (/b2(0,:,:)/)
62 data3(m,n,:,:) = (/b3(0,:,:)/)
66 ; name dimensions and assign coordinate variables
72 data1@long_name = b1@long_name
73 data1@units = b1@units
74 data1@_FillValue = b1@_FillValue
80 data2@long_name = b2@long_name
81 data2@units = b2@units
82 data2@_FillValue = b2@_FillValue
88 data3@long_name = b3@long_name
89 data3@units = b3@units
90 data3@_FillValue = b3@_FillValue
94 date@long_name = "date: yyyymm"
96 c->COL_FIRE_CLOSS = data1