ccnact 0.0.5

KISS SciPy-based CCN activation model

☁️ KISS SciPy-based CCN activation model

📌 overview

ccnact is a simple, yet complete, adiabatic/hydrostatic air-parcel framework employing moving-sectional/particle-resolved aerosol-cloud microphysics, featuring:

• 🧮 integration using SciPy's LSODA interface
• 📝 ODE system based on Arabas & Shima 2017 (extended to polydisperse aerosol size spectrum)
• 🏁 wet radii equilibration for input dry-size distribution using SciPy's elementwise scalar optimisation
• 🌪️ capable of resolving aerosol activation, deactivation, drop growth, evaporation and ripening
• ⚙️ single-function interface allowing to modify every single constants, and returning a tuple of:
  • concentration of activated droplets (at STP)
  • maximal supersaturation
• 🧩 effective interfacing options for Matlab, Julia, etc
• 📈 mulit-modal lognormal (using SciPy's stats routines) spectrum specification (with concentration interpretted as at STP)
• ⚖️ implemeted using Pint dimensional analysis (physical units consistency checks) enabled for tests only
• 🚀 subsecond execution times for common parameter settings
• 🔗 KISS design: SciPy, NumPy & Pint are the only dependencies; single ~500 LOC file (physics + setup + tests)

The last four points motivated the development of this package - the project originated from a search for a simple, lightweight (in dependencies) and fast CCN activation air-parcel model with concise code, automated testing and no hardcoded constants.

💻 notes for users

To install the package, try: pip install ccnact

Using from Python:

from ccnact import parcel
help(parcel)
n_act, s_max = parcel(...)

Interfacing from Matlab (using the built-in Python bridge):

ccnact = py.importlib.import_module('ccnact');
ccnact.parcel(pyargs(...
   'MAC', 1,...
   'n_bins', int32(100),...
   'p', 101300,...
   'T', 300,...
   'RH', .99,...
   'dt', 1,...
   'nt', int32(100),...
   'w', 2,...
   'sigma', 0.072,...
   'kappa', py.tuple({1}),...
   'meanr', py.tuple({3e-8}),...
   'gstdv', py.tuple({1.5}),...
   'n_tot', py.tuple({1e9})...
))

⚙ notes for developers

To execute the tests after checking out from git: pip install -e .[dev]; pytest ccnact.py

To set-up pre-commit: pip install pre-commit; pre-commit install