<div dir="ltr"><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><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-width:1px;border-left-style:solid;border-left-color: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-width:1px;border-left-style:solid;border-left-color: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>