Install With pip: $ pip install titan-aerosols Or directly from the source files: $ git clone https://github.com/seignovert/python-titan-aerosols.git $ cd python-titan-aerosols ; python setup.py install Python usage >>> from aerosols import index_tholins >>> nr, ni = index_tholins(338e-9) (1.6358611317270018, 0.16996428003261624) >>> from aerosols import mie_tholins >>> wvln = 338e-9 # Wavelength (m) >>> rm = 50e-9 # Monomer radius (m) >>> qsct, qext, qabs, gg, theta, P = mie_tholins(wvln, rm) (2.0166342056786466e-15, 5.220659440190628e-15, 3.2040252345119816e-15, 0.18998404961211593, array([0., 0.01745329, ..., 3.14159265]), array([2.22554999, ..., 0.89675182])) >>> from aerosols import fractals_tholins >>> Df = 2.0 >>> N = 266 >>> qsct, qext, qabs, gg, theta, P = fractals_tholins(wvln, rm, Df, N) (1.6372818590819733e-12, 2.403009993357711e-12, 7.657281342757378e-13, None, array([0, ..., 3.14159265]), array([135.84917753, ..., 0.16136205])) >>> N = 3000 >>> qsct, qext, qabs, gg, theta, P = fractals_tholins(wvln, rm, Df, N) ValueError: Model tested only for N = 2 - 1024 (received N=3000) >>> qsct, qext, qabs, gg, theta, P = fractals_tholins(wvln, rm, Df, N, force=True) (2.0008021769173844e-11, 2.7831577064458474e-11, 7.823555295284632e-12, None, array([0, ..., 3.14159265]), array([1.20379507e+03, ..., 1.28640874e-01])) >>> from aerosols import mie >>> qsct, qext, qabs, gg, theta, P = mie_tholins(wvln, nr, ni, rm) (...) >>> from aerosols import fractals >>> qsct, qext, qabs, gg, theta, P = fractals(wvln, nr, ni, rm, Df, N) (...) CLI usage $ fractal_tholins --help usage: fractal_tholins [-h] [--phase-function] [--nang NANG] [--fractal-dimension FRACTAL_DIMENSION] [--force] wvln rm N Fractals cross-sections and phase function for tholin aggregate. Use default tholins indexes (CVD) and Tomasko et al. 2008. positional arguments: wvln Wavelength (m) rm Monomer radius (m) N Number of monomers optional arguments: -h, --help show this help message and exit --phase-function, -p Display the phase function --nang NANG Number of angles for the phase function (0 -> pi/2) --fractal-dimension FRACTAL_DIMENSION, -df FRACTAL_DIMENSION Fractal dimension --force, -f Bypass validity checks $ fractal_tholins 338e-9 60e-9 266 # Cross sections: Scattering: 2.906e-12 m^-2 Absorption: 1.281e-12 m^-2 Extinction: 4.187e-12 m^-2 $ fractal_tholins -p 338e-9 60e-9 266 # Phase function 0.0 1.86e+02 1.0 1.78e+02 ... 179.0 1.17e-01 180.0 1.17e-01 $ fractal_tholins -p --nang 10 338e-9 60e-9 266 # Phase function 0.0 2.52e+02 10.0 2.36e+01 ... 170.0 1.58e-01 180.0 1.59e-01 $ fractal_tholins -df 2.3 338e-9 60e-9 266 Model tested only for Df = 2 (received Df=2.30) $ fractal_tholins 338e-9 60e-9 266 -df 2.3 --force # Cross sections: Scattering: 2.960e-12 m^-2 Absorption: 1.152e-12 m^-2 Extinction: 4.112e-12 m^-2 This package is an early attempt to model Titan's aerosols scattering based on Tomasko et al. 2008 paper (doi:10.1016/j.pss.2007.11.019)