[ncl-talk] Fwd: Problem with contour map shading

Mary Haley haley at ucar.edu
Fri Apr 21 09:45:50 MDT 2017


Perhaps you can fix this with a "draw order" resource. Try:

res at mpFillDrawOrder = "PostDraw"



On Thu, Apr 20, 2017 at 11:21 PM, Priyanka Silva <
priyankakeraminiya6 at gmail.com> wrote:

> Still I have the same problem.
>
>
> ---------- Forwarded message ----------
> From: Barry Lynn <barry.h.lynn at gmail.com>
> Date: Thu, Apr 20, 2017 at 4:59 PM
> Subject: Re: [ncl-talk] Problem with contour map shading
> To: Priyanka Silva <priyankakeraminiya6 at gmail.com>
> Cc: ncl-talk <ncl-talk at ucar.edu>
>
>
> Hi Priyanka:
>
> Try the "mask" routine.
>
> range_only  = snowh_rate_dom                               ; trick to keep
> cv's and atts
>
>   range_only  = mask(snowh_rate_dom,(snowh_rate_dom.ge.0.1),True)
>
> I think you can add a conditional about whether you are over land and/or
> ocean.
>
> Barry
>
> On Thu, Apr 20, 2017 at 10:26 AM, Priyanka Silva <
> priyankakeraminiya6 at gmail.com> wrote:
>
>> Hi all,
>>
>> I want to shade the areas exceeding 95% confidence level in my
>> correlation plot only for the ocean, but when I plot I could see the
>> shading on the land too. How can I remove the shading on the land?
>>
>>  My script and plot are as follows.
>>
>>
>>
>> Greatly appreciate your help.
>>
>> Thank you.
>>
>>
>>
>> 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/contributed.ncl"
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
>>
>>
>>
>> begin
>>
>> ;***********************************************************
>> ********************
>>
>>     ;Seasonal mean precipitation averaged over the domain
>>
>> ;***********************************************************
>> ******************
>>
>> f    = "mon.precip.nc"
>>
>> fin  = addfile(f, "r")
>>
>> time = fin->time
>>
>> ymd = cd_calendar(time,2)
>>
>> iTim = ind (ymd.ge.19790101.and.ymd.lt.20130101)
>>
>> preci = fin->precip(iTim,{6:10},{79:82})
>>
>> ;printVarSummary(precip)
>>
>> ;-----------------------------------------------------------
>> --------------
>>
>> p = month_to_season(preci,"JJA")
>>
>> printVarSummary(p)
>>
>> ;-----------------------------------------------------------
>> ---------------
>>
>> index =dim_avg_Wrap (dim_avg_Wrap (p (time|:,lat|:,lon|:)))
>>
>> printVarSummary(index)
>>
>> ;************************************************************************
>>
>>   ; correlation between the precipitation and the global SST
>>
>> ;************************************************************************
>>
>> f1 = "HadISST_sst.nc"
>>
>> fin1 = addfile(f1, "r")
>>
>> time1 = fin1->time
>>
>> ymd1 = cd_calendar(time1,2)
>>
>> iTim1 = ind (ymd1.ge.19790101.and.ymd1.lt.20130101)
>>
>> sst = fin1->sst(iTim1,:,:)
>>
>> ;;----------------------------------------------------------
>> --------------
>>
>> s = month_to_season(sst,"JJA")
>>
>> ;printVarSummary(s)
>>
>> ;-----------------------------------------------------------
>> --------------
>>
>> ccr = escorc(index,s(latitude|:,longitude|:,time|:))
>>
>> copy_VarCoords(sst(0,:,:),ccr)
>>
>> printVarSummary(ccr)
>>
>> ;------------------------------------------------------------------------
>>
>> siglvl = 0.05
>>
>> Prob = rtest(ccr(latitude|:,longitude|:),34,0)
>>
>>
>>
>> Prob!0 = "lat"
>>
>> Prob!1 = "lon"
>>
>> Prob&lat = ccr&latitude
>>
>> Prob&lon = ccr&longitude
>>
>> printVarSummary(Prob)
>>
>> ;-----------------------------------------------------------
>> --------------
>>
>> wks  = gsn_open_wks("x11","Correlation")
>>
>> gsn_define_colormap(wks,"BlWhRe")
>>
>> ;-----------------------------------------------------------
>> --------------
>>
>> res                      = True
>>
>> res at cnFillOn             = True
>>
>> res at cnLinesOn            = False
>>
>> res at cnLevelSelectionMode = "ManualLevels"
>>
>> res at cnMinLevelValF       = -.8
>>
>> res at cnMaxLevelValF       =  .8
>>
>> res at cnLevelSpacingF      = 0.1
>>
>> res at gsnSpreadColors = True
>>
>> res at mpPerimOn = False
>>
>> res at mpCenterLonF =180
>>
>> res at tiMainString = ""
>>
>> res at gsnLeftString         = "JJA"
>>
>> res at mpFillOn  = False
>>
>> res at gsnDraw = False
>>
>> res at gsnFrame = False
>>
>> ;------------------------------------------------------------------------
>>
>>
>>
>> res1 = True
>>
>> res1 at gsnAddCyclic = False
>>
>> res1 at gsnDraw = False
>>
>> res1 at gsnFrame = False
>>
>> res1 at cnLevelSelectionMode = "ManualLevels"
>>
>> res1 at cnMinLevelValF = 0.00
>>
>> res1 at cnMaxLevelValF = 1.05
>>
>> res1 at cnLevelSpacingF = 0.05
>>
>> res1 at cnInfoLabelOn = False
>>
>> res1 at cnLinesOn = False
>>
>> res1 at cnLineLabelsOn = False
>>
>> res1 at cnFillScaleF = 0.6 ; add extra density
>>
>> res1 at mpLandFillColor ="gray"
>>
>> res1 at gsnShadeDotSizeF =3
>>
>> ;-----------------------------------------------------------
>> ------------------
>>
>> plot = gsn_csm_contour_map(wks,ccr,res)
>>
>>
>>
>> opt = True
>>
>> opt at gsnShadeFillType = "pattern"       ; pattern fill
>>
>> opt at gsnShadeLow = 17                   ; use pattern 17
>>
>>
>>
>> plot1 = gsn_csm_contour(wks,Prob, res1)
>>
>> plot1 = gsn_contour_shade(plot1,0.05,0.75,opt)
>>
>>
>>
>> overlay (plot, plot1)
>>
>> draw(plot)
>>
>> frame (wks)
>>
>>
>>
>> end
>>
>>
>>
>>  [image: Inline image 1]
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108 <(914)%20432-3108>
>
>
> _______________________________________________
> 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/20170421/e60c3c41/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 97171 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170421/e60c3c41/attachment.png 


More information about the ncl-talk mailing list