The PythonicDISORT package is a Discrete Ordinates Solver for the (1D) Radiative Transfer Equation in a plane-parallel, horizontally homogeneous atmosphere. It is coded entirely in Python 3 and is a reimplementation instead of a wrapper. While PythonicDISORT has been optimized for speed, it will naturally be slower than similar FORTRAN algorithms. On the other hand, PythonicDISORT should be easier to install, use, and modify than FORTRAN-based Discrete Ordinates Solvers.
PythonicDISORT is based on Stamnes' FORTRAN DISORT (see References, in particular [2, 3, 8]) and has its main features: multi-layer solver, delta-M scaling, Nakajima-Tanaka (NT) corrections, only flux option, direct beam source, isotropic internal source (blackbody emission), Dirichlet boundary conditions (diffuse flux boundary sources), Bi-Directional Reflectance Function (BDRF) for surface reflection, and more. In addition, we added a subroutine to calculate actinic fluxes to satisfy a user request and we are open to further feature requests as well as feedback.
You may contact me, Dion, through dh3065@columbia.edu.
Our GitHub repository: https://github.com/LDEO-CREW/Pythonic-DISORT also includes our F2PY-wrapped Stamnes DISORT (version 4.0.99) in the disort4.0.99_f2py directory.
The original was downloaded from http://www.rtatmocn.com/disort/. Our wrapper is inspired by https://github.com/kconnour/pyRT_DISORT.
https://pythonic-disort.readthedocs.io/en/latest/
Also see the accompanying Jupyter Notebook Pythonic-DISORT.ipynb in the docs directory
of our GitHub repository.
The Jupyter Notebook provides comprehensive documentation, suggested inputs, explanations,
mathematical derivations and verification tests.
We highly recommend reading the non-optional parts of sections 1 and 2 before use.
Not only do we have verification tests in the Jupyter Notebook, we have also recreated most of the test problems in Stamnes' disotest.f90;
disotest.f90 is included in the disort4.0.99_f2py directory of our GitHub repository. In these tests, the solutions from PythonicDISORT are compared
against solutions from our F2PY-wrapped Stamnes' DISORT (version 4.0.99). With PyTest installed, execute the console command pytest
in the pydisotest directory to run these tests. The pydisotest directory also contains Jupyter Notebooks, one for each test,
to show how the tests were implemented. These notebooks double up as examples of how to use PythonicDISORT.
- From PyPI:
pip install PythonicDISORT - From Conda-forge: (TODO: we need to first publish on Conda-forge)
- By cloning repository:
pip install .in thePythonic-DISORTdirectory;pip install -r all_optional_dependencies.txtto install all optional dependencies (see Requirements to run PythonicDISORT)
- Python 3.8+
numpy >= 1.8.0scipy >= 1.8.0- (OPTIONAL)
pytest >= 6.2.5(Required to use the commandpytest, see PyTest and examples of how to use PythonicDISORT)
autograd >= 1.5jupyter > 1.0.0notebook > 6.5.2
In addition, our F2PY-wrapped Stamnes' DISORT (in the disort4.0.99_f2py directory), or equivalent,
must be set up to run the last section (section 6).
The PythonicDISORT package should be system agnostic given its minimal dependencies and pure Python code. We do not guarantee that our Jupyter Notebook and F2PY-wrapped Stamnes' DISORT will work on other systems though. The latter will almost certainly need user edits to work and in any case it requires FORTRAN compilers which are not included in our GitHub repository. Everything in the repository was built and tested on Windows 11.
I acknowledge funding from NSF through the Learning the Earth with Artificial intelligence and Physics (LEAP) Science and Technology Center (STC) (Award #2019625) under which this package was initially created.
-
S. Chandrasekhar. 1960. Radiative Transfer.
-
Knut Stamnes and S-Chee Tsay and Warren Wiscombe and Kolf Jayaweera. 1988. Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media. http://opg.optica.org/ao/abstract.cfm?URI=ao-27-12-2502.
-
Stamnes, S.. 1999. LLLab disort website. http://www.rtatmocn.com/disort/.
-
Knut Stamnes and Paul Conklin. 1984. A new multi-layer discrete ordinate approach to radiative transfer in vertically inhomogeneous atmospheres. https://www.sciencedirect.com/science/article/pii/0022407384900311.
-
W. J. Wiscombe. 1977. The Delta–M Method: Rapid Yet Accurate Radiative Flux Calculations for Strongly Asymmetric Phase Functions. https://journals.ametsoc.org/view/journals/atsc/34/9/1520-0469_1977_034_1408_tdmrya_2_0_co_2.xml.
-
J. H. Joseph and W. J. Wiscombe and J. A. Weinman. 1976. The Delta-Eddington Approximation for Radiative Flux Transfer. https://journals.ametsoc.org/view/journals/atsc/33/12/1520-0469_1976_033_2452_tdeafr_2_0_co_2.xml.
-
Sykes, J. B.. 1951. Approximate Integration of the Equation of Transfer. https://doi.org/10.1093/mnras/111.4.377.
-
Stamnes, Knut and Tsay, Si-Chee and Wiscombe, Warren and Laszlo, Istvan and Einaudi, Franco. 2000. General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT.
-
Z. Lin and S. Stamnes and Z. Jin and I. Laszlo and S.-C. Tsay and W.J. Wiscombe and K. Stamnes. 2015. Improved discrete ordinate solutions in the presence of an anisotropically reflecting lower boundary: Upgrades of the DISORT computational tool. https://www.sciencedirect.com/science/article/pii/S0022407315000679.
-
Trefethen, L. N.. 1996. Finite difference and spectral methods for ordinary and partial differential equations. https://people.maths.ox.ac.uk/trefethen/pdetext.html.
-
T. Nakajima and M. Tanaka. 1988. Algorithms for radiative intensity calculations in moderately thick atmospheres using a truncation approximation. https://www.sciencedirect.com/science/article/pii/0022407388900313.
-
Connour, Kyle and Wolff, Michael. 2020. pyRT_DISORT: A pre-processing front-end to help make DISORT simulations easier in Python. https://github.com/kconnour/pyRT_DISORT.
