<div dir="ltr">Hi Kai,<div><br></div><div>In your printVarSummary, your 'ps' variable is in units of Pa and this needs to be hPa, so please multiply by .01 before passing it to the wrf_cape_2d function and see if that helps. (However, keep your 'p' variable in units of Pa). Also, many of your variables have a _FillValue attribute which indicates that there might be missing/fill values in your arrays. You can check to see if you have any missing values in these arrays by doing:</div><div><br></div><div><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal">if(<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/any.shtml" style="margin:0px;padding:0px;text-decoration:none;color:rgb(133,45,133);font-family:verdana,sans-serif"><strong style="margin:0px;padding:0px">any</strong></a>(<strong style="margin:0px;padding:0px">ismissing</strong>(your_array))) then</pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"> print ("found missing value")</pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal">end if</pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><br></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><font face="arial, helvetica, sans-serif">If none of that helps, and your data file isn't gigantic, you can FTP following these instructions </font><span style="font-family:arial,helvetica,sans-serif;font-size:13.3333px">(taken from <a href="http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP">http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP</a>)</span><span style="font-family:arial,helvetica,sans-serif;font-size:13.3333px">:</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><p style="margin:10px 0px 20px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px;font-variant-ligatures:normal;white-space:normal"><b style="margin:0px;padding:0px">You can ftp larger datasets and files to:</b></p><p style="margin:10px 0px 20px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px;font-variant-ligatures:normal;white-space:normal"></p><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;font-variant-ligatures:normal">ftp <a href="http://ftp.cgd.ucar.edu">ftp.cgd.ucar.edu</a>
anonymous
<use your email address for the password>
cd incoming
put ...
put ...
.
.
.
quit
</pre><p style="margin:10px 0px 20px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px;font-variant-ligatures:normal;white-space:normal"><font color="red" style="margin:0px;padding:0px">Note that for security reasons, you cannot list the contents of this directory, and neither can we. We need to know the <i style="margin:0px;padding:0px">exact</i> names of the files in order to retrieve them.</font></p></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><font face="arial, helvetica, sans-serif">After you've uploaded the file, send me an email directly to tell me the name of the file. If for some reason the file is too big to upload, </font><span style="font-family:arial,helvetica,sans-serif;font-size:13.3333px">you can upload to dropbox, or some other file sharing service and email me the link.</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><span style="font-family:arial,helvetica,sans-serif;font-size:13.3333px">Regards,</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><span style="font-family:arial,helvetica,sans-serif;font-size:13.3333px">Bill</span></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><br></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><br></pre><pre style="margin-top:0px;margin-bottom:0px;padding:0px;font-size:13.3333px;line-height:16px;font-family:courier;border-style:solid;border-width:0px 0px 15px;border-color:transparent;color:rgb(0,0,0);font-variant-ligatures:normal"><br></pre></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 12, 2016 at 8:19 PM, Kai Zhang <span dir="ltr"><<a href="mailto:kzkaizhang@gmail.com" target="_blank">kzkaizhang@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">Dear Bill,</div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large"><br>Thanks so much for your response!</div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large"><br></div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">I've check in several ways that, all of the units are correct. And non missing value, non strange value. The problem still exists.</div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">Here are the result for printVarSummary on each of my input variables:</div><div class="gmail_default">
<p><span><font face="arial, helvetica, sans-serif">Variable: p</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 595200 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 148800 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 4</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> hour: [0..21]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lev: [10000..100000]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lat: [-48.31578947368421..8.<wbr>526315789473685]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lon: [ 270.. 330]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>Pa</font></span></p>
<p><font face="arial, helvetica, sans-serif"><span></span><br></font></p>
<p><span><font face="arial, helvetica, sans-serif">Variable: T</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 595200 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 148800 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 4</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> hour: [0..21]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lev: [10000..100000]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lat: [-48.31578947368421..8.<wbr>526315789473685]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lon: [ 270.. 330]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 5</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> day :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> year :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>K</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> long_name :<span> </span>Temperature</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> _FillValue :<span> </span>9.96921e+36</font></span></p>
<p><font face="arial, helvetica, sans-serif"><span></span><br></font></p>
<p><span><font face="arial, helvetica, sans-serif">Variable: Q</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 595200 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 148800 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 4</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> hour: [0..21]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lev: [10000..100000]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lat: [-48.31578947368421..8.<wbr>526315789473685]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lon: [ 270.. 330]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 5</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> day :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> year :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>kg/kg</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> long_name :<span> </span>Specific humidity</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> _FillValue :<span> </span>9.96921e+36</font></span></p>
<p><font face="arial, helvetica, sans-serif"><span></span><br></font></p>
<p><span><font face="arial, helvetica, sans-serif">Variable: Z</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 595200 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 148800 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 4</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[hour | 8] x [lev | 24] x [lat | 31] x [lon | 25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> hour: [0..21]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lev: [10000..100000]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lat: [-48.31578947368421..8.<wbr>526315789473685]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> lon: [ 270.. 330]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 5</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> day :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> year :<span> </span>1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>m</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> long_name :<span> </span>Geopotential Height (above sea level)</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> _FillValue :<span> </span>9.96921e+36</font></span></p>
<p><font face="arial, helvetica, sans-serif"><span></span><br></font></p>
<p><span><font face="arial, helvetica, sans-serif">Variable: zs</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 24800 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 6200 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 3</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[8] x [31] x [25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 2</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>m</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> _FillValue :<span> </span>1e+15</font></span></p>
<p><font face="arial, helvetica, sans-serif"><span></span><br></font></p>
<p><span><font face="arial, helvetica, sans-serif">Variable: ps</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Type: float</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Total Size: 24800 bytes</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 6200 values</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number of Dimensions: 3</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Dimensions and sizes:<span> </span>[8] x [31] x [25]</font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Coordinates: </font></span></p>
<p><span><font face="arial, helvetica, sans-serif">Number Of Attributes: 1</font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> units :<span> </span>Pa</font></span></p><span class="">
<p><span><font face="arial, helvetica, sans-serif"> capecalc3d: Outside of lookup table bounds. prs,thte= 925.000000000000 </font></span></p>
<p><span><font face="arial, helvetica, sans-serif"> 856.034700407579 </font></span></p><p><span><font face="arial, helvetica, sans-serif"><br></font></span></p></span><p><span><font face="times new roman, serif" size="4">Thanks again!</font></span></p><p><span><font face="times new roman, serif" size="4">Best,</font></span></p><p><span><font face="times new roman, serif" size="4">Kai</font></span></p></div><img src="http://t.sidekickopen65.com/e1t/o/5/f18dQhb0S7ks8dDMPbW2n0x6l2B9gXrN7sKj6v4fhhMW1p7Y9C7dSGy-W4X9YsK1pctGFW8NqsSQ1k1H6H0?si=4754807383916544&pi=02e16246-a903-4cde-eb8d-b3de50929a20" style="display:none!important" height="1" width="1"></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Mon, Sep 12, 2016 at 10:41 AM, Bill Ladwig <span dir="ltr"><<a href="mailto:ladwig@ucar.edu" target="_blank">ladwig@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Kai,<div><br></div><div>wrf_cape_2d requires the pressure argument (first argument) to be in units of Pa and the surface pressure argument (6th argument) to be in units of hPa, so double check that. Temperature needs to be in Kelvin. Height is in meters. Water vapor mixing ratio is in units of kg/kg. The error message you are getting is indicating that theta-e is being calculated to 856.03 K, which is obviously nonsensical, so one of your input arrays is in the wrong units, contains missing values, or has a bogus number in it. This algorithm relies on lookup tables, which is why things can go wrong in a hurry if something is not right. If you're still having problems, please reply with the result for printVarSummary on each of your input variables.</div><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Fri, Sep 9, 2016 at 8:34 PM, Kai Zhang <span dir="ltr"><<a href="mailto:kzkaizhang@gmail.com" target="_blank">kzkaizhang@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">Dear NCL support,</div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large"><br></div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">I am currently using function wrf_cape_2d to calculate CIN and CAPE. But I cannot stop getting errors:</div><div class="gmail_default" style="font-family:"times new roman",serif;font-size:large">
<p><span>"capecalc3d: Outside of lookup table bounds. prs,thte= 925.000000000000 </span>856.034700407579"</p><p>I am sure that I got units for each variable correctly. </p><p>It is weird that, if I change "p (unit: Pa)" to "p (unit: hPa)", I can get it run smoothly. But the CIN and CAPE calculated are obviously wrong.</p><p>Do you know why?</p><p><br></p><p>Thanks a lot for your help!</p><p>Best,</p><p>Kai</p></div>
<img src="http://t.sidekickopen65.com/e1t/o/5/f18dQhb0S7ks8dDMPbW2n0x6l2B9gXrN7sKj6v4fhhMW1p7Y9C7dSGy-W4X9YsK1pctGFW8NqsSQ1k1H6H0?si=4754807383916544&pi=5a3a3138-7c1f-4942-8f23-e57478f7b21f" style="display:none!important" height="1" width="1"></div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><div>------------------------------<wbr>------------------------------<wbr>-------------------</div><div>Kai Zhang</div><div>The University of Texas at Austin, TX (2012 - )</div><div>Graduate Research Assistant,</div><div>Jackson School of Geosciences</div><div><div>BS, Peking University, China (2008 - 2012)</div><div>Dept. Atmospheric and Oceanic Science, School of Physics</div></div><div><div>Cellphone: <a href="tel:%2B1%20512%20520%206902" value="+15125206902" target="_blank">+1 512 520 6902</a></div><div><a href="https://sites.google.com/site/kaizhanghomepage/" target="_blank">https://sites.google.com/site/<wbr>kaizhanghomepage/</a><br></div><div>Email: <a href="mailto:zhangkai@utexas.edu" style="color:rgb(17,85,204)" target="_blank">kzkaizhang@utexas.edu</a></div></div><div> <a href="mailto:kzkaizhang@gmail.com" target="_blank">kzkaizhang@gmail.com</a></div><div><br></div></span></div></div></div></div>
</font></span></div>
</blockquote></div><br></div>