Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media.

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.

[1]  J. Sykes,et al.  Approximate integration of the equation of transger , 1951 .

[2]  H. V. Hulst.  Multiple scattering in planetary atmospheres , 1971 .

[3]  Kuo-Nan Liou,et al.  A Numerical Experiment on Chandrasekhar's Discrete-Ordinate Method for Radiative Transfer: Applications to Cloudy and Hazy Atmospheres , 1973 .

[4]  J. Hansen,et al.  Light scattering in planetary atmospheres , 1974 .

[5]  Warren J. Wiscombe,et al.  Extension of the doubling method to inhomogeneous sources , 1976 .

[6]  W. Wiscombe.  The Delta–M Method: Rapid Yet Accurate Radiative Flux Calculations for Strongly Asymmetric Phase Functions , 1977 .

[7]  C. Bohren,et al.  An introduction to atmospheric radiation , 1981 .

[8]  V. Klema.  LINPACK user's guide , 1980 .

[9]  Knut Stamnes,et al.  A New Look at the Discrete Ordinate Method for Radiative Transfer Calculations in Anisotropically Scattering Atmospheres , 1981 .

[10]  K. Stamnes,et al.  A New Look at the Discrete Ordinate Method for Radiative Transfer Calculations in Anisotropically Scattering Atmospheres. II: Intensity Computations , 1981 .

[11]  Knut Stamnes,et al.  On the computation of angular distributions of radiation in planetary atmospheres , 1982 .

[12]  Knut Stamnes,et al.  Reflection and transmission by a vertically inhomogeneous planetary atmosphere , 1982 .

[13]  W. Wiscombe.  Atmospheric radiation: 1975–1983 , 1983 .

[14]  Knut Stamnes,et al.  A new multi-layer discrete ordinate approach to radiative transfer in vertically inhomogeneous atmospheres , 1984 .

[15]  Wayne R. Cowell,et al.  Sources and development of mathematical software , 1984 .

[16]  Boris M. Smirnov,et al.  Reference Data on Atoms, Molecules, and Ions , 1985 .

[17]  J. Lenoble.  Radiative transfer in scattering and absorbing atmospheres: Standard computational procedures , 1985 .

[18]  Knut Stamnes,et al.  The theory of multiple scattering of radiation in plane parallel atmospheres , 1986 .

[19]  Teruyuki Nakajima,et al.  Matrix formulations for the transfer of solar radiation in a plane-parallel scattering atmosphere. , 1986 .

[20]  Teruyuki Nakajima,et al.  COMPUTATION OF EIGENVALUES AND EIGENVECTORS FOR THE DISCRETE ORDINATE AND MATRIX OPERATOR METHODS IN RADIATIVE TRANSFER , 1988 .