write_fits_cut4


This routine writes a cut sky HEALPix map into 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 write_fits_cut4( filename, np, pixel, signal , n_obs, serror[, header, coord, nside, order, units] )


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 written in the file
pixel(0:np-1) I4B IN index of observed (or valid) pixels
signal(0:np-1) SP IN value of signal in each observed pixel
n_obs(0:np-1) I4B IN number of observation per pixel
serror(0:np-1) SP IN rms of signal in pixel, for white noise, this is $ \propto$ 1/$ \sqrt{{\rm n\_obs}}$.
header(LEN=80)(1:) (OPTIONAL) CHR IN FITS extension header
coord(LEN=1) (OPTIONAL) CHR IN astrophysical coordinates ('C' or 'Q' Celestial/eQuatorial, 'G' for Galactic, 'E' for Ecliptic)
nside (OPTIONAL) I4B IN HEALPix resolution parameter of data set
order (OPTIONAL) I4B IN HEALPix ordering scheme, 1: RING, 2: NESTED
header(LEN=80) (OPTIONAL) CHR IN FITS header to be included in the FITS file
units(LEN=20) (OPTIONAL) CHR IN maps units (applies only to Signal and Serror)
      Note: the information relative to Nside, Order and Coord has to be given, either thru these keyword or via the FITS Header.


MODULES & ROUTINES

This section lists the modules and routines used by write_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 write_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
read_fits_cut4
subroutine to read a HEALPix cut sky map from a FITS file

Eric Hivon 2003-02-07