     polynomials over the interval [-1, 1] is 2/(2L+1). Note that the
       derivative of the Legendre polynomials is calculated with respect to
       its arguement Z, and not latitude or colatitude. If Z=cos(theta), where
       theta is the colatitude, then it is only necessary to multiply DP by
       -sin(theta) to obtain the derivative with respect to theta.


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

                 An array of Schm