<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none"><!-- P { margin-top: 0px; margin-bottom: 0px; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Thanks, Mary and Rick. I expected some missing values where large water bodies existed, but the zoom-in Mary provided helped diagnose the issue. I had a typo in the SnowDays function that looped over every other lat/lon instead of every lat/lon pair,
lines 43-44 in subsetSNODAS_avgSNWDAYS.ncl should have been:<br>
</p>
<p><br>
</p>
<p>do nl=0,nlat-1,1<br>
</p>
<p> do ml=0,mlat-1,1<br>
</p>
<p><br>
</p>
<p>instead of the original:<br>
</p>
<p><br>
</p>
<p style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 16px; background-color: rgb(255, 255, 255);">
do nl=0,nlat-1,2<br>
</p>
<p style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 16px; background-color: rgb(255, 255, 255);">
do ml=0,mlat-1,2<br>
</p>
<p><br>
</p>
<p>The plot now looks as I would expect, mottling is gone. Thanks!<br>
</p>
<p><br>
</p>
<p>Liz.<br>
</p>
<p><br>
</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font face="Corbel, sans-serif">Dr. Elizabeth A Burakowski</font>
<div><font face="Corbel, sans-serif">Research Assistant Professor</font></div>
<div><font face="Corbel, sans-serif">Earth Systems Research Center & </font></div>
<div><font face="Corbel, sans-serif">Institute for the Study of Earth, Oceans, and Space</font></div>
<div><font face="Corbel, sans-serif">University of New Hampshire</font></div>
<div><font face="Corbel, sans-serif">8 College Road Durham, NH 03824</font></div>
<div><font face="Corbel, sans-serif">603-862-1796</font></div>
<div><font face="Corbel, sans-serif">elizabeth.burakowski@unh.edu</font></div>
</div>
</div>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Mary Haley <haley@ucar.edu><br>
<b>Sent:</b> Thursday, August 17, 2017 12:58 PM<br>
<b>To:</b> Rick Brownrigg<br>
<b>Cc:</b> Burakowski, Elizabeth; ncl-talk@ucar.edu<br>
<b>Subject:</b> Re: [ncl-talk] Mottled output from gsn_csm_contour_map</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hi Elizabeth,</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">I was working on a plot to show the very issue that Rick pointed out. He's right that the data is pretty sparse.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Here's a modified version of your script that zooms in on the map a little, and draws filled dots in red where the data is missing, and in black otherwise. You can see that each valid point is surrounded by
missing points.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">I had to zoom in on the map because the data is so dense.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">--Mary</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Aug 17, 2017 at 10:31 AM, Rick Brownrigg <span dir="ltr">
<<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="ltr">
<div>
<div>
<div>Hi,<br>
<br>
</div>
I don't fully understand how the data are being computed, but the output file appears to be rather sparse, and that is likely the cause of the mottled appearance. Attached is a "CellFill" version of the plot -- CellFill more nearly shows the grid cells.<br>
<br>
And If I print sections of the ShowDays variable, its very typically like below; lots of missing-values interspersed with actual values.<br>
<br>
</div>
I'm not much help here, but hopefully this provides some insight...<br>
</div>
Rick<br>
<div>
<div>
<div><br>
285,250) -9999<br>
(286,0) 6.2<br>
(286,1) -9999<br>
(286,2) 7.4<br>
(286,3) -9999<br>
(286,4) 8<br>
(286,5) -9999<br>
(286,6) 7.9<br>
(286,7) -9999<br>
(286,8) 10.22222<br>
(286,9) -9999<br>
(286,10) 8.7<br>
(286,11) -9999<br>
(286,12) 9<br>
(286,13) -9999<br>
(286,14) 8<br>
(286,15) -9999<br>
(286,16) 10.1<br>
(286,17) -9999<br>
(286,18) 8.4<br>
(286,19) -9999<br>
(286,20) 11<br>
(286,21) -9999<br>
(286,22) 8.454545<br>
(286,23) -9999<br>
(286,24) 11.5<br>
(286,25) -9999<br>
(286,26) 9.363636<br>
(286,27) -9999<br>
(286,28) 9<br>
(286,29) -9999<br>
(286,30) 9.5<br>
(286,31) -9999<br>
(286,32) 9.272727<br>
(286,33) -9999<br>
(286,34) 9.3<br>
(286,35) -9999<br>
(286,36) 8.5<br>
(286,37) -9999<br>
(286,38) 8.2<br>
(286,39) -9999<br>
(286,40) 7.636364<br>
(286,41) -9999<br>
<br>
</div>
</div>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">
<div>
<div class="h5">On Thu, Aug 17, 2017 at 9:37 AM, Burakowski, Elizabeth <span dir="ltr">
<<a href="mailto:Elizabeth.Burakowski@unh.edu" target="_blank">Elizabeth.Burakowski@unh.edu</a>></span> wrote:<br>
</div>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div>
<div class="h5">
<div dir="ltr" style="font-size:12pt; color:#000000; background-color:#ffffff; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi ncl-talk, <br>
</p>
<p><br>
</p>
<p>I am trying to create a map showing the number of days with snow depth greater than 200 mm from SNODAS files. I have two scripts attached:</p>
<p><br>
</p>
<p>(1) subsetSNODAS_avgSNWDAYS.ncl<br>
</p>
<p><span style="font-size:12pt">- subsets a smaller region from nation SNODAS domain, </span><span style="font-size:12pt"> </span></p>
<p><span style="font-size:12pt">- calculates the number of days with snow depth greater than user defined threshold</span></p>
<p><span style="font-size:12pt">- calculates 2004-2017 average number of days with snow depth greater than threshold</span></p>
<p><span style="font-size:12pt">- creates .nc file with 2004-2017 average</span></p>
<p><span style="font-size:12pt"><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">I have uploaded a few .nc source SNODAS files that serve as input to the subset scrip</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">t
and the resulting 2004-2017 .nc file to the n</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; background-color:rgb(255,255,255)">cl-talk ftp site.</span><br>
</span></p>
<p><br>
</p>
<p>(2) plotNESNODAS_snowdays.ncl<br>
</p>
<p>- plots the output .nc file produced in subsetSNODAS_avgSNWDAYS.ncl<br>
</p>
<p>- this is a separate script because the subset script takes a while to run on the full time series for large regions. <br>
</p>
<p><br>
</p>
<p>Thanks, <br>
</p>
<p><br>
</p>
<p>Liz.<br>
</p>
<p><br>
</p>
<p><br>
</p>
<div id="m_-7870137323674382634m_-7333007136264988580Signature">
<div name="divtagdefaultwrapper"><font face="Corbel, sans-serif">Dr. Elizabeth A Burakowski</font>
<div><font face="Corbel, sans-serif">Research Assistant Professor</font></div>
<div><font face="Corbel, sans-serif">Earth Systems Research Center & </font></div>
<div><font face="Corbel, sans-serif">Institute for the Study of Earth, Oceans, and Space</font></div>
<div><font face="Corbel, sans-serif">University of New Hampshire</font></div>
<div><font face="Corbel, sans-serif">8 College Road Durham, NH 03824</font></div>
<div><font face="Corbel, sans-serif"><a href="tel:(603)%20862-1796" value="+16038621796" target="_blank">603-862-1796</a></font></div>
<div><font face="Corbel, sans-serif"><a href="mailto:elizabeth.burakowski@unh.edu" target="_blank">elizabeth.burakowski@unh.edu</a></font></div>
</div>
</div>
</div>
<br>
</div>
</div>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=c6MrceVCY5m5A_KAUkrdoA&r=JW2rHHzYnVYi6ZDBI-2JFMNSreDPLZO-wwOrczEnOlg&m=mbqQACvpKD_SdEZEWAtkUeePV-5K9I1V6PKH-xNjRwU&s=AIY6vClsvgpiKhjoJgxyskPzPhzWQCrSMMnljn9s23s&e=" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</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="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=c6MrceVCY5m5A_KAUkrdoA&r=JW2rHHzYnVYi6ZDBI-2JFMNSreDPLZO-wwOrczEnOlg&m=mbqQACvpKD_SdEZEWAtkUeePV-5K9I1V6PKH-xNjRwU&s=AIY6vClsvgpiKhjoJgxyskPzPhzWQCrSMMnljn9s23s&e=" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</body>
</html>