<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hello,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think you need to check your input units again.  You didn't indicate exactly how you are calling wrf_cape_3d, but if you read the documentation at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_cape_3d.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_cape_3d.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">you will see that that the model pressure values need to be in units of Pa, while the units of  surface pressure must be hPa. I'm not sure why the function was written this way, but you need to make sure you have the correct units.</div><div class="gmail_default"><br></div><div class="gmail_default">If you continue to have problems, please respond back to ncl-talk with this information. It would help if you could provide a short test program that shows the values you are inputting and the incorrect values that are coming from the function. Also, we need to know what version of NCL you're running. </div><div class="gmail_default"><br></div><div class="gmail_default">Thanks,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Jan 8, 2019 at 9:23 PM nafiseh pegahfar <<a href="mailto:pegahfar@inio.ac.ir">pegahfar@inio.ac.ir</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>
<div style="font-size:13.3333px;font-family:"Times New Roman";color:rgb(0,0,0);font-weight:400;font-style:normal;background:none 0% 0%/auto repeat scroll padding-box border-box rgba(0,0,0,0)">
        <pre style="white-space:pre-wrap">To whom it may concern,

I used "wrf_cape_3d" function to calculate Cape. for this aim I 
used 3D vand 
2D values. unfortunately in some points (based on lat/lon)  NaN value for 
cape is seen. I am sure that all data have correct units. Input values for a 

point is reported as below:

cape=9.969209968386869e+36

level=92500
TK=300.6299140816899
rv=0.008510922866573187
GPH=723.139880629395
z=0.7325572010323265
p=100078.4522848937
I think where rv value is less that 0.01, this nan value appears.
please help me to overcome to these nan values

</pre>
</div>
<div style="font-family:"Times New Roman";color:rgb(0,0,0);font-size:13.3333px;font-weight:400;font-style:normal;background:none 0% 0%/auto repeat scroll padding-box border-box rgba(0,0,0,0)">
        Best Regards<br>
        =================================<br>
        Nafiseh Pegahfar<br>
        Assistant Professor<br>
        Iranian National Institute for Oceanography and Atmospheric Science<br>
        (<a href="http://www.inio.ac.ir" target="_blank">http://www.inio.ac.ir</a>)<br>
        Phone: (0098)21- 66944873-5 Ext. 224<br>
        Fax: (0098)21- 66944869<br>
        Email: (<a href="mailto:pegahfar@inio.ac.ir" target="_blank">pegahfar@inio.ac.ir</a>)<br>
        (<a href="mailto:pegahfar@ut.ac.ir" target="_blank">pegahfar@ut.ac.ir</a>)<br>
        =================================</div>
</div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>