TY - JOUR
T1 - Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media
AU - Stamnes, Knut
AU - Tsay, S. Chee
AU - Wiscombe, Warren
AU - Jayaweera, Kolf
PY - 1988/6
Y1 - 1988/6
N2 - We summarize an advanced, thoroughly documented, and quite general purpose discrete ordinate algorithm for time-independent transfer calculations in vertically inhomogeneous, nonisothermal, plane-parallel media. Atmospheric applications ranging from the UV to the radar region of the electromagnetic spectrum are possible. The physical processes included are thermal emission, scattering, absorption, and bidirectional reflection and emission at the lower boundary. The medium may be forced at the top boundary by parallel or diffuse radiation and by internal and boundary thermal sources as well. We provide a brief account of the theoretical basis as well as a discussion of the numerical implementation of the theory. The recent advances made by ourselves and our collaborators-advances in both formulation and numerical solution-are all incorporated in the algorithm. Prominent among these advances are the complete conquest of two illconditioning problems which afflicted all previous discrete ordinate implementations: (1) the computation of eigenvalues and eigenvectors and (2) the inversion of the matrix determining the constants of integration. Copies of the FORTRAN program on microcomputer diskettes are available for interested users.
AB - We summarize an advanced, thoroughly documented, and quite general purpose discrete ordinate algorithm for time-independent transfer calculations in vertically inhomogeneous, nonisothermal, plane-parallel media. Atmospheric applications ranging from the UV to the radar region of the electromagnetic spectrum are possible. The physical processes included are thermal emission, scattering, absorption, and bidirectional reflection and emission at the lower boundary. The medium may be forced at the top boundary by parallel or diffuse radiation and by internal and boundary thermal sources as well. We provide a brief account of the theoretical basis as well as a discussion of the numerical implementation of the theory. The recent advances made by ourselves and our collaborators-advances in both formulation and numerical solution-are all incorporated in the algorithm. Prominent among these advances are the complete conquest of two illconditioning problems which afflicted all previous discrete ordinate implementations: (1) the computation of eigenvalues and eigenvectors and (2) the inversion of the matrix determining the constants of integration. Copies of the FORTRAN program on microcomputer diskettes are available for interested users.
UR - http://www.scopus.com/inward/record.url?scp=84966148288&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84966148288&partnerID=8YFLogxK
U2 - 10.1364/AO.27.002502
DO - 10.1364/AO.27.002502
M3 - Article
AN - SCOPUS:84966148288
SN - 1559-128X
VL - 27
SP - 2502
EP - 2509
JO - Applied Optics
JF - Applied Optics
IS - 12
ER -