This is an implementation of a simple, adiabatic cloud parcel model for use in aerosol-cloud interaction studies. It is based on the model used by Nenes et al (2001), but with several key modifications:
- Implementation of
$\kappa$ -Kohler theory for condensation physics (Petters and Kreidenweis, 2007) - Extension of model to handle arbitrary sectional representations of aerosol populations, based on user-controlled empirical or parameterized size distributions
- Improved, modular numerical framework for integrating the model, including bindings to several different stiff integrators:
-
lsoda- scipy ODEINT wrapper -
vode, lsode*, lsoda*- ODEPACK via odespy -
cvode- SUNDIALS via Assimulo
among other details. It also includes a library of droplet activation routines and scripts/notebooks for evaluating those schemes against equivalent calculations done with the parcel model.
Updated code can be found the project github repository. If you'd like to use this code or have any questions about it, please contact the author. In particular, if you use this code for research purposes, be sure to carefully read through the model and ensure that you have tweaked/configured it for your purposes (i.e., modifying the accomodation coefficient); other derived quantities).
Detailed documentation is available, including a scientific description, installation details, and a basic example which produces a figure like the plot at the top of this page.
Required
Optional
The following packages are used for better numerics (ODE solving), handling data in memory, or saving output.
The easiest way to satisfy the basic requirements for building and running the model is to use the Anaconda scientific Python distribution.
http://github.com/darothen/parcel_model
Please fork this repository if you intend to develop the model further so that the code's provenance can be maintained.
All scientific code should be licensed. This code is released under the New BSD (3-clause) license.