Hi all,
I am trying to shade an area over a colored contour plot (lat/lon).
I have succeed to reproduce the shading over the map without a color
table (without the line " res_at_cnFillOn = True ", see the script
bellow).
But when I add this line, I obtain a map in Black and White with the
contour + the shading + the color table BUT no colored map.
I don't figure out what is the problem in that script? Is the function
Shade appropriate for my plot?
Your help is much appreciated.
Thanks,
Gwenaelle
;***********************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
f = addfile("toto.nc","r")
t = f->TS(0,:,:)
t = t -273.15
wks = gsn_open_wks ("ps", "toto" ) ; open workstation
res = True ; plot mods desired
res_at_gsnDraw = False ; do not draw
res_at_gsnFrame = False ; do not advance
frame
res_at_cnFillOn = True ; Turn on color
plot = gsn_csm_contour_map(wks, t, res )
plot = ShadeLtContour(plot,10.,17)
draw(plot)
frame(wks)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 10 2009 - 15:36:02 MDT
This archive was generated by hypermail 2.2.0 : Thu Mar 12 2009 - 14:10:21 MDT