forrest@0: ;*************************************************
forrest@0: ; ce_1.ncl
forrest@0: ;************************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
forrest@0: ;************************************************
forrest@0: begin
forrest@0: ;************************************************
forrest@0: ; read in observed data
forrest@0: ;************************************************
forrest@0:   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/"
forrest@0:   fili  = "LAI_2000-2005_mean_T42.nc"
forrest@0:   f     = addfile(diri+fili,"r")
forrest@0: 
forrest@0:   y     = f->LAI
forrest@0:   printVarSummary(y)
forrest@0: 
forrest@0:   print (min(y)+"/"+max(y))
forrest@0: 
forrest@0: ;************************************************
forrest@0: ; read in model data
forrest@0: ;************************************************
forrest@0:   fili2 = "i01.03cn_1545-1569_ANN_climo.nc"
forrest@0:   g     = addfile(diri+fili2,"r")
forrest@0:   x     = g->LAISUN
forrest@0: 
forrest@0:   delta = 0.0000001
forrest@0:   x0    = x(0,:,:)
forrest@0:   y = where(ismissing(y).and.(ismissing(x0).or.(x0.lt.delta)),0.,y)
forrest@0: ; y = where((y.gt.1000.),1.e+36,y)
forrest@0:   printVarSummary(y)
forrest@0:   print (min(y)+"/"+max(y))
forrest@0: ;************************************************
forrest@0: ; create default plot
forrest@0: ;************************************************
forrest@0:   
forrest@0: ; setvalues NhlGetWorkspaceObjectId()
forrest@0: ;   "wsMaximumSize" : 199999999
forrest@0: ; end setvalues
forrest@0: 
forrest@0:   wks = gsn_open_wks("ps","xy")         ; open a ps file
forrest@0: ; gsn_define_colormap(wks,"wgne15")          ; choose colormap
forrest@0:   gsn_define_colormap(wks,"gui_default")     ; choose colormap
forrest@0: 
forrest@0:   res                     = True             ; Use plot options
forrest@0:   res@cnFillOn            = True             ; Turn on color fill
forrest@0:   res@gsnSpreadColors      = True            ; use full colormap
forrest@0: ; res@cnFillMode          = "RasterFill"     ; Turn on raster color
forrest@0: ; res@lbLabelAutoStride   = True
forrest@0:   res@cnLinesOn           = False            ; Turn off contourn lines
forrest@0:   res@mpFillOn            = False            ; Turn off map fill
forrest@0:   res@tiMainString        = "MODIS MOD 15A2 2000-2005"
forrest@0: 
forrest@0:   res@gsnSpreadColors     = True             ; use full colormap
forrest@0:   res@cnLevelSelectionMode = "ManualLevels"  ; Manual contour invtervals
forrest@0:   res@cnMinLevelValF       = 0.              ; Min level
forrest@0:   res@cnMaxLevelValF       = 10.             ; Max level
forrest@0:   res@cnLevelSpacingF      = 1.              ; interval
forrest@0:   pres                            = True     ; panel plot mods desired
forrest@0:   pres@gsnMaximize                = True     ; fill the page
forrest@0: 
forrest@0:   plot = gsn_csm_contour_map_ce(wks,y,res)   ; for observed
forrest@0:   gsn_panel(wks,plot,(/1,1/),pres)           ; create panel plot
forrest@0:   system("convert xy.ps xy.png")
forrest@0: end