<p>Hi Mary,</p><p><br></p><p>Thank you for you reply. That is useful.</p><p>I have checked the description page of wrf_user_getvar, it has been updated. However, there is still a mistake about the unit of wet-bulb temp. On the page, it shows the unit is C, but I think it should be K. I have tried to get the wet-bulb temp by using "twb=wrf_user_getvar(a,"twb",0)" and print the summary of variable "twb". The result is:</p><p><img src="cid:part1.549b613b1a16a1"></p><p>please refer to it. Thank you very much!</p><p><br></p><p>Hopefully, my second question in my previous email could be answer. Someone may check and correct it, because it may mislead some new men in NCL.</p><p><br></p><p>Whatever, I love NCL, the wonderful software! : )</p><p><br></p><p>Merry Christmas!</p><p>Shengxuan</p><br><div id="">--------------------------------<br></div><div id=""><br></div><br><div id="origbody"><div style="background: #f2f2f2;">----- 原始邮件 -----<br>发件人：Mary Haley <haley@ucar.edu><br>收件人：shengxuan_ji@sina.com<br>抄送人：ncl-talk <ncl-talk@ucar.edu>, wrfhelp NCAR-MMM <wrfhelp@ucar.edu><br>主题：Re: [ncl-talk] I might find some bugs about calculation of wet-bulb temperature<br>日期：2014年12月24日 23点49分<br></div><br><div dir="ltr"><div class="" style="font-size:small">Hi,</div><div class="" style="font-size:small"><br></div><div class="" style="font-size:small">I do not know the answer to your second question about the correct units of "p" for wrf_wetbulb, but wrfhelp should be able to answer that.</div><div class="" style="font-size:small"><br></div><div class="" style="font-size:small">The documentation for "wrf_user_getvar" says to use "tw", but it should be "twb". I don't know if this is an error in the documentation or the WRFUserARW.ncl script, but I've updated the documentation to reflect this.</div><div class="" style="font-size:small"><br></div><div class="" style="font-size:small">--Mary</div><div class="" style="font-size:small"><br></div></div><div class=""><br><div class="">On Tue, Dec 23, 2014 at 12:26 AM, Mr.J <span dir="ltr"><<a target="_blank" href="mailto:shengxuan_ji@sina.com">shengxuan_ji@sina.com</a>></span> wrote:<br><blockquote class="" 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 target="_blank" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</div>