SHTOOLS --- Tools for working with spherical harmonics

> Home > Documentation > Fortran 95 > Global Spectral Analysis

SHPowerSpectrumDensity

Compute the power spectral density of a real function.

Usage

call SHPowerSpectrumDensity (cilm, lmax, pspectrum)

Parameters

cilm : input, real*8, dimension (2, lmaxin+1, lmaxin+1)
The real 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 spectral density of the function.

Description

SHPowerSpectrumDensity will calculate the power spectral density of a function expressed in real spherical harmonics. For a given spherical harmonic degree l, this is explicitly calculated as:

pspectrum(l) = Sum_{i=1}^2 Sum_{m=0}^l cilm(i, l+1, m+1)**2 / (2l + 1).

See also

shpowerl, shpowerdensityl, shcrosspowerl, shcrosspowerdensityl, shpowerspectrum, shcrosspowerspectrum, shcrosspowerspectrumdensity, shadmitcorr

> Home > Documentation > Fortran 95 > Global Spectral Analysis

Institut de Physique du Globe de Paris University of Sorbonne Paris Cité © 2016 SHTOOLS