Spectral spherical harmonics discrete ordinate method