<div dir="ltr">Scott:<div><br></div><div>I have attached two programs.</div><div><br></div><div>The second simply maps the cross-section on a horizontal map, so you can see what line/cross-section you've chosen.</div><div><br></div><div>I separated this program from the one that makes vertical cross-sections because I didn't want to mix horizontal and vertical cross-section maps in the same output.</div><div><br></div><div>The first creates the cross section within specified lat/lon and vertical heights using wrf-zoom, for instance.</div><div><br></div><div>It is set up for non-standard radar output, but you can make any cross-sections of WRF output that you want. It only makes one direction at a time (you would have to choose a different IP values in the loop).</div><div><br></div><div>Barry </div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jan 21, 2018 at 8:42 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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5"><div dir="ltr"><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Scott Capps</b> <span dir="ltr"><<a href="mailto:scapps@atmosdatasolutions.com" target="_blank">scapps@atmosdatasolutions.com</a><wbr>></span><br>Date: Fri, Jan 19, 2018 at 5:12 PM<br>Subject: Re: [ncl-talk] WRF Cross-section Plot - Terrain Depiction and Lowest Level Data Contours<br>To: Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>><br>Cc: Bill Ladwig <<a href="mailto:ladwig@ucar.edu" target="_blank">ladwig@ucar.edu</a>>, NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br><br><br><div dir="ltr">I modify the placement and width by passing in the lat/lon of the two points. I would be interested in viewing your code and can send mine as well.</div><div class="m_-6092361436124769203HOEnZb"><div class="m_-6092361436124769203h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 19, 2018 at 4:22 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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Scott:<div><br></div><div>I am curious if you are able to set the cross section to appear where you want it, or if you are able to modify the width (cross-sectional extent)?</div><div><br></div><div>If not, I can send you my cross-section code that does all that.</div><div><br></div><div>Barry Lynn</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 19, 2018 at 1:55 AM, Bill Ladwig via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@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 Scott,<div><br></div><div>My first guess is that the vertical levels for interpolation are too large (not your fault). The interpolation algorithm automatically chooses vertical levels at .01*z_max. In your case z_max is 20251.8 m, which means your interpolation levels are at 202.5 m per grid box. This is probably too large for what you're interested in looking at. <br></div><div><br></div><div>As a first attempt to fix, let's edit your WRFUserARW.ncl.</div><div><br></div><div>At line 383, change this:</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<p class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851m_-5331414785227793666gmail-p1" style="margin:0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Monaco"><span class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851m_-5331414785227793666gmail-s1" style="text-decoration:underline">dz</span> = 0.01 * z_max</p></blockquote><div><br></div><div>to this:</div><div><br></div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851m_-5331414785227793666gmail-s1" style="font-family:Monaco;font-size:11px;text-decoration:underline">dz</span><span style="font-family:Monaco;font-size:11px"> = 0.0025 * z_max</span></blockquote><div> </div><div>If this doesn't help your situation, I'll need to dig in to this a bit more.</div><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851h5">On Thu, Jan 18, 2018 at 3:04 PM, Scott Capps via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.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="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851h5"><div dir="ltr">Greetings,<div><br></div><div>I have some questions about the cross-section plots from WRF output. I basically started from the example scripts at: <a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/CROSS_SECTION/wrf_CrossSection_smooth4.htm" target="_blank">http://www2.mmm.ucar.edu/<wbr>wrf/OnLineTutorial/Graphics/NC<wbr>L/Examples/CROSS_SECTION/wrf_C<wbr>rossSection_smooth4.htm</a></div><div><br></div><div>The whitespace along the terrain height is much too large and the terrain is too coarse (this is a 667m resolution WRF run) or the lowest level atmospheric data is not being contoured correctly. The cross-section plotting routine in GRaDs seems to show true the resolution of terrain and terrain following data. Is there a way to improve this in NCL?</div><div><br></div><div>I have attached some code (cross_sxn_plot_test.ncl and modified WRFUserARW.ncl) and have ftp'd an example wrfout file to the incoming directory as: wrfout_d05_2017-10-09_00:00:<a href="http://00.nc" target="_blank">00<wbr>.nc</a></div><div><br></div><div><br></div><div>NCL v6.3.0</div><div><a href="http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/CROSS_SECTION/wrf_CrossSection_smooth4.htm" target="_blank">http://www2.mmm.ucar.edu/wrf/O<wbr>nLineTutorial/Graphics/NCL/Exa<wbr>mples/CROSS_SECTION/wrf_CrossS<wbr>ection_smooth4.htm</a><span class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851m_-5331414785227793666HOEnZb"><font color="#888888"><br clear="all"><div><br></div>-- <br><div class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851m_-5331414785227793666m_8360096949610619015gmail_signature"><div dir="ltr"><div><div dir="ltr"><b><font size="4">Scott Capps</font></b><div>Principal</div><div><a href="http://www.atmosphericdatasolutions.com" target="_blank">Atmospheric Data Solutions</a></div><div><a href="mailto:scapps@atmosdatasolutions.com" target="_blank">scapps@atmosdatasolutions.com</a><div>cell: <a href="tel:(949)%20910-4385" value="+19499104385" target="_blank">(949) 910-4385</a></div><div><br></div></div></div></div></div></div>
</font></span></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>
<br>______________________________<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><br clear="all"><div><br></div>-- <br><div class="m_-6092361436124769203m_-9203516629450714248m_6599770701295368851gmail_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>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_-6092361436124769203m_-9203516629450714248gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><b><font size="4">Scott Capps</font></b><div>Principal</div><div><a href="http://www.atmosphericdatasolutions.com" target="_blank">Atmospheric Data Solutions</a></div><div><a href="mailto:scapps@atmosdatasolutions.com" target="_blank">scapps@atmosdatasolutions.com</a><div>cell: <a href="tel:(949)%20910-4385" value="+19499104385" target="_blank">(949) 910-4385</a></div><div><br></div></div></div></div></div></div>
</div>
</div></div></div><br><br clear="all"><div><br></div>-- <br><div class="m_-6092361436124769203gmail_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></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_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 914 432 3108<br></div></div></div>
</div>