<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Dear NCL experts,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Here is my script. Can you please take a look and advice what happens? Thank you very much!<o:p></o:p></p>
<p class="MsoNormal">=============================================<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<o:p></o:p></p>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">begin<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> files1 = systemfunc("ls ./T2_IV_DJF_50_rs*.nc")<o:p></o:p></p>
<p class="MsoNormal"> a = addfiles(files1,"r")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> f = addfile("T2_IV_DJF_12to50_16m.nc","r")<o:p></o:p></p>
<p class="MsoNormal"> c = addfile("/global/cscratch1/sd/wangjl/RCP85_IV_RUNS/50km_1995_IV_runs/wrfout_d01_1996-03-01_00_00_00"+".nc","r")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> lat = c->XLAT(0,:,:)<o:p></o:p></p>
<p class="MsoNormal"> lon = c->XLONG(0,:,:)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> ListSetType (a, "join")<o:p></o:p></p>
<p class="MsoNormal"> IV_10 = a[:]->IV<o:p></o:p></p>
<p class="MsoNormal"> IV_16 = f->IV<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> rain = new(dimsizes(IV_10),typeof(IV_10))<o:p></o:p></p>
<p class="MsoNormal">do r = 0, 199<o:p></o:p></p>
<p class="MsoNormal"> rain(r,:,:) = IV_10(r,:,:) - IV_16(:,:)<o:p></o:p></p>
<p class="MsoNormal">end do<o:p></o:p></p>
<p class="MsoNormal"> printVarSummary (rain)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> std_r = dim_stddev_n(rain,0)<o:p></o:p></p>
<p class="MsoNormal"> avg_r = dim_avg_n(rain,0)<o:p></o:p></p>
<p class="MsoNormal"> alpha = abs(avg_r/(std_r+0.00001))<o:p></o:p></p>
<p class="MsoNormal">; printMinMax(alpha,0)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">ff = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")<o:p></o:p></p>
<p class="MsoNormal">lsdata = ff->LSMASK<o:p></o:p></p>
<p class="MsoNormal">lsm = landsea_mask(lsdata,lat,lon)<o:p></o:p></p>
<p class="MsoNormal">avg_r = mask(avg_r, lsm .eq.0 .or. lsm.eq.2, False)<o:p></o:p></p>
<p class="MsoNormal">alpha = mask(alpha, lsm .eq.0 .or. lsm.eq.2, False)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;contour plot<o:p></o:p></p>
<p class="MsoNormal"> type = "ps"<o:p></o:p></p>
<p class="MsoNormal"> wks = gsn_open_wks(type,"T2_IV_DJF_50_vs_12_1995") ;;;no 1983 for sum, aut, win<o:p></o:p></p>
<p class="MsoNormal"> gsn_define_colormap(wks,"ViBlGrWhYeOrRe")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> res = True<o:p></o:p></p>
<p class="MsoNormal"> pltres = True<o:p></o:p></p>
<p class="MsoNormal"> mpres = True<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> mpres@mpGeophysicalLineColor = "Black"<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpNationalLineColor = "Black"<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpUSStateLineColor = "Black"<o:p></o:p></p>
<p class="MsoNormal"> mpres@mpGridAndLimbOn = False<o:p></o:p></p>
<p class="MsoNormal"> mpres@tiMainString = "IV 50 vs. 12"<o:p></o:p></p>
<p class="MsoNormal"> mpres@tiMainFontHeightF = 0.015<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> ; Plotting options for Precipitation<o:p></o:p></p>
<p class="MsoNormal"> opts_r = res<o:p></o:p></p>
<p class="MsoNormal"> opts_r@cnFillMode = "RasterFill"<o:p></o:p></p>
<p class="MsoNormal">opts_r@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels<o:p></o:p></p>
<p class="MsoNormal">opts_r@cnMinLevelValF = -3. ; set min contour level<o:p></o:p></p>
<p class="MsoNormal">opts_r@cnMaxLevelValF = 3. ; set max contour level<o:p></o:p></p>
<p class="MsoNormal">opts_r@cnLevelSpacingF = .3 ; set contour spacing<o:p></o:p></p>
<p class="MsoNormal"> opts_r@cnInfoLabelOn = False<o:p></o:p></p>
<p class="MsoNormal"> opts_r@cnConstFLabelOn = False<o:p></o:p></p>
<p class="MsoNormal"> opts_r@cnFillOn = True<o:p></o:p></p>
<p class="MsoNormal"> opts_r@lbLabelFont = 25<o:p></o:p></p>
<p class="MsoNormal"> opts_r@pmLabelBarOrthogonalPosF= -0.012<o:p></o:p></p>
<p class="MsoNormal"> contour_tot = wrf_contour(c,wks,avg_r,opts_r)<o:p></o:p></p>
<p class="MsoNormal"> delete(opts_r)<o:p></o:p></p>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
<p class="MsoNormal">; Plotting options for Significant level<o:p></o:p></p>
<p class="MsoNormal"> opts = res<o:p></o:p></p>
<p class="MsoNormal"> opts@cnFillOn = False<o:p></o:p></p>
<p class="MsoNormal"> opts@cnLinesOn = False<o:p></o:p></p>
<p class="MsoNormal"> opts@cnLineLabelsOn = False<o:p></o:p></p>
<p class="MsoNormal"> contour_sgl = wrf_contour(c,wks,alpha,opts)<o:p></o:p></p>
<p class="MsoNormal"> contour_sgl = ShadeGtContour(contour_sgl,4.,17)<o:p></o:p></p>
<p class="MsoNormal"> delete(opts)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> ; MAKE PLOTS<o:p></o:p></p>
<p class="MsoNormal"> plot = wrf_map_overlays(c,wks,(/contour_tot,contour_sgl/),pltres,mpres)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">end<o:p></o:p></p>
<p class="MsoNormal">=======================<o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Jiali Wang, Ph. D<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Environmental Science Division<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Argonne National Laboratory<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif">Tel: 630-252-2848<o:p></o:p></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">http://www.evs.anl.gov/about-evs/staff/detail/index.cfm?/Wang/Jiali</a><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>