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

Version 3.31, 2017-01-06