<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=ks_c_5601-1987">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Andy & Dennis --
<div><br>
</div>
<div>Andy, thanks so much for sharing this script!! It is quite helpful indeed.</div>
<div><br>
</div>
<div>Dennis, I'll also play around with the script you've included! Thank you so much!</div>
<div><br>
</div>
<div>Have a great weekend,</div>
<div>Kacie</div>
<div><br>
</div>
<div>---</div>
<div>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Courier New";">Kacie N. Shourd</span><span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"></span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Malgun Gothic", sans-serif;">µå¸²</span><span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"></span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"> </span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Courier New"; color: rgb(34, 34, 34);">Graduate Research Assistant</span><span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"></span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Courier New"; color: rgb(34, 34, 34);">Desert Research Institute</span><span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"></span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Courier New"; color: rgb(34, 34, 34);"><a href="mailto:Kacie.Shourd@dri.edu" target="_blank" style="color: rgb(149, 79, 114);"><span style="color: blue;">Kacie.Shourd@dri.edu</span></a></span><span style="font-size: 9.5pt; font-family: Arial, sans-serif; color: rgb(34, 34, 34);"></span></p>
<p class="MsoNormal" style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background: white;">
<span style="font-size: 10pt; font-family: "Courier New"; color: rgb(34, 34, 34);">Work: (775) 674-7049</span></p>
</div>
<div><br>
</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF391094" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> Dennis Shea [shea@ucar.edu]<br>
<b>Sent:</b> Saturday, August 19, 2017 9:55 AM<br>
<b>To:</b> Andrew Penny<br>
<b>Cc:</b> Kacie Shourd; ncl-talk@ucar.edu<br>
<b>Subject:</b> Re: [ncl-talk] DCAPE scripts?<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>Kacie & Andy ....<br>
<br>
</div>
Andy ... THX!<br>
<br>
</div>
A <b>very minor</b> comment on Andy's script.<br>
<br>
delete(ind_cape)<br>
delete(rev_PA)<br>
delete(pseudo_PA)<br>
delete(rev_NA)<br>
delete(pseudo_NA)<br>
delete(rev_CAPE)<br>
delete(pseudo_CAPE)<br>
delete(DCAPE) <br>
<br>
</div>
can be consolidated into one line (personal preference)<br>
<br>
</div>
delete( [/ ind_cape, rev_PA, pseudo_PA, rev_NA, pseudo_NA, rev_CAPE, pseudo_CAPE, DCAPE/] )<br>
<br>
====<br>
</div>
I looked at my CAPE archives. I have a number of f77 CAPE codes from 1997 onward !!!
<br>
Also, one IDL code which I have attached 'for fun' !!! <br>
====<br>
I have attached an <b>untested version</b> of Kerry Emmanual's code (dated 21 Jan 1997) which I was going to include in NCL as a built-in function. I am not sure why I did not do it!!! I made a minor change and it compiles under WRAPIT<br>
</div>
<br>
</div>
<div><br>
%> WRAPIT cape_sub_ke.f<br>
COMPILING cape_sub_ke.f<br>
LINKING<br>
END WRAPIT<br>
<br>
</div>
<div>The idea was that this would have a hidden 'C-driver' (written by NCL developers) that would work on one (lev) or more individual soundings (ntim,lev) or (nstation,lev) or gridded data (time,lev,lat,lon).<br>
<br>
</div>
<div> cape_results = cape_ke(p, t, q) ; p(levl) or p(time,lev,lat,lon), t/q (time,lev,lat,lon)<br>
<br>
</div>
<div> cape_results would contain cape_pos(lev), cape_neg, dcape(lev) and the integrated (vertically weighted) sums.<br>
<br>
</div>
<div>Alas, the best laid plans did not result in any function being added. :-(<br>
</div>
<div><br>
</div>
<div>=====<br>
</div>
As I recall, results from several different CAPE codes gave very different answers. My recollection is that CAPE is quite sensitive to where the parcel ascent is started.<br>
<br>
</div>
Cheers <br>
</div>
D<br>
<div>
<div><br>
<br>
<div><br>
<div><br>
<br>
</div>
</div>
</div>
</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Sat, Aug 19, 2017 at 9:38 AM, Andrew Penny <span dir="ltr">
<<a href="mailto:andybpenny@gmail.com" target="_blank">andybpenny@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 style="word-wrap:break-word">
<div>Hi Kacie,</div>
<div><br>
</div>
<div>Attached is a script with a procedure (¡°computecape") that calls a fortran program to compute DCAPE (among other things). The fortran is adapted from Kerry Emanuel¡¯s code that is available here: </div>
<div><br>
</div>
<div><a href="ftp://texmex.mit.edu/pub/emanuel/BOOK/" target="_blank">ftp://texmex.mit.edu/pub/<wbr>emanuel/BOOK/</a></div>
<div><br>
</div>
<div>I haven¡¯t used the code in a while, but hopefully it will be helpful.</div>
<div><br>
</div>
<div>Andy</div>
<span class="HOEnZb"><font color="#888888">
<div><br>
</div>
<div>-- <br>
<div class="m_-1853359923273109499gmail_signature">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"><span style="font-family:arial,helvetica,sans-serif"><font size="2"><span style="color:rgb(102,102,102)"><span style="background-color:rgb(255,255,255)">Andrew Penny<br>
Hurricane Model Diagnostician<br>
SRG / National Hurricane Center<br>
11691 SW 17th Street<br>
Miami, FL 33165-2149<br>
</span></span></font></span></div>
<div><span style="font-family:arial,helvetica,sans-serif"><font size="2"><span style="color:rgb(102,102,102)"><span style="background-color:rgb(255,255,255)">phone:
<a href="tel:(305)%20229-4457" value="+13052294457" target="_blank">305.229.4457</a><br>
</span></span></font></span></div>
<div dir="ltr"><span style="font-family:arial,helvetica,sans-serif"><font size="2"><span style="color:rgb(102,102,102)"><span style="background-color:rgb(255,255,255)">email: <a href="mailto:andrew.penny@noaa.gov" target="_blank">andrew.penny@noaa.gov</a></span></span></font></span></div>
<div dir="ltr"><br>
</div>
<div dir="ltr"></div>
</div>
</div>
</div>
</div>
</font></span></div>
<br>
<div style="word-wrap:break-word">
<div>
<div class="m_6745979242578083254gmail_signature">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"></div>
</div>
</div>
</div>
</div>
</div>
<br>
<div style="word-wrap:break-word">
<div>
<div class="m_6612807224869386573gmail_signature">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr"></div>
</div>
</div>
</div>
</div>
<div><br>
</div>
<div>
<blockquote type="cite">Hi everyone!<br>
<br>
I'm just wondering if anyone has a NCL script to calculate DCAPE that they would be willing to share?<br>
<br>
Thanks!<br>
<br>
Best<br>
Kacie<br>
<br>
---<br>
Kacie N. Shourd<br>
??<br>
<br>
Graduate Research Assistant<br>
Desert Research Institute<br>
<a href="mailto:Kacie.Shourd@dri.edu" target="_blank">Kacie.Shourd@dri.edu</a><<a href="mailto:Kacie.Shourd@dri.edu" target="_blank">mailto:<wbr>Kacie.Shourd@dri.edu</a>><br>
Work: <a href="tel:(775)%20674-7049" value="+17756747049" target="_blank">(775) 674-7049</a></blockquote>
</div>
</div>
<br>
______________________________<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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
<br>
<font face="Arial" size="2">PUBLIC RECORDS NOTICE: In accordance with NRS Chapter 239, this email and responses, unless otherwise made confidential by law, may be subject to the Nevada Public Records laws and may be disclosed to the public upon request.</font>
</body>
</html>