<div dir="ltr"><div>[0] wrfhelp would likely be the best place to ask about this.<br>---<br>[1] Attached are some fortran subroutines to calculate various quantities using WRF data. <br>   wrf_curl_div.f: <a href="http://forum.wrfforum.com/viewtopic.php?f=32&t=2267">http://forum.wrfforum.com/viewtopic.php?f=32&t=2267</a><br>   wrf_vort.f:       <a href="http://mailman.ucar.edu/pipermail/ncl-talk/2014-November/001483.html">http://mailman.ucar.edu/pipermail/ncl-talk/2014-November/001483.html</a><br><br></div><div>These can be invoked from an NCL script by using WRAPIT:<br>   <a href="https://www.ncl.ucar.edu/Document/Tools/WRAPIT.shtml">https://www.ncl.ucar.edu/Document/Tools/WRAPIT.shtml</a><br></div><div>Likely, you can ignore warning messages.<br></div><div><br>%> WRAPIT wrf_curl_div.f<br>%> WRAPIT wrf_vort.f<br><br></div><div><br></div><div>NOTE: these were *not written or tested* by NCL developers. Hence, they are *not supported*.<br><br></div><div>==================<br></div><div>[2] An alternative that may be adequate if the divergence is being calculated for visual purposes is:<br><br>   <a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml">https://www.ncl.ucar.edu/Applications/ESMF.shtml</a><br><br>   (a) interpolate the zonal and meridional wind components to a rectilinear grid; <br>  
(b) Use the <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/uv2dv_cfd.shtml"><strong>uv2dv_cfd</strong></a> function to compute divergence on the
rectilinear grid; <br>   (c) interpolate 
the derived divergence (b) onto the original WRF grid.
<br><br></div><div>Good Luck<br><br></div><div><br><br><br></div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 22, 2018 at 1:44 AM, 古枫玥瑟 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>Dear ncl-talk,</div><div>There are 3 functions to c<span style="font-family:verdana,sans-serif;font-size:13.3333px">ompute the divergence. They are </span><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/uv2dvF-1.shtml" style="margin:0px;padding:0px;text-decoration:none;color:rgb(133,45,133);font-family:verdana,sans-serif;font-size:13.3333px" target="_blank">uv2dvF</a>, <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/uv2dvF-1.shtml" style="margin:0px;padding:0px;text-decoration:none;color:rgb(133,45,133);font-family:verdana,sans-serif;font-size:13.3333px" target="_blank">uv2dv</a>G,<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/uv2dv_cfd.shtml" style="margin:0px;padding:0px;text-decoration:none;color:rgb(133,45,133);font-family:verdana,sans-serif;background:rgb(205,232,251);font-size:13.3333px" target="_blank">uv2dv_cfd</a>.</div><div>I would like to know which function should I use to c<span style="font-family:verdana,sans-serif;font-size:13.3333px">ompute the divergence given WRF output file (eg. </span><font face="verdana, sans-serif"><span style="font-size:13.3333px">wrfout_d03_2016-08-11_00:<wbr>00:00</span></font><span style="font-family:verdana,sans-serif;font-size:13.3333px">) because the number of dimensions of the lat/lon is 2 ?</span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px">Then how to write the code? Thank you!</span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px">                              </span><span style="white-space:pre-wrap;font-family:"lucida Grande",Verdana,"Microsoft YaHei";background-color:rgb(255,231,198)">Many thanks in advance</span></div><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word;font-family:"lucida Grande",Verdana,"Microsoft YaHei""><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word"><font style="line-height:47.6px" face="verdana, sans-serif"><span style="font-size:13.3333px;line-height:16px;background-color:rgb(255,231,198)"><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;font-size:14px;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word;font-family:"lucida Grande",Verdana,"Microsoft YaHei""><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word"><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer" style="background-color:rgb(255,255,255);line-height:23.8px"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word"><div class="m_-6994153169287473787mailcontainer" id="m_-6994153169287473787qqmail_mailcontainer" style="line-height:23.8px"><div id="m_-6994153169287473787mainmail" style="margin-bottom:12px"><div id="m_-6994153169287473787contentDiv" class="m_-6994153169287473787body" style="line-height:1.7;height:auto;padding:15px 15px 10px;zoom:1"><div id="m_-6994153169287473787qm_con_body"><div id="m_-6994153169287473787mailContentContainer" class="m_-6994153169287473787qmbox m_-6994153169287473787qm_con_body_content m_-6994153169287473787qqmail_webmail_only" style="padding:0px;height:auto;min-height:100px;word-wrap:break-word"><blockquote type="cite" style="line-height:23.8px"><div class="m_-6994153169287473787WordSection1"><pre style="line-height:23.8px"><span lang="EN-US">Regards<u></u><u></u></span></pre><pre style="line-height:23.8px"><span lang="EN-US">Yan<u></u><u></u></span></pre><pre style="line-height:23.8px"><span lang="EN-US">Shanghai Institute of Meteorology,China.</span></pre></div></blockquote></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></span></font></div></div></div></div></div></div></div></div></div></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px"><br></span></div><div><span style="font-family:verdana,sans-serif;font-size:13.3333px">             </span></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>