The HEALPix package contains the Fortran90 facility synfast
which takes as input a power spectrum Cl and generates a realisation of
on the HEALPix grid. The convention for power spectrum input into
synfast is straightforward: each Cl is just the expected
variance of the alm at that l.
Example: The spherical harmonic coefficient a00 is the integral of theNote that this definition implies the standard result that the total power at the angular wavenumber l is (2l+1)Cl, because there are 2l+1 modes for each l.over the sphere. To obtain realisations of functions which have a00 distributed as a Gaussian with zero mean and variance 1, set C0 to 1. The value of the synthesised function at each pixel will be Gaussian distributed with mean zero and variance
. As required, the integral of
over the full
solid angle of the sphere has zero mean and variance
.
This defines unambiguously how the Cl have to be defined given the units of the physical quantity f. In cosmic microwave background research, popular choices for simulated maps are
Eric Hivon 2005-08-31