Program s2dw_synthesis

program s2dw_synthesis

        ! Uses
    use s2dw_types_mod
    use s2dw_error_mod
    use s2dw_fileio_mod
    use s2dw_core_mod
    use s2_sky_mod

        ! Calls
    call parse_options
    call s2_sky_compute_map
    call s2_sky_free
    call s2_sky_write_map_file
    call s2dw_core_free_wavdyn
    call s2dw_core_init_directionality
    call s2dw_core_init_kernels
    call s2dw_core_synthesis_wav2flm_dynamic
    call s2dw_error
    call s2dw_fileio_fits_wav_read
    call s2dw_fileio_matlab_wav_read

        ! Variables
    complex (kind=dpc), allocatable, dimension (:,:) :: flm
    complex (kind=spc), allocatable, dimension (:,:) :: flm_temp
    real (kind=dp), allocatable, dimension (:,:) :: K_gamma
    real (kind=dp), allocatable, dimension (:) :: Phi2
    complex (kind=dpc), allocatable, dimension (:,:) :: Slm
    real (kind=dp), allocatable, dimension (:) :: admiss
    type (s2dw_wav_abg), allocatable, dimension (:) :: wavdyn
    complex (kind=dpc), allocatable, dimension (:,:) :: scoeff
    integer :: J
    integer :: B
    integer :: N
    integer :: bl_scoeff
    real (kind=dp) :: alpha
    integer :: nside
    logical :: nside_in = .false.
    logical :: admiss_pass
    type (s2_sky) :: sky
    integer :: fail = 0
    character (len=STRING_LEN) :: filename_in
    character (len=STRING_LEN) :: filename_out
    character (len=*), parameter :: FILE_TYPE_FITS = 'fits'
    character (len=*), parameter :: FILE_TYPE_MAT = 'm'
    character (len=STRING_LEN) :: file_type = FILE_TYPE_FITS

        ! Subroutines and functions
    subroutine parse_options ()

end program s2dw_synthesis

Reconstructs a Healpix sky map from S2DW wavelet and scaling coefficients.

Notes:

The Healpix nside resolution parameter of the output sky map is set by nside=B/2.
The spherical harmonic transform and inverse provided by Healpix is not exact, hence these reduce the numerical precision of the `exact' reconstruction of the original real space map from its wavelet coefficients.

Usage: s2dw_synthesis

[-help]: Displays usage information.
[-inp filename_in]: Name of input S2DW formatted fits/matlab file containing wavelet and scaling coefficients.
[-out filename_out]: Name of output Healpix sky fits map.
[-file_type file_type (fits; m)]: String specifying type of input S2DW file to read (fits or matlab m file) [default=fits].
[-nside nside]: Resolution of output Healpix map (if not specified then set to B/2).

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

Version: 0.1 - November 2007

Description of Variables

flm

complex (kind=dpc), allocatable, dimension (:,:) :: flm

flm_temp

complex (kind=spc), allocatable, dimension (:,:) :: flm_temp

K_gamma

real (kind=dp), allocatable, dimension (:,:) :: K_gamma

Phi2

real (kind=dp), allocatable, dimension (:) :: Phi2

Slm

complex (kind=dpc), allocatable, dimension (:,:) :: Slm

admiss

real (kind=dp), allocatable, dimension (:) :: admiss

wavdyn

type (s2dw_wav_abg), allocatable, dimension (:) :: wavdyn

scoeff

complex (kind=dpc), allocatable, dimension (:,:) :: scoeff

J

integer :: J

B

integer :: B

N

integer :: N

bl_scoeff

integer :: bl_scoeff

alpha

real (kind=dp) :: alpha

nside

integer :: nside

nside_in

logical :: nside_in = .false.

admiss_pass

logical :: admiss_pass

sky

type (s2_sky) :: sky

fail

integer :: fail = 0

filename_in

character (len=STRING_LEN) :: filename_in

filename_out

character (len=STRING_LEN) :: filename_out

FILE_TYPE_FITS

character (len=*), parameter :: FILE_TYPE_FITS = 'fits'

FILE_TYPE_MAT

character (len=*), parameter :: FILE_TYPE_MAT = 'm'

file_type

character (len=STRING_LEN) :: file_type = FILE_TYPE_FITS

Description of Subroutines and Functions

parse_options

subroutine parse_options ()
    ! Calls: getArgument
end subroutine parse_options

Parses the options passed when program called. ! @author J. D. McEwen (mcewen@mrao.cam.ac.uk)

Version: 0.1 - November 2007