Python package for Titan''s aerosols models
Project description
Python package for Titan’s aerosols models
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
Note
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)
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
titan-aerosols-0.1.1.tar.gz
(27.1 kB
view hashes)
Built Distribution
Close
Hashes for titan_aerosols-0.1.1-py2.py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | eccc088e65364ab4f4d93870dba779c05ed968f15d05b627d0bc75f24352700b |
|
MD5 | c29c65154648fc7123cf57dbc67c85d2 |
|
BLAKE2b-256 | bde7cdf9462710efd32e6c78743ab8e0e2b580c793327da61b81649dc03b930f |