LMAX. The index corresponds to L*(L+1)/2 + M +
                 1, which can be calculated by a call to <B>PlmIndex</B>.

       LMAX      (input) INTEGER

                 The maximum degree of the associated Legendre functions to be
                 computed. If LMAX is -1, allocated memory will be deallocated
                 (see notes below).

       Z         (input) REAL*8

                 The argument of the associated Legendre functions.

       CSPHASE   (input) INTEGER, OPTIONAL

               