<div dir="ltr"><div class="gmail_default" style="font-size:small">Barry,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It is hard to debug a problem like this, especially given the fact that you are calling a Fortran code to get the ij indices, and we have no idea if this code is returning the correct values. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It looks to me like your mpRight/LeftCorner settings are probably incorrect, and that the values you are using for the lower right corner is outside the range of your and this could be due to GET_IJ being incorrect. Is it possible that GET_IJ is returning indices from 1 to N, and not 0 to N-1? If so, then you may need to subtract 1 from the return values:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace"> x_start = i_loc-1</font></div><div class="gmail_default"><font face="monospace, monospace"> y_start = j_loc-1</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace"> . . .</font></div><div class="gmail_default"><font face="monospace, monospace"> x_end = i_loc-1</font></div><div class="gmail_default"><font face="monospace, monospace"> y_end = j_loc-1</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div>If this still doesn't look correct, then I recommend plotting your data without zooming in on it, and add the locations of (lat_dom(0,0),lon(0,0)) and (lat_dom(njj,nii),lon_dom(njj,nii)) as markers on your plot (use gsn_add_polymarker) to see if they fall within the range of your lat/lon domain.</div><div class="gmail_default"><br></div><div class="gmail_default">I've included a sample script that you can use as a start. You should be able to just change the name of the WRF output file and the x_start/x_end/y_start/y_end values to match what you have in your own script, and see if you get the expected results.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jan 14, 2017 at 12:53 PM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@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:<div><br></div><div>I have attached both the output and ncl program that reads WRF data and plots using gsn.</div><div><br></div><div>It sets the grid as "Lambert Conformal" and reads the true latitudes and standard longitude from the WRF file.</div><div><br></div><div>It also overlays data using an Opacity factor.</div><div><br></div><div>Lastly, it uses "mask" to plot land vs ocean topography over which the snow forecast is overlaid.</div><div><br></div><div>Thank you for your help.</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 12, 2017 at 7:15 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@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">Hi Barry,<div>I am guessing that the data is being drawn to the boundary of the data coordinates of your overlay plot. You will not be able to get NCL to draw anything outside of the data coordinates. Thus, short of interpolating your data to a different grid, I believe your only choice would be to try different map projections and map limits. This will only truly be done via trial-and-error. You will want to set mpProjection along with various mpLimitMode related resources. (You are already setting mpMinLatF/mpMaxLatF/mpMinLonF<wbr>/mpMaxLonF with mpLimitMode = "LatLon" (the default)., and plotting on a cylindrical equidistant projection (again, the default).)</div><div><br></div><div>The NCL examples website has numerous examples. <a href="http://www.ncl.ucar.edu/Applications/" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/</a>, specifically see the Map Projections examples. As you stated that your data is on a Lambert grid, you might want to start with the Lambert Conformal examples. </div><div><br></div><div>Hope that helps. If not, please respond to the ncl-talk email list.</div><div>Adam</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_6939181481719391749h5">On Mon, Jan 9, 2017 at 9:17 AM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@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_6939181481719391749h5"><div dir="ltr">Hi:<div><br></div><div>I am excerpting a part of the WRF domain to make a snow accumulation plot.</div><div><br></div><div> snowh_rate_dom = (/snowh_rate(y_start:y_end,x_s<wbr>tart:x_end)/)</div><div><br clear="all"><div>Perhaps because the WRF domain (WPS) setting was Lambert, I get a blank area on the far right lower side (the overlay of snow does not extend all the way to the boundary).</div><div><br></div><div>is there a way to force/cut-out the overlay so that it falls to the boundary?</div><div><br></div><div>Barry</div>-- <br><div class="m_6939181481719391749m_8624959363627858631m_656014754122191997gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div></div>
<br></div></div><span>______________________________<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></span></blockquote></div><span class="m_6939181481719391749HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_6939181481719391749m_8624959363627858631gmail_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/asphill<wbr>i/</a> </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </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>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_6939181481719391749gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</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></div>