Abstract
The difficulties inherent in the conventional numerical implementation of the discrete ordinate method (following Chandrasekhar's prescription) for solving the radiative transfer equation are discussed. A matrix formulation is developed to overcome these difficulties, and it is specifically shown that the order of the algebraic eigenvalue problem can be reduced by a factor of 2. An expression for the source function is derived and used to obtain angular distributions. By appealing to the reciprocity principle, it is shown that substantial computational shortcuts are possible if only integrated quantities such as albedo and transmissivity are required. Comparison of fluxes calculated by the present approach with those obtained by other methods shows that low-order discrete ordinate approximations yield very accurate results. Thus, the present approach offers an efficient and reliable computational scheme that lends itself readily to the solution of a variety of radiative transfer problems in realistic planetary atmospheres.
