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 ;************************************************
10 ;************************************************
12 ;************************************************
13 diri = "/fis/cgd/cseg/people/jeff/clamp/class/"
14 fili = "T42.clamp.surface-data.060412.nc"
15 f = addfile(diri+fili,"r")
20 diri = "/fis/cgd/cseg/people/jeff/surface_data/"
22 g = addfile(diri+fili,"r")
24 landmask = g->landmask
38 x@long_name = "Model PFT Classes"
42 x(j,i) = maxind(y(:,j,i))
48 x = where(landmask .lt. 1.,x@_FillValue,x)
50 ;************************************************
52 ;************************************************
54 wks = gsn_open_wks("ps","class") ; open a ps file
55 gsn_define_colormap(wks,"gui_default") ; choose colormap
57 res = True ; Use plot options
58 res@cnFillOn = True ; Turn on color fill
59 res@gsnSpreadColors = True ; use full colormap
60 res@cnLinesOn = False ; Turn off contourn lines
61 res@mpFillOn = False ; Turn off map fill
62 res@tiMainString = "Class"
64 res@gsnSpreadColors = True ; use full colormap
65 res@cnLevelSelectionMode = "ManualLevels" ; Manual contour invtervals
66 res@cnMinLevelValF = 0. ; Min level
67 res@cnMaxLevelValF = 16. ; Max level
68 res@cnLevelSpacingF = 1. ; interval
70 plot = gsn_csm_contour_map_ce(wks,x,res)