<div dir="ltr">My guess would be that some locations do not meet the criteria within the code for calculating CAPE -- these points are ignored. For instance, if the iteration fails within.<div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 10, 2021 at 3:11 AM Dennis Shea via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</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 dir="ltr"><div>NCL developers did not write CAPE_{2D/3D) fortran subroutines called by NCL. I have no experience with these functions.</div><div><br></div><div>I have attached the fortran code invoked by NCL. Perhaps the following subroutines will provide some insight.</div><div><br></div><div>DCAPECALC3D</div><div><div>DCAPECALC2D</div><div><br></div><div>Note that at the top of the code a disclaimer is issued:</div><div><br></div><div>"! Also, be advised that missing data values are not checked during the computation."</div><div>===</div><div>Perhaps some input variable has _FillValue</div><div><br></div><div>Good Luck</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 9, 2021 at 10:21 AM Wang, Jiali via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</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 lang="EN-US">
<div>
<p class="MsoNormal"><span style="font-size:11pt">Dear NCL users,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt">I have been using wrf_get_uservar to get cape_2d and cape_3d for my research domain, over the Great Lakes region. There are water bodies and land. I noticed there are many missing values in both cape_2d and
3d outcomes. See some examples attached. I have also tried to use wrf_cape_2d to calculate cape_2d, and I made sure the PSFC is converted to hPa. The outcome is almost the same, with lots of missing values. Missing values are located over different places
at different times. <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt">However, it seems the missing value is not expected?
<a href="https://www.ncl.ucar.edu/Support/talk_archives/2013/2422.html" target="_blank">https://www.ncl.ucar.edu/Support/talk_archives/2013/2422.html</a>
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt">Can you advice what is the correct way to calculate CAPE or whether I am getting reasonable results?<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt">Thank you!<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt">Jiali <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><img style="width: 5.6562in; height: 5.302in;" id="gmail-m_6346225659095524968gmail-m_2808232562085047435Picture_x0020_1" src="cid:17d075e87784cff311" width="543" height="509" border="0"><img style="width: 5.7291in; height: 5.3333in;" id="gmail-m_6346225659095524968gmail-m_2808232562085047435Picture_x0020_2" src="cid:17d075e87785b16b22" width="550" height="512" border="0"></span><span style="font-size:11pt"><u></u><u></u></span></p>
</div>
</div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><br></div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior <font color="#000000">Scientist, Lecturer,</font></div><div><div><font color="#000000">The Institute of Earth Sciences, <br>The Hebrew University of Jerusalem, <br>Givat Ram, Jerusalem 91904, Israel <br></font></div><font color="#000000">Tel: 972 547 231 170<br>Fax: (972)-25662581</font></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">https://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>