[ncl-talk] gsn_contour_shade

Mary Haley haley at ucar.edu
Tue Aug 11 10:18:49 MDT 2015


Laura,

I think the problem is simply that the min/max of your data is 0.0111 and
0.69999. You are telling it to shade all values above 0.7, but there are no
values above 0.7, so you end up with no shading. If you use a smaller value
like 0.6, you will get some yellow contours.

--Mary



On Fri, Aug 7, 2015 at 1:57 PM, Laura Fowler <laura at ucar.edu> wrote:

> Hello:
>
> I am trying to use gsn_contour_shade to plot a variable that varies
> betwen 0.1 and 0.7, following the bottom of the script coneff_13.ncl.
> Unfortunately, I cannot make the area of values greater than 0.7
> filled with color 150 and I am not sure why. My script is in
>
> /glade/p/work/laura/GFconvection/MPAS.GFconvection-v3.3.mynn/rotation.50_3.overSouthAmerica
>
> but basically, it contains the lines:
>
> ;SET RESOURCES:
>
> res = True
> res at gsnDraw             = False
> res at gsnFrame            = False
> res at gsnMaximize         = False
> res at gsnSpreadColors     = False
> res at gsnPaperOrientation = "landscape"
>
> res at mpProjection        = "CylindricalEquidistant"
> res at mpDataBaseVersion   = "MediumRes"
> res at mpCenterLatF        = 0.
> res at mpCenterLonF        = 0.
> res at mpMinLatF           =  -90.
> res at mpMaxLatF           =   90.
> res at mpMinLonF           = -150.
> res at mpMaxLonF           =   30.
>
> diri = "./"
> f = addfile("50_3.grid.overSouthAmerica.nc","r")
>
> .
> .
> .
> .
>
> res at cnLevelSelectionMode        = "ManualLevels"
> res at cnMinLevelValF              = 0.1
> res at cnMaxLevelValF              = 0.7
> res at cnLevelSpacingF             = 0.1
>
> ;res at mpOceanFillColor       = "lightcyan"
> ;res at mpInlandWaterFillColor = "lightcyan"
> ;res at mpLandFillColor        = "gray75"
>
> wks = gsn_open_wks("pdf","schematicRefineMesh.limited")
> gsn_define_colormap(wks,"wh-bl-gr-ye-re")
>
> plot = gsn_csm_contour_map(wks,ConvF_meshD,res)
>
> opt = True
> opt at gsnShadeHigh = 150
> plot = gsn_contour_shade(plot,ConvF_meshD at _FillValue,0.7,opt)
>
> draw(plot)
> frame(wks)
>
>
> A few hints would be greatly appreciated. Thanks,
> Laura
>
>
> --
>
> !-------------------------------------------------------------------------------------------------------------
> Laura D. Fowler
> Mesoscale and Microscale Meteorology Division (MMM)
> National Center for Atmospheric Research
> P.O. Box 3000, Boulder CO 80307-3000
>
> e-mail: laura at ucar.edu
> phone: 303-497-1628
>
>
> !-------------------------------------------------------------------------------------------------------------
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150811/0659e2cb/attachment.html 


More information about the ncl-talk mailing list