read_fits_cut4

This routine reads a cut sky HEALPix map from a FITS file. The format used for the FITS file follows the one used for Boomerang98 and is adapted from COBE/DMR

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

FORMAT

call read_fits_cut4( filename, np, pixel, [signal, n_obs, serror, header, units, extno] )

Arguments appearing in italic are optional.

ARGUMENTS

name&dimensionality	kind	in/out	description
filename(LEN=filenamelen)	CHR	IN	FITS file to be read from, containing a cut sky map
np	I4B	IN	number of pixels to be read from the file
pixel(0:np-1)	I4B	OUT	index of observed (or valid) pixels
signal(0:np-1)(OPTIONAL)	SP	OUT	value of signal in each observed pixel
n_obs(0:np-1)	I4B	OUT	number of observation per pixel
serror(0:np-1)	SP	OUT	rms of signal in pixel. (For white noise, this would be $\propto 1/\sqrt{{\rm n\_obs}}$)
header(LEN=80)(1:)	CHR	OUT	FITS extension header
units(LEN=20)	CHR	OUT	maps units (applies only to Signal and Serror, which are assumed to have the same units)
extno	I4B	IN	extension number (0 based) for which map is read. Default = 0 (first extension).

MODULES & ROUTINES

This section lists the modules and routines used by read_fits_cut4.

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_fits_cut4

anafast

executable that reads a HEALPix map and analyses it.

synfast

executable that generate full sky HEALPix maps

getsize_fits

routine to know the size of a FITS file and its type (eg, full sky vs cut sky)

input_map

all purpose routine to input a map of any kind from a FITS file

output_map

subroutine to write a FITS file from a HEALPix map

write_fits_cut4

subroutine to write a cut sky map into a FITS file