[ncl-talk] DCAPE scripts?

[Dennis Shea]

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> 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 <(305)%20229-4457>
> email: 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 <Kacie.Shourd at dri.edu>>
> Work: (775) 674-7049
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170819/3bbcb7e9/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cape_sound.pro
Type: application/octet-stream
Size: 5585 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170819/3bbcb7e9/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cape_sub_ke.f
Type: application/octet-stream
Size: 13853 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170819/3bbcb7e9/attachment-0001.obj