skymodel.simulate_atmosphere

skymodel.simulate_atmosphere(freq, nside_out=4096, lmax=None, beam_FWHM=None, unit='cmb', white=False, s=1)

Computes an all-sky atmospheric noise map at a given frequency and nside based on the SO noise model presented by the SO Collaboration (2019) and using the model parameters provided by Choi et al. (2019).

Parameters

freq: float or float array

Frequency of the output map in Hz. Must be a valid SO or CCAT-prime central band frequency, i.e. 27, 39, 93, 145, 225, 279, 220, 280, 350, 405, or 860 GHz.

nside_out: float, optional

Healpix nside parameter of the output map. Must be a valid value for nside. Default: 4096

lmax: float, optional

Maximum value of the multipolemoment at which the atmospheric power spectrum wil be computed. Default: 3*nside_out-1

beam_FWHM: bool, optional

If set, the output will be convolved with a gaussian. The FWHM of the Gaussian in units of arcmin is given by the provided value. Default: None

unit: bool, optional

Determines the units of the output map. The available units are ‘mjy’ –> MJy/sr (specific intensity), ‘cmb’ –> K_CMB (thermodynamic temperature), and ‘rj’ –> K_RJ (brightness temperature). Default: ‘cmb’.

white: bool, optional

If True, both the white and the red noise component are simulated. Default: False

s: integer

Is the seed, for each new see a new realization of the noise is made.

Returns

——-

noise_map: float array

Healpix all-sky map of the atmospheric noise at the specified frequency.