Definitions: Real normalized harmonics
Any real squareintegrable function can be expressed as a series of spherical harmonic functions
\begin{equation} f\left(\theta,\phi\right) = \sum_{l=0}^{\infty} \sum_{m=l}^l f_{lm} \, Y_{lm}\left(\theta,\phi \right), \label{eq:f} \end{equation}
where is the spherical harmonic coefficient, is the corresponding spherical harmonic function, is colatitude, is longitude, and and are the spherical harmonic degree and order, respectively. The real spherical harmonics are defined as
where the normalized associated Legendre functions for use with the normalized spherical harmonic functions are given by
and where is the Kronecker delta function. The unnormalized associated Legendre functions are derived from the standard Legendre polynomials using the relations \begin{equation} P_{lm}(\mu) = \left( 1\mu^2\right)^{m/2} \frac{d^m}{d\mu^m} P_l(\mu) \end{equation} and
\begin{equation} P_l(\mu) = \frac{1}{2^l l!}\frac{d^l}{d\mu^l}\left(\mu^21\right)^l. \end{equation}
The normalized associated Legendre functions are orthogonal for a given value of ,
and the spherical harmonic functions are orthogonal for all degrees and orders
\begin{equation} \int_\Omega Y_{lm}(\theta,\phi) \, Y_{l’m’}(\theta,\phi) \, d\Omega = 4 \pi \, \delta_{ll’}\, \delta_{mm’}, \end{equation} where is the differential surface area on the unit sphere, . By multiplying equation \eqref{eq:f} by and integrating over all space, it is straightforward to show that the spherical harmonic coefficients of a function can be calculated by the integral \begin{equation} f_{lm} = \frac{1}{4\pi} \int_\Omega f(\theta,\phi) \, Y_{lm}(\theta,\phi) \, d\Omega. \end{equation}
Power spectrum
Parseval’s theorem in Cartesian geometry relates the integral of a function squared to the sum of the squares of the function’s Fourier coefficients. This relation is easily extended to spherical geometry using the orthogonality properties of the spherical harmonic functions. Defining power to be the integral of the function squared divided by the area it spans, the total power of a function is equal to a sum over its power spectrum \begin{equation} \frac{1}{4\pi} \int_\Omega f^2(\theta,\phi) \, d\Omega = \sum_{l=0}^{\infty} S_{ff}(l), \end{equation} where the power spectrum is related to the spherical harmonic coefficients by
Similarly, the cross power of two functions and is given by \begin{equation} \frac{1}{4\pi} \int_\Omega f(\theta,\phi) \, g(\theta,\phi)\, d\Omega = \sum_{l=0}^{\infty} S_{fg}(l), \end{equation} with
The power spectrum is unmodified by a rotation of the coordinate system. Furthermore, the numerical values of the power spectrum are independent of the normalization convention used for the spherical harmonic functions (though the mathematical formulae will be different, as given below). If the functions and have a zero mean, and represent the contribution to the variance and covariance, respectively, as a function of degree .
is the total power of the function at spherical harmonic degree , which in SHTOOLS is called the power per degree . Alternatively, one can calculate the average power per coefficient at spherical harmonic degree , which in SHTOOLS is referred to as the power per . Since there are spherical harmonic coefficients at degree , this is \begin{equation} \mbox{power per $lm$} = \frac{S(l)}{(2l+1)}. \end{equation} One can also calculate the power from all angular orders over an infinitesimal logarithmic spherical harmonic degree band , where is the logarithmic base. In SHTOOLS, this is referred to as the power per , which is given by \begin{equation} \mbox{power per $d\log_a l$} = S(l)\, l \, \ln a. \end{equation} Finally, SHTOOLS defines the energy of a function as the integral of its square. The energy spectrum is thus equal to the power spectrum multiplied by .
CondonShortley phase factor
The above definitions of the Legendre functions and spherical harmonic functions do not include the CondonShortley phase factor of that is often employed in the physics and seismology communities [Varshalovich et al. 1988, Dahlen and Tromp 1998]. Nevertheless, this phase can be included in most SHTOOLS routines by specifying the optional parameter
csphase = 0
: exclude the CondonShortley phase factor (default)csphase = 1
: append the CondonShortley phase factor to the Legendre functions.
The choice of the CondonShortley phase factor does not affect the numerical value of the power spectrum.
Supported normalizations
SHTOOLS supports the use of normalized, Schmidt seminormalized, orthonormalized, and unnormalized spherical harmonic functions. To specify which normalization should be used, it is only necessary to specify the optional parameter norm
in the Fortran 95 routines or normalization
in the Python routines:
norm = 1
,normalization = '4pi'
: normalized (default, unless stated otherwise)norm = 2
,normalization = 'schmidt'
: Schmidt seminormalizednorm = 3
,normalization = 'unnorm'
: Unnormalizednorm = 4
,normalization = 'ortho'
: Orthonormalized.
Each of these normalizations has slightly different definitions for the normalized Legendre functions, the orthogonality conditions of the Legendre functions and spherical harmonic functions, and the power spectrum. These equations are provided below.
normalized
Schmidt seminormalized
Orthonormalized
Unnormalized
References

Dahlen, F. A. and J. Tromp, “Theoretical Global Seismology,” Princeton University Press, Princeton, New Jersey, 1025 pp., 1998.

Varshalovich, D. A., A. N. Moskalev, and V. K. Khersonskii, “Quantum theory of angular momentum,” World Scientific, Singapore, 1988.