J/A+A/616/A147 Python Mie Doubling-Adding Programme (Rossi+, 2018)
PyMieDAP : a Python-Fortran tool for computing fluxes and polarization
signals of (exo)planets.
Rossi L., Berzosa-Molina J., Stam D.M.
<Astron. Astrophys. 616, A147 (2018)>
=2018A&A...616A.147R 2018A&A...616A.147R (SIMBAD/NED BibCode)
ADC_Keywords: Models ; Exoplanets ; Polarization
Keywords: planets and satellites: atmospheres - polarization -
radiative transfer
Abstract:
PyMieDAP (the Python Mie Doubling-Adding Programme) is a Python-based
tool for computing the total linearly and circularly polarized fluxes
of incident unpolarized sunlight or starlight that is reflected by
solar system planets or moons, respectively, or by exoplanets at a
range of wavelengths. The radiative transfer computations are based on
an doubling-adding Fortran algorithm and fully include polarization
for all orders of scattering. The model (exo)planets are described by
a model atmosphere composed of a stack of homogeneous layers
containing gas and/or aerosol and/or cloud particles bounded below by
an isotropically depolarizing surface (that is optionally black). The
reflected light can be computed spatially resolved and/or
disk-integrated. Spatially resolved signals are mostly representative
for observations of solar system planets (or moons), while
disk-integrated signals are mostly representative for exoplanet
observations. PyMieDAP is modular and flexible, and allows users to
adapt and optimize the code according to their needs. PyMieDAP keeps
options open for connections with external programs and for future
additions and extensions. In this paper, we describe the radiative
transfer algorithm that PyMieDAP is based on and the principal
functionalities of the code. We also provide benchmark results of
PyMieDAP that can be used for testing its installation and for
comparison with other codes. PyMieDAP is available online under the
GNU GPL license at http://gitlab.com/loic.cg.rossi/pymiedap
Description:
These tables are showing benchmark tests to compare results of the
PyMieDAP programme (http://gitlab.com/loic.cg.rossi/pymiedap) with
similar doubling adding method from de Haan et al. (1987, A&A, 183,
371). Users can use them to check if their installation of Pymiedap is
correct.
File Summary:
--------------------------------------------------------------------------------
FileName Lrecl Records Explanations
--------------------------------------------------------------------------------
ReadMe 80 . This file
tablea1.dat 88 12 Locally reflected Stokes elements for model 1
described in the paper
tablea2.dat 93 12 Locally reflected Stokes elements for model 2
described in the paper
tablea3.dat 35 37 Disk-integrated flux and polarization for a
clear sky planet
tablea4.dat 35 37 Disk-integrated flux and polarization for an
atmosphere with aerosols D
--------------------------------------------------------------------------------
See also:
http://gitlab.com/loic.cg.rossi/pymieda : PyMieDAP programme
Byte-by-byte Description of file: tablea1.dat tablea2.dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1- 3 F3.1 --- mu0 Cosine of the solar zenith angle (mu_0)
5- 7 F3.1 --- mu Cosine of the viewing angle (mu)
9- 12 F4.1 deg phi-phi0 Azimuth angle (phi-phi_0)
14- 21 F8.6 W/m+2 IdeHaan Local stokes element I, as computed by
de Haan et al. (1987, A&A, 183, 371)
(IdeHaan)
23- 30 F8.6 W/m+2 Ipymiedap Local stokes element I, as computed by
Pymiedap (I_pymiedap)
32- 40 F9.6 W/m+2 QdeHaan Local stokes element Q, as computed by
de Haan et al. (1987, A&A, 183, 371)
(QdeHaan)
42- 50 F9.6 W/m+2 Qpymiedap Local stokes element Q, as computed by
Pymiedap (Q_pymiedap)
52- 60 F9.6 W/m+2 UdeHaan Local stokes element U, as computed by
de Haan et al. (1987, A&A, 183, 371)
(UdeHaan)
62- 70 F9.6 W/m+2 Upymiedap Local stokes element U, as computed by
Pymiedap (U_pymiedap)
72- 79 F8.6 W/m+2 VdeHaan Local stokes element V, as computed by
de Haan et al. (1987, A&A, 183, 371)
(VdeHaan)
81- 93 F13.11 W/m+2 Vpymiedap Local stokes element V, as computed by
Pymiedap (V_pymiedap)
--------------------------------------------------------------------------------
Byte-by-byte Description of file: tablea3.dat tablea4.dat
--------------------------------------------------------------------------------
Bytes Format Units Label Explanations
--------------------------------------------------------------------------------
1- 5 F5.1 deg alpha Planetary phase angle (alpha)
7- 12 F6.4 --- FStam Reflected flux for Stam et al.
2006A&A...452..669S 2006A&A...452..669S, normalized to planet
geometric albedo (F_Stam)
14- 19 F6.4 --- Fpymiedap Reflected flux for PyMieDAP normalized to
planet geometric albedo (F_pymiedap)
21- 27 F7.4 --- PsStam Reflected degree of polarization for
Stam et al. 2006A&A...452..669S 2006A&A...452..669S (Ps_Stam)
29- 35 F7.4 --- Pspymiedap Reflected degree polarization for PyMieDAP
(Ps_pymiedap)
--------------------------------------------------------------------------------
Acknowledgements:
Loic Rossi, l.c.g.rossi(at)tudelft.nl
(End) Patricia Vannier [CDS] 05-Jun-2018