<div dir="ltr"><div><div><div>Hello, Your question/post has been forwarded to <a href="mailto:wrfhelp@ucar.edu">wrfhelp@ucar.edu</a><br><br></div>This is the holiday season so many people are taking vacation.<br><br></div>Regards<br></div>ncl-talk<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 23, 2014 at 12:26 AM, Mr.J <span dir="ltr"><<a href="mailto:shengxuan_ji@sina.com" target="_blank">shengxuan_ji@sina.com</a>></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,"tw",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:["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 ", I also found some incorrect stuffs.On the </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 = 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 </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>