Arguments appearing in italic are optional.
call generate_beam( fwhm_arcmin, lmax, beam[, beam_file] )
|name & dimensionality||kind||in/out||description|
|fwhm_arcmin||DP||IN||fwhm size of the gaussian beam in arcminutes.|
|lmax||I4B||IN||maximum l value of the window function.|
|beam(0:lmax,1:p)||DP||OUT||beam window function generated. The second index runs form 1:1 for temperature only, and 1:3 for polarisation. In the latter case, 1=T, 2=E, 3=B.|
|beam_file(LEN=filenamelen) (OPTIONAL)||CHR||IN||name of the file containing the (non necessarily gaussian) window function Bl of a circular beam. If present, it will override the argument fwhm_arcmin. If fewer columns than requested are found in the file, missing colums will duplicate the existing ones (based on the assumption that Bl is the same in T, E and B). Supports the fitsio 'Extended Filename Syntax' (see examples below).|
gb0 will contain the window function of a gaussian beam of FWHM = 5 arcmin, for .
b1 will contain the first 3 columns (if available) of file.fits. If the file contains only two columns, then b1(:,3) = b1(:,2), and if it contains a single column, then b1(:,3) = b1(:,2) = b1(:,1).
b2 will be based on a virtual FITS file containg only the first column of file.fits, and we will have b2(:,3) = b2(:,2) = b2(:,1).
Finally b3 will read a virtual FITS file in which the first column is the same as in file.fits, while the columns 2 and 3 are set to 0. Therefore b3(:,3) = b3(:,2) = 0.
Version 3.31, 2017-01-06