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.

Version 3.82, 2022-07-28