Program s2_gcmb

program s2_gcmb

        ! Uses
    use s2_types_mod
    use s2_cmb_mod
    use s2_sky_mod
    use s2_pl_mod
    use s2_error_mod
    use s2_vect_mod

        ! Calls
    call parse_options
    call s2_cmb_free
    call s2_error
    call s2_pl_conv
    call s2_pl_free
    call s2_sky_free
    call s2_sky_write_file
    call system_clock

        ! Variables
    type (s2_cmb) :: cmb
    integer :: nside
    integer :: lmin
    integer :: lmax
    integer :: ncomment
    integer :: seed
    integer :: fail = 0
    character (len=S2_STRING_LEN) :: filename_cl
    character (len=S2_STRING_LEN) :: filename_out = 'out.fits'
    logical :: scale_cl = .true.
    logical :: seed_set = .false.
    real (kind=s2_sp) :: noise_var = 0.0e0
    real (kind=s2_sp) :: beam_fwhm = 0.0e0
    real (kind=s2_sp), allocatable, dimension (:) :: npl
    logical :: beam_present = .false.
    logical :: noise_present = .false.
    type (s2_pl) :: beam_pl
    type (s2_pl) :: noise_pl
    type (s2_pl) :: background_pl
    type (s2_pl) :: background_pl_tmp
    type (s2_sky) :: cmb_sky
    character (len=*), parameter :: MAP_FILE = 'map'
    character (len=*), parameter :: ALM_FILE = 'alm'
    character (len=*), parameter :: SKY_FILE = 'sky'
    integer :: file_type_out = S2_SKY_FILE_TYPE_MAP
    character (len=S2_STRING_LEN) :: file_type_out_str = MAP_FILE
    logical :: compute_map = .false.

        ! Subroutines and functions
    subroutine parse_options ()

end program s2_gcmb

Usage:

• [-help]: Display usage information.
• [-inp filename_cl]: Name of cl input ascii file.
• [-scale_cl scale_cl (logical)]: Logical to specify whether to scale the input cl values by 2*pi/(l(l+1)).
• [-out filename_out]: Name of output cmb map file.
• [-file_type_out file_type_out_str]: String specifying output file type.
• [-nside nside]: Healpix nside resolution to generate map for.
• [-lmax lmax]: Maximum hamonic l to consider in power spectrum when generating map.
• [-lmin lmin]: Minimum l of power spectrum contained in input file (others are set to zero). Often this is 2.
• [-ncomment ncomment]: Number of comment lines in input file to ignore.
• [-seed seed]: Integer seed for random number generator used to computed Gaussian realised CMB. If not specified then use system clock.

Author: J. D. McEwen (mcewen@mrao.cam.ac.uk)

Version: 0.1 - July 2005

Description of Variables

type (s2_cmb) :: cmb

integer :: nside

integer :: lmin

integer :: lmax

integer :: ncomment

integer :: seed

integer :: fail = 0

character (len=S2_STRING_LEN) :: filename_cl

character (len=S2_STRING_LEN) :: filename_out = 'out.fits'

logical :: scale_cl = .true.

logical :: seed_set = .false.

real (kind=s2_sp) :: noise_var = 0.0e0

real (kind=s2_sp) :: beam_fwhm = 0.0e0

real (kind=s2_sp), allocatable, dimension (:) :: npl

logical :: beam_present = .false.

logical :: noise_present = .false.

type (s2_pl) :: beam_pl

type (s2_pl) :: noise_pl

type (s2_pl) :: background_pl

type (s2_pl) :: background_pl_tmp

type (s2_sky) :: cmb_sky

character (len=*), parameter :: MAP_FILE = 'map'

character (len=*), parameter :: ALM_FILE = 'alm'

character (len=*), parameter :: SKY_FILE = 'sky'

integer :: file_type_out = S2_SKY_FILE_TYPE_MAP

character (len=S2_STRING_LEN) :: file_type_out_str = MAP_FILE

logical :: compute_map = .false.

Description of Subroutines and Functions

subroutine parse_options ()
    ! Calls: getArgument
end subroutine parse_options

Author: J. D. McEwen (mcewen@mrao.cam.ac.uk)

Version: 0.1 - November 2004