Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;*************************************************
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 diro = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
13 filo = "class_pft_T42.nc"
14 c = addfile(diro+filo,"c")
18 diri = "/fis/cgd/cseg/people/jeff/surface_data/"
19 fili = "T42.clamp.surface-data.060412.nc"
20 f = addfile(diri+fili,"r")
25 ; read model grid data
27 diri = "/fis/cgd/cseg/people/jeff/surface_data/"
29 g = addfile(diri+fili,"r")
31 landmask = g->landmask
45 x@long_name = "Model PFT Classes"
49 x(j,i) = maxind(y(:,j,i))
55 x = where(landmask .lt. 1.,x@_FillValue,x)