write_plm

This routine creates a double precision binary FITS-file from a given array. The routine is used by the HEALPix facility plmgen to store precomputed $P_{\ell m}(\theta)$.

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


FORMAT

call write_plm( plm, nplm, nhar, header, nlheader, filename, nsmax, nlmax )


ARGUMENTS

name&dimensionality kind in/out description
       
plm(0:nplm-1,1:nhar) DP IN the array with the precomputed $P_{\ell m}(\theta)$ values.
nplm I4B IN Number of $P_{\ell m}$ values to store.
nhar I4B IN 1 for scalar $P_{\ell m}$ only and 3 for tensor harmonics.
header(LEN=80) (1:nlheader) CHR IN The header for the FITS-file.
nlheader I4B IN number of header lines to write to the file.
filename(LEN=filenamelen) CHR IN the precomputed $P_{\ell m}(\theta)$ values are written to this file.
nsmax I4B IN $N_{\mathrm{side}}$ for the precomputed $P_{\ell m}\!\!$ s.
nlmax I4B IN maximum $\ell$ value for the precomputed $P_{\ell m}\!\!$ s.

EXAMPLE:

call write_plm (plm, 65*66*32, 1, header, 120, `plm_32.fits', 32, 64)
Makes a double precision binary FITS-file called `plm_32.fits' from the precomputed $P_{\ell m}(\theta)$ in the array plm(0:65*66*32-1,1:1). The number 65*66*32 corresponds to the number of precomputed $P_{\ell m}\!\!$ s needed for a $N_{\mathrm{side}}=32$ HEALPix map synthesis/analysis. The header for the FITS-file is given in the string array header and the number of lines in the header is 120.


MODULES & ROUTINES

This section lists the modules and routines used by write_plm.

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 write_plm

read_dbintab, read_bintab
routines which reads a file created by write_plm.
map2alm, alm2map
routines using precomputed $P_{\ell m}(\theta)$.

Version 3.82, 2022-07-28