read_conbintab*

This routine reads a FITS file containing $a_{\ell m}$ values for constained realisation. The FITS file is supposed to contain one integer column with $index=\ell^2+\ell+m+1$ and 2 or 4 single (or double) precision columns with real/imaginary $a_{\ell m}$ values and real/imaginary standard deviation on these $a_{\ell m}$. It is supposed to contain either 1 or 3 extension(s) containing respectively the $a_{\ell m}$ for T and if applicable E and B.

Location in HEALPix directory tree: src/f90/mod/fitstools.F90 

FORMAT

call read_conbintab*( filename, alms, nalms[, units, extno] )

ARGUMENTS

name&dimensionality	kind	in/out	description
filename(LEN=filenamelen)	CHR	IN	filename of FITS file containing $a_{\ell m}$.
nalms	I4B	IN	Number of $a_{\ell m}$ values to read from the file.
alms(0:nalms-1,1:6)	SP/ DP	OUT	the $a_{\ell m}$ read from the file. alms(i,1) and alms(i,2) contain the $\ell$ and $m$ values for the ith $a_{\ell m}$ . alms(i,3) and alms(i,4) contain the real and imaginary value of the ith $a_{\ell m}$ . Finally, the standard deviation for the ith $a_{\ell m}$ is contained in alms(i,5) (real) and alms(i,6) (imaginary).
units(len=20)(1:) (OPTIONAL)	CHR	OUT	character string containing the units of the $a_{\ell m}$
extno (OPTIONAL)	I4B	IN	extension (0 based) of the FITS file to be read

EXAMPLE:

call read_conbintab ('alms.fits',alms,65*66/2)

Read 65*66/2 (the number of $a_{\ell m}$ needed to fill the whole range from l=0 to l=64) $a_{\ell m}$ values from the file `alms.fits' into the array alms(0:65*66/2-1,1:6).

MODULES & ROUTINES

This section lists the modules and routines used by read_conbintab*.

fitstools

module, containing:

printerror

routine for printing FITS error messages.

cfitsio

library for FITS file handling.

RELATED ROUTINES

This section lists the routines related to read_conbintab*

alms2fits, dump_alms

routines to write $a_{\ell m}$ to a FITS-file

fits2alms

has the same function as read_conbintab but is more general.

number_of_alms, getsize_fits

can be used to find out the number of $a_{\ell m}$ available in the file.