<div dir="ltr">Hi Jiali,<div>ShadeGtContour will shade everything starting with the first contour less than the given input contour (in your case, 4). However, as you are not setting the contours and letting NCL select them, NCL is choosing contours to go from 40 to 400 by 40. Thus, I think the solution is to set the contour levels yourself when you call <span style="font-size:12.8px">wrf_contour. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Try this:</span></div><div><p class="MsoNormal" style="font-size:12.8px">opts = res<u></u><u></u></p><p class="MsoNormal" style="font-size:12.8px"> opts@cnFillOn = False<u></u><u></u></p><p class="MsoNormal" style="font-size:12.8px"> opts@cnLinesOn = False<u></u><u></u></p><p class="MsoNormal" style="font-size:12.8px"> opts@cnLineLabelsOn = False<u></u><u></u></p><div><span style="font-size:12.8px">opts@cnLevelSelectionMode = "ExplicitLevels"</span><br></div><div><span style="font-size:12.8px">opts@cnLevels = (/2.,4./) ; note that a contour level of 4 is not really necessary</span></div><p class="MsoNormal" style="font-size:12.8px"> contour_sgl = wrf_contour(c,wks,alpha,opts)<u></u><u></u></p><p class="MsoNormal" style="font-size:12.8px"> contour_sgl = ShadeGtContour(contour_sgl,2.1,<wbr>17)<u></u><u></u></p><p class="MsoNormal" style="font-size:12.8px"> delete(opts)</p><p class="MsoNormal" style="font-size:12.8px"><br></p><p class="MsoNormal" style="font-size:12.8px">If that does not give you what you expect I would take a close look at your contour_sgl array (if you have not done so already) to make sure the values are what you think they are.</p><p class="MsoNormal" style="font-size:12.8px"><br></p><p class="MsoNormal" style="font-size:12.8px">Finally, note that ShadeGtContour is deprecated. It is suggested you use gsn_contour_<span style="font-size:12.8px">shade</span><span style="font-size:12.8px"> instead:</span></p><p class="MsoNormal"><span style="font-size:12.8px"><a href="https://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_contour_shade.shtml">https://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_contour_shade.shtml</a></span><br></p><p class="MsoNormal"><span style="font-size:12.8px"><br></span></p><p class="MsoNormal"><span style="font-size:12.8px">If you have any further questions please respond to the ncl-talk email list.</span></p><p class="MsoNormal"><span style="font-size:12.8px">Adam</span></p></div><div><span style="font-size:12.8px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 18, 2017 at 4:22 PM, Wang, Jiali <span dir="ltr"><<a href="mailto:jialiwang@anl.gov" target="_blank">jialiwang@anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="m_-722483701460875027WordSection1">
<p class="MsoNormal">Dear NCL experts,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am trying to add stipples (or hatchings, doesn’t matter) on a contour map when my “alpha >2”. I plot alpha, which is not always larger than 2, meaning that I should have some places where I don’t have stipples. However, I get stipples
everywhere even when I set the threshold to 4. Attached are two maps, one is alpha, and the other is something else but with sitpples when alpha>2.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Here is my script. Can you please take a look and advice what happens? Thank you very much!<u></u><u></u></p>
<p class="MsoNormal">==============================<wbr>===============<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/wrf/WRFUserARW.ncl"<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl"<u></u><u></u></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">begin<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> files1 = systemfunc("ls ./T2_IV_DJF_50_rs*.nc")<u></u><u></u></p>
<p class="MsoNormal"> a = addfiles(files1,"r")<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> f = addfile("T2_IV_DJF_12to50_16m.<wbr>nc","r")<u></u><u></u></p>
<p class="MsoNormal"> c = addfile("/global/cscratch1/sd/<wbr>wangjl/RCP85_IV_RUNS/50km_<wbr>1995_IV_runs/wrfout_d01_1996-<wbr>03-01_00_00_00"+".nc","r")<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> lat = c->XLAT(0,:,:)<u></u><u></u></p>
<p class="MsoNormal"> lon = c->XLONG(0,:,:)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> ListSetType (a, "join")<u></u><u></u></p>
<p class="MsoNormal"> IV_10 = a[:]->IV<u></u><u></u></p>
<p class="MsoNormal"> IV_16 = f->IV<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> rain = new(dimsizes(IV_10),typeof(IV_<wbr>10))<u></u><u></u></p>
<p class="MsoNormal">do r = 0, 199<u></u><u></u></p>
<p class="MsoNormal"> rain(r,:,:) = IV_10(r,:,:) - IV_16(:,:)<u></u><u></u></p>
<p class="MsoNormal">end do<u></u><u></u></p>
<p class="MsoNormal"> printVarSummary (rain)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> std_r = dim_stddev_n(rain,0)<u></u><u></u></p>
<p class="MsoNormal"> avg_r = dim_avg_n(rain,0)<u></u><u></u></p>
<p class="MsoNormal"> alpha = abs(avg_r/(std_r+0.00001))<u></u><u></u></p>
<p class="MsoNormal">; printMinMax(alpha,0)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">ff = addfile("$NCARG_ROOT/lib/<wbr>ncarg/data/cdf/<a href="http://landsea.nc" target="_blank">landsea.nc</a>","r"<wbr>)<u></u><u></u></p>
<p class="MsoNormal">lsdata = ff->LSMASK<u></u><u></u></p>
<p class="MsoNormal">lsm = landsea_mask(lsdata,lat,lon)<u></u><u></u></p>
<p class="MsoNormal">avg_r = mask(avg_r, lsm .eq.0 .or. lsm.eq.2, False)<u></u><u></u></p>
<p class="MsoNormal">alpha = mask(alpha, lsm .eq.0 .or. lsm.eq.2, False)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;contour plot<u></u><u></u></p>
<p class="MsoNormal"> type = "ps"<u></u><u></u></p>
<p class="MsoNormal"> wks = gsn_open_wks(type,"T2_IV_DJF_<wbr>50_vs_12_1995") ;;;no 1983 for sum, aut, win<u></u><u></u></p>
<p class="MsoNormal"> gsn_define_colormap(wks,"<wbr>ViBlGrWhYeOrRe")<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> res = True<u></u><u></u></p>
<p class="MsoNormal"> pltres = True<u></u><u></u></p>
<p class="MsoNormal"> mpres = True<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> mpres@mpGeophysicalLineColor = "Black"<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpNationalLineColor = "Black"<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpUSStateLineColor = "Black"<u></u><u></u></p>
<p class="MsoNormal"> mpres@mpGridAndLimbOn = False<u></u><u></u></p>
<p class="MsoNormal"> mpres@tiMainString = "IV 50 vs. 12"<u></u><u></u></p>
<p class="MsoNormal"> mpres@tiMainFontHeightF = 0.015<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> ; Plotting options for Precipitation<u></u><u></u></p>
<p class="MsoNormal"> opts_r = res<u></u><u></u></p>
<p class="MsoNormal"> opts_r@cnFillMode = "RasterFill"<u></u><u></u></p>
<p class="MsoNormal">opts_r@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels<u></u><u></u></p>
<p class="MsoNormal">opts_r@cnMinLevelValF = -3. ; set min contour level<u></u><u></u></p>
<p class="MsoNormal">opts_r@cnMaxLevelValF = 3. ; set max contour level<u></u><u></u></p>
<p class="MsoNormal">opts_r@cnLevelSpacingF = .3 ; set contour spacing<u></u><u></u></p>
<p class="MsoNormal"> opts_r@cnInfoLabelOn = False<u></u><u></u></p>
<p class="MsoNormal"> opts_r@cnConstFLabelOn = False<u></u><u></u></p>
<p class="MsoNormal"> opts_r@cnFillOn = True<u></u><u></u></p>
<p class="MsoNormal"> opts_r@lbLabelFont = 25<u></u><u></u></p>
<p class="MsoNormal"> opts_r@<wbr>pmLabelBarOrthogonalPosF= -0.012<u></u><u></u></p>
<p class="MsoNormal"> contour_tot = wrf_contour(c,wks,avg_r,opts_<wbr>r)<u></u><u></u></p>
<p class="MsoNormal"> delete(opts_r)<u></u><u></u></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;;;;;;<u></u><u></u></p>
<p class="MsoNormal">; Plotting options for Significant level<u></u><u></u></p>
<p class="MsoNormal"> opts = res<u></u><u></u></p>
<p class="MsoNormal"> opts@cnFillOn = False<u></u><u></u></p>
<p class="MsoNormal"> opts@cnLinesOn = False<u></u><u></u></p>
<p class="MsoNormal"> opts@cnLineLabelsOn = False<u></u><u></u></p>
<p class="MsoNormal"> contour_sgl = wrf_contour(c,wks,alpha,opts)<u></u><u></u></p>
<p class="MsoNormal"> contour_sgl = ShadeGtContour(contour_sgl,4.,<wbr>17)<u></u><u></u></p>
<p class="MsoNormal"> delete(opts)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> ; MAKE PLOTS<u></u><u></u></p>
<p class="MsoNormal"> plot = wrf_map_overlays(c,wks,(/<wbr>contour_tot,contour_sgl/),<wbr>pltres,mpres)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">end<u></u><u></u></p>
<p class="MsoNormal">=======================<u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Jiali Wang, Ph. D<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Environmental Science Division<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Argonne National Laboratory<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Tel: <a href="tel:(630)%20252-2848" value="+16302522848" target="_blank">630-252-2848</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif"><a href="http://www.evs.anl.gov/about-evs/staff/detail/index.cfm?/Wang/Jiali" target="_blank">http://www.evs.anl.gov/about-<wbr>evs/staff/detail/index.cfm?/<wbr>Wang/Jiali</a><u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>