Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;************************************************
3 ; output fire(COL_FIRE_CLOSS), NPP, NEE
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"
21 diri = "/ptmp/jeff/"+model_name+"/"
28 filo = model_name + "_Fire_C_"+year_start+"-"+year_end+"_monthly.nc"
29 c = addfile(diro+filo,"c")
30 filedimdef(c,"year",-1,True)
33 nyear = year_end - year_start + 1
36 ; data1 = new((/nyear,nmonth,nlat,nlon/),float)
37 data2 = new((/nyear,nmonth,nlat,nlon/),float)
38 data3 = new((/nyear,nmonth,nlat,nlon/),float)
39 date = new((/nyear,nmonth/),integer)
46 text = "mss:/FORREST/CCSM/results/"+model_name+"/lnd/hist/*"+year+"* "
49 system("msrcp "+text+diri)
55 date(m,n) = year*100 +nn
59 fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
61 fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
66 g = addfile (diri+fili,"r")
68 ; b1 = g->COL_FIRE_CLOSS
69 ; data1(m,n,:,:) = (/b1(0,:,:)/)
72 data2(m,n,:,:) = (/b2(0,:,:)/)
75 data3(m,n,:,:) = (/b3(0,:,:)/)
79 ; system("rm "+diri+"*")
83 ; name dimensions and assign coordinate variables
89 ; data1@long_name = b1@long_name
90 ; data1@units = b1@units
91 ; data1@_FillValue = b1@_FillValue
97 data2@long_name = b2@long_name
98 data2@units = b2@units
99 data2@_FillValue = b2@_FillValue
105 data3@long_name = b3@long_name
106 data3@units = b3@units
107 data3@_FillValue = b3@_FillValue
111 date@long_name = "date: yyyymm"
113 ; c->COL_FIRE_CLOSS = data1