Abstract: The recently developed matrix method to solve the discrete ordinate approximation to the radiative transfer equation in plane parallel geometry (Stamnes and Swanson, 1981) is extended to compute the full azimuthal dependence ...