Extension of the discrete-ordinate algorithm and efficient radiative transfer calculation

As an accurate and efficient algorithm, the discrete-ordinate method (DOM) has been used to solve the radiative transfer problem of plane-parallel scattering atmosphere illuminated by a parallel beam, an idealized case of the sun, from above the atmosphere. In this paper, we extend this algorithm so that radiative problems of more general sources, such as parallel surface sources that illuminate with a parallel beam in any direction from any vertical position, and general surface sources that illuminate continuously in a hemisphere, can be solved. For a problem where intensity distributions are sought for a number of different sources within the same atmosphere-surface system, the intrinsic properties of DOM are used so that the time required for the solution for extra sources is reduced to a substantially small amount. In the case of parallel surface sources, numerical testing has shown that the amount can be reduced to as little as 15% of a full solution. Examples of applications are presented.