Program s2dw_wav2sky

program s2dw_wav2sky

        ! 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_free
    call s2_sky_write_map_file
    call s2dw_core_free_wavdyn
    call s2dw_error
    call s2dw_fileio_fits_wav_read
    call s2dw_fileio_matlab_wav_read

        ! Variables
    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
    character (len=STRING_LEN) :: filename_in
    character (len=STRING_LEN) :: filename_out
    type (s2_sky) :: sky
    integer :: jj
    integer :: gg
    integer :: bl_hi
    integer :: nside
    integer :: pix_scheme
    logical :: interp
    logical :: all
    logical :: use_default_nside
    integer :: S2DW_FITS_FILENAME_EXT_LEN = 4
    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_wav2sky

Converts wavelet coefficients read from a S2DW formatted fits/matlab file (equi-angular pixelisation) to a sky Healpix fits file. Either the sky corresponding to a specified dilation and orientation index may be written only, or skies corresponding to all dilations and orientations may be written.

Usage: s2dw_wav2sky

[-help]: Displays usage information.
[-inp filename_in]: Name of input wavelet coefficient S2DW formatted fits/matlab file.
[-file_type file_type (fits; m)]: String specifying type of input S2DW file to read (fits or matlab m file) [default=fits].
[-out filename_out]: Name of output sky coefficient map (only required if all=.false.)
[-jj jj]: Scale index of coefficients to write (only required if all=.false.)
[-gg gg]: Orientation index of coefficients to write (only required if all=.false.)
[-nside healpix_nside]: Optional Healpix nside resolution of written sky. If this option is not specified then an appropriate nside value is set for each level jj.
[-interp interpolation_status]: Logical to specify whether to perform linear interpolation (if true) or else simply take nearest pixel (default=false).
[-all all_status]: Logical to specify whether skies corresponding to all dilations and orientations are to be written (default=false).

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

Version: 0.1 - November 2007

Description of Variables

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

filename_in

character (len=STRING_LEN) :: filename_in

filename_out

character (len=STRING_LEN) :: filename_out

sky

type (s2_sky) :: sky

jj

integer :: jj

gg

integer :: gg

bl_hi

integer :: bl_hi

nside

integer :: nside

pix_scheme

integer :: pix_scheme

interp

logical :: interp

all

logical :: all

use_default_nside

logical :: use_default_nside

S2DW_FITS_FILENAME_EXT_LEN

integer :: S2DW_FITS_FILENAME_EXT_LEN = 4

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