SCIAPOL:discrete ordinates FORTRAN codefor computing Stokes parameters of light
reflected and transmitted by a finite homogeneous slab with a given optical
thickness, single scattering albedo, and phase
matrix: