<div id="origbody"><div style="background: #f2f2f2;"><span style="background-color: rgb(255, 255, 255);">Dear Sir / Madam</span></div><p><span class="" style="white-space:pre">        </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&nbsp;Wet bulb temperature is available to be calculated in one step by wrf_user_getvar(file,"tw",time). However, when I used the way</p><p>to get it, I got a error massage:&nbsp;</p><p><br></p><p><span class="" style="white-space:pre">                </span>fatal:["Execute.c":6321]:variable (tw) is not in file (nc_file)</p><p><br></p><p>I think it means the function "wrf_user_getvar" is trying to get the "tw" variable from the wrfout file instead of calculating the diagnostics "we-</p><p>t bulb temperature".</p><p><br></p><p>Secondly, I tried to use the other way to calculate the "Tw" by using the function " wrf_wetbulb ", &nbsp;I also found some incorrect stuffs.On the&nbsp;</p><p>description page of the "wrf_wetbulb", it shows that we need to put argument "p" in mb (hPa). However in the given example 1, the script</p><p>uses "p &nbsp; = wrf_user_getvar(a,"p",0)" 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 "wrf_wetbulb", 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&nbsp;</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&nbsp;</p><p>you very very much!</p><p><br></p><p>Yours sincerely,</p><p>Shengxuan JI</p><div id="">--------------------------------<br></div><div id=""><br></div></div>