<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.gmail-apple-converted-space
{mso-style-name:gmail-apple-converted-space;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Dear Dave, Barry and Dennis,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you all for the detail notes, they are extremely helpful for me to continue exploring my problem. While I now understand the missing values are normal, I did test the difference in CAPE calculation between different versions as Dave
suggested. I do find differences between v6.5.0 and v6.6.2, and I also found an even bigger difference between v6.3.0 and v6.6.2. In v 6.3.0 (the only version I tested), there is a bug in CAPE calculation in WRFARWUser.ncl, the PSFC should be converted to
hPa. This is fixed in 6.5.0 and after (again these are the versions I tested, it seems v6.4.0 is also good from previous ncl-talk). See images below for JJA averaged CAPE_2D (maximum CAPE) in 2018, left v6.6.2 vs 6.3.0; right v6.6.2 vs 6.5.0.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><img width="545" height="540" style="width:5.677in;height:5.625in" id="Picture_x0020_1" src="cid:image001.png@01D7D6E9.B39CCA00"><img width="550" height="547" style="width:5.7291in;height:5.6979in" id="Picture_x0020_2" src="cid:image002.png@01D7D6E9.B39CCA00"><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">Dave Allured - NOAA Affiliate <dave.allured@noaa.gov><br>
<b>Date: </b>Tuesday, November 9, 2021 at 10:38 PM<br>
<b>To: </b>Barry Lynn <barry.h.lynn@gmail.com><br>
<b>Cc: </b>"Dennis edu>" <shea@ucar.edu>, Jiali Wang <jialiwang@anl.gov>, "ncl-talk@mailman.ucar.edu" <ncl-talk@mailman.ucar.edu><br>
<b>Subject: </b>Re: [ncl-talk] CAPE calculation<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal">Output missing values are <b>normal</b> for the WRF CAPE and CIN functions. Please refer to NCL documentation for the
<b>specific</b> functions <b>wrf_cape_2d</b> and <b>wrf_cape_3d</b>. The descriptions are confusing. They are trying to say in essence, "Missing values are returned for out-of-range values". It is also not clear what "out of range" means. With a quick look
at the fortran code, it seems that "out of range" means very low energy values by a combination of two different tests.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I recently used <b>wrf_cape_2d</b> with global grids from a climate model. Spot checking gave these results which are similar to Jiali's graphical results.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">CAPE: 64% to 71% missing <b>globally</b>. Minimum value =
<b>0.10000 J/kg</b>.<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">CIN: 77% to 80% missing<span class="gmail-apple-converted-space"> </span><b>globally</b>. Minimum value =<span class="gmail-apple-converted-space"> </span><b>0.10000 J/kg</b>.<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The large difference in proportion missing between Jiali's case and mine may be explained because Jiali's region is central US, a land region with a reputation for high convective activity. Presumably a distribution of global grid points
including oceans and arctic regions averages to much less convective energy per grid point, therefore more out-of-range low values. Also I suspect that the proportion missing has a seasonal variation, but I did not do enough checks to confirm this.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">On a final note, there was an undocumented change in the code for the CAPE and CIN functions, somewhere between NCL versions 6.5.0 and 6.6.2. The two versions produce significantly different numeric results! The fortran 90 CAPE code in
6.6.2 is obviously more modern than the corresponding fortran 77 code in 6.5.0. I recommend that only NCL version 6.6.2 or later should be used to calculate CAPE and CIN.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Tue, Nov 9, 2021 at 8:39 PM Barry Lynn via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">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.<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Barry<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Wed, Nov 10, 2021 at 3:11 AM Dennis Shea via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal">NCL developers did not write CAPE_{2D/3D) fortran subroutines called by NCL. I have no experience with these functions.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I have attached the fortran code invoked by NCL. Perhaps the following subroutines will provide some insight.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">DCAPECALC3D<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">DCAPECALC2D<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Note that at the top of the code a disclaimer is issued:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">"! Also, be advised that missing data values are not checked during the computation."<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">===<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Perhaps some input variable has _FillValue<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Good Luck<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">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:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Dear NCL users,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">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. <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">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> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Can you advice what is the correct way to calculate CAPE or whether I am getting reasonable results?<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thank you!<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Jiali<o:p></o:p></p>
</div>
</div>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</div>
</div>
</div>
</body>
</html>