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 ;************************************************
8 ;************************************************
10 ;************************************************
11 diri = "/fis/cgd/cseg/people/jeff/clamp_data/"
12 fili = "i01.03cn_1545-1569_ANN_climo.nc"
13 ; fili = "i01.04casa_1605-1629_ANN_climo.nc"
15 g = addfile(diri+fili,"r")
20 ;************************************************
22 ;************************************************
24 wks = gsn_open_wks("ps","xy") ; open a ps file
25 ; gsn_define_colormap(wks,"wgne15") ; choose colormap
26 gsn_define_colormap(wks,"gui_default") ; choose colormap
28 res = True ; Use plot options
29 res@cnFillOn = True ; Turn on color fill
30 res@gsnSpreadColors = True ; use full colormap
31 ; res@cnFillMode = "RasterFill" ; Turn on raster color
32 ; res@lbLabelAutoStride = True
33 res@cnLinesOn = False ; Turn off contourn lines
34 res@mpFillOn = False ; Turn off map fill
35 ; res@tiMainString = "Observed MODIS MOD 15A2"
36 res@tiMainString = "Model i01.03cn"
37 ; res@tiMainString = "Model i01.04casa"
39 res@gsnSpreadColors = True ; use full colormap
40 res@cnLevelSelectionMode = "ManualLevels" ; Manual contour invtervals
41 res@cnMinLevelValF = 0. ; Min level
42 res@cnMaxLevelValF = 10. ; Max level
43 res@cnLevelSpacingF = 1. ; interval
45 plot = gsn_csm_contour_map_ce(wks,x(0,:,:),res) ; for model
46 system("convert xy.ps xy.png")