Program s2_nmask

program s2_nmask

        ! Uses
    use s2_types_mod
    use s2_error_mod
    use s2_sky_mod
    use s2_distn_mod
    use s2_ylm_mod

        ! Calls
    call parse_options
    call s2_error
    call s2_sky_free
    call s2_sky_get_map
    call s2_sky_io_fits_write
    call s2_sky_write_map_file

        ! Variables
    real (kind=s2_sp), parameter :: ZERO_TOL = 1e-5
    character (len=*), parameter :: MAP_FILE = 'map'
    character (len=*), parameter :: SKY_FILE = 'sky'
    integer :: file_type = S2_SKY_FILE_TYPE_MAP
    integer :: ext = 1
    character (len=S2_STRING_LEN) :: file_type_str = MAP_FILE
    character (len=S2_STRING_LEN) :: filename_mask
    character (len=S2_STRING_LEN) :: filename_out
    type (s2_sky) :: mask
    type (s2_sky) :: nmask
    real (kind=s2_sp), allocatable, dimension (:) :: mask_map
    real (kind=s2_sp), allocatable, dimension (:) :: nmask_map
    type (s2_sky) :: ylm1_real
    type (s2_sky) :: ylm1_imag
    type (s2_sky) :: ylm2_real
    type (s2_sky) :: ylm2_imag
    real (kind=s2_sp), allocatable, dimension (:) :: ylm1_real_map
    real (kind=s2_sp), allocatable, dimension (:) :: ylm1_imag_map
    real (kind=s2_sp), allocatable, dimension (:) :: ylm2_real_map
    real (kind=s2_sp), allocatable, dimension (:) :: ylm2_imag_map
    complex (kind=s2_spc), allocatable, dimension (:,:,:,:) :: cov
    integer :: ipix
    integer :: npix
    integer :: nside
    integer :: pix_scheme
    integer :: lmax
    integer :: fail
    integer :: seed
    integer :: el1
    integer :: m1
    integer :: el2
    integer :: m2
    real (kind=s2_sp) :: nsigma = 1e0
    complex (kind=s2_dpc), parameter :: I = cmplx (0d0, 1d0)

        ! Subroutines and functions
    subroutine parse_options ()

end program s2_nmask

Generate noisy masks and computed harmonic space noise covariance matrix (not yet outputted) from a binary mask.

Usage: s2_nmask

[-help]: Display usage information.
[-mask filename_mask]: Name of input file containined mask construct noisy mask from.
[-lmax lmax]: Harmonic lmax for computing covariance matrix.
[-nsigma nsigma]: Noise standard deviation.
[-out filename_out]: Name of output noisy mask.
[-file_type file_type_str (map or sky)]: String specifying file types.
[-ext ext (optional)]: File extension for map files.

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

Version: 0.1 - February 2008

Description of Variables

ZERO_TOL

real (kind=s2_sp), parameter :: ZERO_TOL = 1e-5

MAP_FILE

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

SKY_FILE

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

file_type

integer :: file_type = S2_SKY_FILE_TYPE_MAP

ext

integer :: ext = 1

file_type_str

character (len=S2_STRING_LEN) :: file_type_str = MAP_FILE

filename_mask

character (len=S2_STRING_LEN) :: filename_mask

filename_out

character (len=S2_STRING_LEN) :: filename_out

mask

type (s2_sky) :: mask

nmask

type (s2_sky) :: nmask

mask_map

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

nmask_map

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

ylm1_real

type (s2_sky) :: ylm1_real

ylm1_imag

type (s2_sky) :: ylm1_imag

ylm2_real

type (s2_sky) :: ylm2_real

ylm2_imag

type (s2_sky) :: ylm2_imag

ylm1_real_map

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

ylm1_imag_map

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

ylm2_real_map

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

ylm2_imag_map

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

cov

complex (kind=s2_spc), allocatable, dimension (:,:,:,:) :: cov

ipix

integer :: ipix

npix

integer :: npix

nside

integer :: nside

pix_scheme

integer :: pix_scheme

lmax

integer :: lmax

fail

integer :: fail

seed

integer :: seed

el1

integer :: el1

m1

integer :: m1

el2

integer :: el2

m2

integer :: m2

nsigma

real (kind=s2_sp) :: nsigma = 1e0

I

complex (kind=s2_dpc), parameter :: I = cmplx (0d0, 1d0)

Description of Subroutines and Functions

parse_options

subroutine parse_options ()
    ! Calls: getArgument
end subroutine parse_options

Parse the options passed when program called.

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

Version: 0.1 - November 2004