<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I do not know the answer to your second question about the correct units of &quot;p&quot; for wrf_wetbulb, but wrfhelp should be able to answer that.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The documentation for &quot;wrf_user_getvar&quot; says to use &quot;tw&quot;, but it should be &quot;twb&quot;.  I don&#39;t know if this is an error in the documentation or the WRFUserARW.ncl script, but I&#39;ve updated the documentation to reflect this.</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 Tue, Dec 23, 2014 at 12:26 AM, Mr.J <span dir="ltr">&lt;<a href="mailto:shengxuan_ji@sina.com" target="_blank">shengxuan_ji@sina.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="background:#f2f2f2"><span style="background-color:rgb(255,255,255)">Dear Sir / Madam</span></div><p><span style="white-space:pre-wrap">        </span></p><p>Thank you for your hard working on the development of NCL software. I can easily manipulate the met data and professional plots now.</p><p><br></p><p>During my usage of wrf related functions, I found a bug, which is about calculation of wet-bulb temperature. I noticed that after the version</p><p>6.2.0, the Wet bulb temperature is available to be calculated in one step by wrf_user_getvar(file,&quot;tw&quot;,time). However, when I used the way</p><p>to get it, I got a error massage: </p><p><br></p><p><span style="white-space:pre-wrap">                </span>fatal:[&quot;Execute.c&quot;:6321]:variable (tw) is not in file (nc_file)</p><p><br></p><p>I think it means the function &quot;wrf_user_getvar&quot; is trying to get the &quot;tw&quot; variable from the wrfout file instead of calculating the diagnostics &quot;we-</p><p>t bulb temperature&quot;.</p><p><br></p><p>Secondly, I tried to use the other way to calculate the &quot;Tw&quot; by using the function &quot; wrf_wetbulb &quot;,  I also found some incorrect stuffs.On the </p><p>description page of the &quot;wrf_wetbulb&quot;, it shows that we need to put argument &quot;p&quot; in mb (hPa). However in the given example 1, the script</p><p>uses &quot;p   = wrf_user_getvar(a,&quot;p&quot;,0)&quot; to get the pressure. It may be wrong, because the unit is Pa calculating in this way. All in all, what unit</p><p>should I use for the pressure.<font color="#cc0000"> (Please refer to test_scrpit1.ncl)</font></p><p><br></p><p>Thirdly, in the description page of function &quot;wrf_wetbulb&quot;, the unit of return value is [C]. When I get the results, these values are larger than</p><p>250. So, I am confused about the return value unit. Is it [C] or [K]?<font color="#cc0000"> (Please refer to test_script2.ncl)</font></p><p><br></p><p>The attachments are my test scripts, please check them. The wrfout file is too large to put it on a website, you may find</p><p>another one to try the scripts.</p><p><br></p><p>My system details:</p><p><b><font color="#cc0000">ncl version: </font></b><b style="color:rgb(102,194,255)">6.2.0</b></p><p><font color="#66c2ff"><b>( I also tried the vesion 6.2.1, the results are the same.)</b></font></p><p><font color="#cc0000"><b><br></b></font></p><p><font color="#cc0000"><b>uname -a </b></font></p><p><font color="#66c2ff"><b>Linux nowcast 2.6.32-358.el6.x86_64 #1 SMP Tue Jan 29 11:47:41 EST 2013 x86_64 x86_64 x86_64 GNU/Linux</b></font></p><p><b><font color="#cc0000"><br></font></b></p><p><b><font color="#cc0000">gcc --version</font></b></p><p><font color="#66c2ff"><b>gcc (GCC) 4.4.7 20120313 (Red Hat 4.4.7-3)</b></font></p><p><font color="#66c2ff"><b></b></font></p><p><font color="#66c2ff"><b>Copyright (C) 2010 Free Software Foundation, Inc.</b></font></p><p><b><font color="#cc0000"><br></font></b></p><p><b><font color="#cc0000">cat /etc/redhat-release</font></b></p><p><font color="#66c2ff"><b>Red Hat Enterprise Linux Server release 6.4 (Santiago)</b></font></p><p><br></p><p><br></p><p>To be honest, I am not sure whether these are bugs or not. Please help me to deal with these questions. Hope your reply ASAP. Thank </p><p>you very very much!</p><p><br></p><p>Yours sincerely,</p><p>Shengxuan JI</p><div>--------------------------------<br></div><div><br></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>