      complex normalization will be used where the integral of the squared
       Legendre functions over the interval [-1, 1] is 2/(2L+1). The default
       is to exclude the Condon-Shortley phase, but this can be modified by
       setting the optional argument CSPHASE to -1.


</PRE>
<H2 class="man">ARGUMENTS</H2 class="man"><PRE>
       P         (output) REAL*8, DIMENSION ((LMAX+1)*(LMAX+2)/2)

                 An array of Schmidt-normalized associated Legendre functions
                 up to degree 