[ncl-talk] DCAPE scripts?

[Kacie Shourd]

Hi Andy & Dennis --

Andy, thanks so much for sharing this script!! It is quite helpful indeed.

Dennis, I'll also play around with the script you've included! Thank you so much!

Have a great weekend,
Kacie

---
Kacie N. Shourd
드림

Graduate Research Assistant
Desert Research Institute
Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu>
Work: (775) 674-7049


________________________________
From: Dennis Shea [shea at ucar.edu]
Sent: Saturday, August 19, 2017 9:55 AM
To: Andrew Penny
Cc: Kacie Shourd; ncl-talk at ucar.edu
Subject: Re: [ncl-talk] DCAPE scripts?

Kacie & Andy ....

Andy ... THX!

A very minor comment on Andy's script.

  delete(ind_cape)
  delete(rev_PA)
  delete(pseudo_PA)
  delete(rev_NA)
  delete(pseudo_NA)
  delete(rev_CAPE)
  delete(pseudo_CAPE)
  delete(DCAPE)

can be consolidated into one line  (personal preference)

   delete( [/ ind_cape, rev_PA, pseudo_PA, rev_NA, pseudo_NA, rev_CAPE, pseudo_CAPE, DCAPE/] )

====
I looked at my CAPE archives. I have a number of f77 CAPE codes from 1997 onward !!!
Also, one IDL code which I have attached 'for fun'  !!!
====
I have attached an untested version 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


%> WRAPIT cape_sub_ke.f
COMPILING cape_sub_ke.f
LINKING
END WRAPIT

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).

   cape_results = cape_ke(p, t, q)    ; p(levl) or p(time,lev,lat,lon), t/q (time,lev,lat,lon)

   cape_results would contain cape_pos(lev), cape_neg, dcape(lev)  and the integrated (vertically weighted) sums.

Alas, the best laid plans did not result in any function being added.  :-(

=====
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.

Cheers
D






On Sat, Aug 19, 2017 at 9:38 AM, Andrew Penny <andybpenny at gmail.com<mailto:andybpenny at gmail.com>> wrote:
Hi Kacie,

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:

ftp://texmex.mit.edu/pub/emanuel/BOOK/

I haven’t used the code in a while, but hopefully it will be helpful.

Andy

--
Andrew Penny
Hurricane Model Diagnostician
SRG / National Hurricane Center
11691 SW 17th Street
Miami, FL 33165-2149
phone: 305.229.4457<tel:(305)%20229-4457>
email: andrew.penny at noaa.gov<mailto:andrew.penny at noaa.gov>




Hi everyone!

I'm just wondering if anyone has a NCL script to calculate DCAPE that they would be willing to share?

Thanks!

Best
Kacie

---
Kacie N. Shourd
??

Graduate Research Assistant
Desert Research Institute
Kacie.Shourd at dri.edu<mailto:Kacie.Shourd at dri.edu><mailto:Kacie.Shourd at dri.edu>
Work: (775) 674-7049<tel:(775)%20674-7049>

_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk



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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170819/870a59ba/attachment.html