PyMieDAP: a Python--Fortran tool to compute fluxes and polarization signals of (exo)planets

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 sun- or starlight that is reflected by, respectively, Solar System planets or moons, or exoplanets at a range of wavelengths. The radiative transfer computations are based on an adding--doubling 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 code's principal functionalities. And we 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