[Pyaos] CAPE and CIN calculation in Python

Thanks Daniel. Your CAPE calculation approach looks robust. I will give it a try with my sounding data. In the mean time, do you plan on putting a function for CIN calculation?


On Sun, Mar 30, 2014 at 10:37 AM, Daniel Rothenberg <darothen@…1166…> wrote:

Hello Gokhan,

I’ve written my own SkewT/LogP code with some lifted parcel calculations - it’s included in a GitHub repo I use to maintain a collection of utilities I’ve written for analyzing output from a cloud resolving model I use. You can find my code for this purpose at https://github.com/darothen/crm-tools/blob/master/vis/soundings.py, along with other useful scripts.

  • Daniel Rothenberg


On Sat, Mar 29, 2014 at 6:32 PM, Gökhan Sever <gokhansever@…287…> wrote:


Lately, I am working on plotting sounding profiles on a SkewT/LogP diagram. The SkewT package which is located at https://github.com/tchubb/SkewT has a nice feature to lift a parcel on dry/moist adiabats. This is very useful to demonstrate the regions of CIN and CAPE overlaid with the full sounding.

However, the package misses these diagnostic calculations. This is the only step holding me back to use Python only (migrating from NCL) for my plotting tasks. I am aware that these calculations are usually performed in fortran. Are there any routines wrapped in Python to calculate CAPE and CIN parameters? Any suggestions or comments would be really appreciated.


Pyaos mailing list