1 ;************************************************
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)
47 ; text = "mss:/FORREST/CCSM/results/"+model_name+"/lnd/hist/*"+year+"* "
50 ; system("msrcp "+text+diri)
56 date(m,n) = year*100 +nn
60 fili = model_name+".clm2.h0."+year+"-"+nn+".nc"
62 fili = model_name+".clm2.h0."+year+"-0"+nn+".nc"
67 g = addfile (diri+fili,"r")
69 ; b1 = g->COL_FIRE_CLOSS
70 ; data1(m,n,:,:) = (/b1(0,:,:)/)
73 data2(m,n,:,:) = (/b2(0,:,:)/)
76 data3(m,n,:,:) = (/b3(0,:,:)/)
80 ; system("rm "+diri+"*")
84 ; name dimensions and assign coordinate variables
90 ; data1@long_name = b1@long_name
91 ; data1@units = b1@units
92 ; data1@_FillValue = b1@_FillValue
98 data2@long_name = b2@long_name
99 data2@units = b2@units
100 data2@_FillValue = b2@_FillValue
106 data3@long_name = b3@long_name
107 data3@units = b3@units
108 data3@_FillValue = b3@_FillValue
112 date@long_name = "date: yyyymm"
114 ; c->COL_FIRE_CLOSS = data1