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 ;************************************************
16 diri = "/fis/cgd/cseg/people/jeff/clamp_data/observed/biomass/"
17 fili = "amazon_mask_"+ grid + ".nc"
18 f = addfile (diri+fili,"r")
24 wks = gsn_open_wks ("ps","test") ; open workstation
25 gsn_define_colormap(wks,"gui_default") ; choose colormap
27 ;-----------------------------------------
31 gRes@txFontHeightF = 0.02
36 area_sum_text = "(mask area = "+sprintf("%.2f", Sum_area)+"(x1.e12 m2))"
38 gsn_text_ndc(wks,area_sum_text,0.50,0.80,gRes)
39 ;-----------------------------------------
41 ;************************************************
43 ;************************************************
45 res = True ; Use plot options
46 res@cnFillOn = True ; Turn on color fill
47 res@gsnSpreadColors = True ; use full colormap
48 ; res@cnFillMode = "RasterFill" ; Turn on raster color
49 ; res@lbLabelAutoStride = True
50 res@cnLinesOn = False ; Turn off contourn lines
51 res@mpFillOn = False ; Turn off map fill
52 res@gsnAddCyclic = False
53 res@tiMainString = "Amazon Mask: grid = "+ grid
55 res@mpMinLatF = -21.1 ; range to zoom in on
57 res@mpMinLonF = 277.28
58 res@mpMaxLonF = 326.43
60 res@gsnSpreadColors = True ; use full colormap
61 res@cnLevelSelectionMode = "ManualLevels" ; Manual contour invtervals
62 ; res@cnMinLevelValF = 1. ; Min level
63 ; res@cnMaxLevelValF = 11. ; Max level
64 ; res@cnLevelSpacingF = 1. ; interval
65 res@cnMinLevelValF = 0. ; Min level
66 res@cnMaxLevelValF = 1. ; Max level
67 res@cnLevelSpacingF = 0.1 ; interval
69 ; pres = True ; panel plot mods desired
70 ; pres@gsnMaximize = True ; fill the page
72 plot = gsn_csm_contour_map_ce(wks,y,res) ; for observed
73 ; gsn_panel(wks,plot,(/1,1/),pres) ; create panel plot
74 ; system("convert xy.ps xy.png")