<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Saji,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I was able to reproduce the problem with NCL V6.4.0, and verified that 6.5.0 (in testing right now) also has the issue.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">One thing I did note is that if you set:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">res@cnFillMode = "RasterFill"</div><div class="gmail_default" style="font-size:small">res@cnRasterSmoothingOn = True</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">you get a very similar plot and it works with 6.3.0 or 6.4.0. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Also, with the original smoothing, if you set cnSmoothingDistanceF to 0.004, then the plot is correct with NCL V6.4.0.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've attached a modified version of your original script that shows the results with raster smoothing and with setting different values for cnSmoothingDistanceF. Hopefully you can use one of these as a work-around, if you need to use NCL V6.4.0.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've reported this issue in ticket NCL-2759. <span style="background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Thanks for providing the files and the script.<span> </span></span></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 Mon, Jun 25, 2018 at 7:23 PM, Saji Hameed <span dir="ltr"><<a href="mailto:saji.nh@gmail.com" target="_blank">saji.nh@gmail.com</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">Hi Mary,<div><br></div><div>Thank you for checking this. I have uploaded two files (<a href="http://saji_rain_jun-sep.nc" target="_blank">saji_rain_jun-sep.nc</a>, <a href="http://saji_landsea_mask.nc" target="_blank">saji_landsea_mask.nc</a>) to <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a>.</div><div>The corresponding NCL script is below (I noted that the issue appears only when using the option cnSmoothingOn = True</div><div><br></div><div>Best,</div><div>saji</div><div><br></div><div><div>lt1 = 5</div><div>lt2 = 25</div><div>ln1 = 60</div><div>ln2 = 90</div><div><br></div><div>fin = addfile("<a href="http://saji_rain_jun-sep.nc" target="_blank">saji_rain_jun-sep.nc</a>"<wbr>,"r")</div><div>fmk = addfile("<a href="http://saji_landsea_mask.nc" target="_blank">saji_landsea_mask.nc</a>"<wbr>,"r")</div><div>rain = fin->var(:,{lt1:lt2},{ln1:ln2}<wbr>)</div><div>msk = fmk->ls_mask({lt1:lt2},{ln1:<wbr>ln2})</div><div>lmsk = conform(rain,msk,(/1,2/))</div><div><br></div><div>rain_masked = rain</div><div>rain_masked = (/where(lmsk.eq.1, rain@_FillValue,rain) /)</div><div>rjja_masked = dim_avg_n_Wrap(rain_masked,0)</div><div><br></div><div>wks = gsn_open_wks("x11","mask_fig")</div><div>gsn_define_colormap(wks,"<wbr>rainbow")</div><div>res = True</div><div>opt = True</div><div>res@mpMinLatF = lt1</div><div>res@mpMaxLatF = lt2</div><div>res@mpMinLonF = ln1</div><div>res@mpMaxLonF = ln2</div><div>res@gsnAddCyclic = False</div><div>res@cnFillOn = True</div><div>res@cnLinesOn = False</div><div>res@cnSmoothingOn = True ; <--- The problem is with this one</div><div>plt = gsn_csm_contour_map(wks,rjja_<wbr>masked,res)</div><div>; ----- the following one also has the strange issue that I reported ---</div><div>plt = gsn_csm_contour_map(wks,rain_<wbr>masked(3,:,:),res) ; only the data for the last time step</div></div><div><br></div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 26, 2018 at 12:14 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">Hi Saji,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">This is indeed strange.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Would it be possible for you to include a script and your data so we can look at the problem here? I especially want to see how the masking is being done, and I can test it using our 6.5.0 beta version.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can use our ftp for this, or put the files on your own site:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/ftp_files.shtml" target="_blank">http://www.ncl.ucar.edu/ftp_fi<wbr>les.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">Thanks,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-8266912830308425495h5">On Sun, Jun 24, 2018 at 8:02 PM, Saji Hameed <span dir="ltr"><<a href="mailto:saji.nh@gmail.com" target="_blank">saji.nh@gmail.com</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="m_-8266912830308425495h5"><div dir="ltr">Dear NCL developers,<div><br></div><div>I would like to report a rather strange issue that appears in version 6.4.0 of NCL.</div><div>The issue is related to shaded contours. For me, the issue happens when using a dataset for which a mask was applied to "undefine" values over a portion of the</div><div>domain. </div><div><br></div><div>I have attached four figures to illustrate the issue (two of these were made with</div><div>NCL version 6.3.0, and the other two with NCL 6.4.0).</div><div><br></div><div>nomask_ncl6.4.png and nomask_ncl6.3.png are plots of rainfall made using</div><div>NCL 6.4.0 and NCL 6.3.0 respectively. There are two panels, with the top</div><div>showing a shaded contour plot, and the bottom a simple contour plot.</div><div>In the cases shown above, all data in the domain were used. There are no issues with both the above mentioned plots (they are also identical to each other, as expected).</div><div><br></div><div>In the next set of plots, I have masked the data over land areas. Once I do this, I get</div><div>unexpected colors (white) for some contour levels (between 10 and 12) only when using NCL6.4.0.</div><div>There is no issue with NCL6.3.0, which plots the data without any apparent</div><div>problems.</div><div><br></div><div>I did note this issue with other plots, but at that time I thought that I was overusing colors on my desktop, and my OS could not allocate enough colors to NCL. But only today, I figured out that this is seemingly related to use of a mask to undefine</div><div>portions of the data domain.</div><div><br></div><div>For the time being, I will use NCL6.3.0, and hope that the source of this issue</div><div>can be found and fixed.</div><div><br></div><div>Best regards,</div><div>saji</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="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>