ccnact 0.0.4

KISS SciPy-based CCN activation model

☁️ KISS SciPy-based CCN activation model

📌 overview

• 🧮 integration using SciPy interface to LSODA ODE solver
• 📝 ODE system based on Arabas & Shima 2017 (extended to polydisperse aerosol size spectrum)
• 🌪️ capable of resolving aerosol activation, deactivation, drop growth, evaporation and ripening
• ⚙️ single-function interface allowing to modify all constants, and returning a tuple of:
  • concentration of activated droplets (at STP) &
  • maximal supersaturation
• 📈 mulit-modal lognormal spectrum specification (with concentration interpretted as at STP)
• ⚖️ implemeted using Pint dimensional analysis (physical units consistency checks) enabled for tests only
• 🔗 KISS desing: SciPy, NumPy and Pint are the only dependencies, model+tests wrapped into a single (and short) .py file
• 🚀 subsecond execution times for common parameter settings

💻 notes for users

To install the package, try: pip install git+https://github.com/open-atmos-krk/ccnact.git

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: pip install -e .[dev]; pytest ccnact.py

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