A NEW LOOK AT THE DISCRETE ORDINATE METHOD FOR RADIATIVE-TRANSFER CALCULATIONS IN ANISOTROPICALLY SCATTERING ATMOSPHERES .2. INTENSITY COMPUTATIONS

The recently developed matrix method to solve the discrete ordinate approximation to the radiative transfer equation in plane parallel geometry is extended to compute the full azimuthal dependence of the intensity. Comparing computed intensities with those obtained by other established methods, it is found that for phase functions typical of atmospheric haze 32 streams are sufficient for better than 1% agreement, while 16 streams yield an accuracy of about 1-5% except for angles close to the forward and backward directions for which the error is about 10-15%. The results of the intensity computations are summarized by presenting three-dimensional ″stack plots″ of the intensity as a function of polar and azimuthal angles. It is shown that for flux calculations four streams suffice to obtain 1% accuracy, while eight streams yield an accuracy better than 0. 1%.