Setup:
$ P = \begin{pmatrix} P1 & P2 & 0 & 0 \\ P2 & P1 & 0 & 0 \\ 0 & 0 & P3 & P4 \\ 0 & 0 & -P4 & P3 \end{pmatrix} $
The corresponding Legendre polynomials (“pmom”) and the single scattering albedo (“ssa”) are also included in the file.
The expansion moments over generalized spherical functions for the phase matrix have been provided by S. Korkin and are available here.
Output format: