 L, order M, and I (1 = cosine, 2 = sine) corresponds to the
       index 1 + L**2 + (I-1)*L + M.


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

                    The input 1-D indexed array of real spherical harmonic
                    coefficients.

       CILM         (output) REAL*8, DIMENSION (2, LMAX+1, LMAX+1)

                    The 3-D arrary of output real spherical harmonic
                    coefficients.

       LMAX 