<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" data-smartmail="gmail_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" data-smartmail="gmail_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" data-smartmail="gmail_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">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>