<div dir="ltr"><div class="gmail_default" style="font-family:garamond,times new roman,serif">Dave, <br></div><div class="gmail_default" style="font-family:garamond,times new roman,serif">Thank you for pointing to this discussion. It's informative for sure - but doesn't answer my basic questions. If I am unable to derive this information from the code, I'll post it on this forum. <br></div><div class="gmail_default" style="font-family:garamond,times new roman,serif"><br></div><div class="gmail_default" style="font-family:garamond,times new roman,serif">Thank you, <br></div><div class="gmail_default" style="font-family:garamond,times new roman,serif">P. <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 5, 2022 at 12:38 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov">dave.allured@noaa.gov</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>Prashanth, please remember to include the mailing list in your replies.</div><div><br></div><div>It turns out that there are two different WRF forums.  The other one seems to be more active.  In particular, there is a modern discussion of CAPE behavior and parameterization in the post processing/RIP category.</div><div><br></div><div>    <a href="https://forum.mmm.ucar.edu/" target="_blank">https://forum.mmm.ucar.edu/</a></div><div><br></div><div>Also note, NCL upgraded the CAPE routine in NCL 6.6.2.  I think that is the version that Dennis sent.  The older version in 6.5.0 was in fortran 77, and the results were somewhat different.  I would recommend the newer version over the older one.</div><div><br></div><div>Last time I checked, late last year, the CAPE routine in WRF was identical to the NCL 6.6.2 version, except for a few "cosmetic" differences.</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 5, 2022 at 7:26 AM Prashanth Bhalachandran <<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank">prashanth.bhalachandran@gmail.com</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="auto"><div dir="ltr"></div><div dir="ltr">Thanks, Dave. I’ve gone through the various WRF forums but I am unable to find what I specifically need. I am a bit weak in Fortran, so I am struggling to find answers to the four basic elements (the parcel, the environment, the levels of free convection, and neutral buoyancy) needed to compute CAPE but I guess with enough persistence, I should extract it from the Fortran code that Dennis sent. </div><div dir="ltr"><br></div><div dir="ltr">Prashanth </div><div dir="ltr"><br><blockquote type="cite">On May 4, 2022, at 2:58 PM, Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>> wrote:<br></blockquote></div><blockquote type="cite"><div dir="ltr">Prashanth, the CAPE routines were copied from the WRF project.  There are active discussions about their software in the WRF forums.  I suggest look in the Graphics and Post-Processing category.<div dir="ltr"><div dir="ltr"><br></div><div dir="ltr">    <a href="https://wrfforum.com" target="_blank">https://wrfforum.com</a><br></div><div dir="ltr"><br></div><div>You can probably use the NCL wrap functions to call your own customized CAPE fortran routine directly from NCL.</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 4, 2022 at 3:26 PM Dennis Shea 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 dir="ltr"><div>NCL developers did NOT develop the underlying fortran codes invoked by 'wrf_cape_3d' <br></div><div>NCL developers created only the C interfaces between NCL and the fortran subroutines.<br></div><div><br></div><div>Hence, we can not answer your questions. I have attached the fortran code that I think is used by NCL.<br></div><div><br></div><div>Good Luck. <br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 4, 2022 at 2:31 PM Sai Prasanth <<a href="mailto:prashanth.bhalachandran@gmail.com" target="_blank">prashanth.bhalachandran@gmail.com</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 style="font-family:garamond,"times new roman",serif">Dear NCL experts, <br></div><div style="font-family:garamond,"times new roman",serif">I am trying to compute the CAPE for the various columns in a cyclonic storm and I'd like to change the definition of the "environment". As a result, I'd like to have control over the various steps in the process and prefer not to directly use the CAPE_3D subroutine. Unfortunately, I can't seem to figure out how to calculate the integral limits viz. the level of free convection and the level of equilibrium. I see that the parcel is defined as the 500 m vertical block within which the various parameters are averaged. Can you please share the code of wrf_cape_3d or explain in algorithmic terms how the environment is defined, and how the two levels that form the integral bounds are computed? <br></div><div style="font-family:garamond,"times new roman",serif"><br></div><div style="font-family:garamond,"times new roman",serif">Thank you, <br></div><div style="font-family:garamond,"times new roman",serif">Prashanth</div></div></blockquote></div></blockquote></div></div>
</div></blockquote></div></blockquote></div></div>
</blockquote></div>