generate_beam
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 B_{l} 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 B_{l} 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