Compute the power spectrum of a real function.
call SHPowerSpectrum (cilm
, lmax
, pspectrum
)
cilm
: input, real*8, dimension (2, lmaxin
+1, lmaxin
+1)
: The function expressed in real spherical harmonics.
lmax
: input, integer
: The maximum spherical harmonic degree used in calculating the power spectrum. This must be less than or equal to lmaxin
.
pspectrum
: output, real*8, dimension (lmax
+1)
: The power spectrum of the function.
SHPowerSpectrum
will calculate the power spectrum of a function expressed in real spherical harmonics. For a given spherical harmonic degree l
, this is calculated explicitly as:
pspectrum(l) = Sum_{i=1}^2 Sum_{m=0}^l cilm(i, l+1, m+1)**2
.
shpowerl, shpowerdensityl, shcrosspowerl, shcrosspowerdensityl, shpowerspectrumdensity, shcrosspowerspectrum, shcrosspowerspectrumdensity, shadmitcorr