Program cswt_analysis

program cswt_analysis

        ! Uses
    use s2_types_mod
    use s2_sky_mod
    use cswt_error_mod
    use cswt_tmpl_mod
    use cswt_swav_mod
    use cswt_tr_mod

        ! Calls
    call cswt_error
    call cswt_swav_free
    call cswt_tr_analysis
    call cswt_tr_free
    call cswt_tr_io_fits_write_wcoeff
    call parse_options
    call s2_sky_free

        ! Variables
    character (len=S2_STRING_LEN), parameter :: WAV_TYPE_FILE = 'file'
    character (len=S2_STRING_LEN), parameter :: WAV_TYPE_MEXHAT = 'mexhat'
    character (len=S2_STRING_LEN), parameter :: WAV_TYPE_BUTTERFLY = 'butterfly'
    character (len=S2_STRING_LEN), parameter :: WAV_TYPE_MORLET = 'morlet'
    character (len=S2_STRING_LEN) :: filename_dilation
    character (len=S2_STRING_LEN) :: filename_in
    character (len=S2_STRING_LEN) :: filename_out
    character (len=S2_STRING_LEN) :: filename_wavelet
    character (len=S2_STRING_LEN) :: filename_isotropic
    character (len=S2_STRING_LEN) :: wavelet_type
    integer :: i_param = 0
    integer :: n_param = 0
    integer :: n_gamma = 1
    integer :: lmax = 0
    integer :: mmax = 0
    integer :: nside
    integer :: fail
    integer :: file_extension = 1
    real (kind=s2_sp), allocatable, dimension (:) :: param
    integer :: method
    character (len=S2_STRING_LEN), parameter :: METHOD_FFT_REAL = 'fft_real'
    character (len=S2_STRING_LEN), parameter :: METHOD_FFT = 'fft'
    character (len=S2_STRING_LEN), parameter :: METHOD_DFT = 'dft'
    character (len=S2_STRING_LEN), parameter :: METHOD_ISOTROPIC = 'isotropic'
    character (len=S2_STRING_LEN), parameter :: METHOD_DIRECT = 'direct'
    character (len=S2_STRING_LEN) :: method_str = METHOD_FFT_REAL
    logical :: message = .true.
    type (s2_sky) :: data
    type (cswt_swav) :: swav_mother
    type (cswt_tr) :: tr
    logical :: time_it = .false.
    logical :: isotropic_out = .false.

        ! Subroutines and functions
    subroutine parse_options ()

end program cswt_analysis

Usage: cswt_analysis

• [-help]: Display usage information.
• [-dil filename_dilation]: Name of dilation file containing dilations to consider.
• [-inp filename_data]: Name of data file containing sky to take transform of.
• [-ext file_extension]: Extension of filename_data to read sky from.
• [-out filename_out]: Name of output wavelet coefficient file.
• [-method method_str]: String specifying cswt algorithm to use.
• [-message message]: Logical specifying whether to display messages.
• [-filename_isotropic filename_isotropic]: Optional name of fits output file of wavelet coefficients in healpix format for isotropic algorithm.
• [-wav wavelet_type]: String specifying wavelet type.
• [-wavfile filename_wavelet]: Name of wavelet file if wavelet is to be read from file (i.e. -wav file)
• [-lmax lmax]: lmax of wavelet transform.
• [-mmax mmax]: mmax of wavelet transform.
• [-ngamma n_gamma]: Number of orientations to consider in transform.
• [-time time_it]: Logical specifying whether to time algorithm.
• [-npar n_param]: Number of wavelet parameters specified.
• [-par param(i) (repeat n_param times)]: Wavelet parameters. Parameters are read sequentially on order appear in command line.

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

Version: 0.1 - November 2004

Description of Variables

character (len=S2_STRING_LEN), parameter :: WAV_TYPE_FILE = 'file'

character (len=S2_STRING_LEN), parameter :: WAV_TYPE_MEXHAT = 'mexhat'

character (len=S2_STRING_LEN), parameter :: WAV_TYPE_BUTTERFLY = 'butterfly'

character (len=S2_STRING_LEN), parameter :: WAV_TYPE_MORLET = 'morlet'

character (len=S2_STRING_LEN) :: filename_dilation

character (len=S2_STRING_LEN) :: filename_in

character (len=S2_STRING_LEN) :: filename_out

character (len=S2_STRING_LEN) :: filename_wavelet

character (len=S2_STRING_LEN) :: filename_isotropic

character (len=S2_STRING_LEN) :: wavelet_type

integer :: i_param = 0

integer :: n_param = 0

integer :: n_gamma = 1

integer :: lmax = 0

integer :: mmax = 0

integer :: nside

integer :: fail

integer :: file_extension = 1

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

integer :: method

character (len=S2_STRING_LEN), parameter :: METHOD_FFT_REAL = 'fft_real'

character (len=S2_STRING_LEN), parameter :: METHOD_FFT = 'fft'

character (len=S2_STRING_LEN), parameter :: METHOD_DFT = 'dft'

character (len=S2_STRING_LEN), parameter :: METHOD_ISOTROPIC = 'isotropic'

character (len=S2_STRING_LEN), parameter :: METHOD_DIRECT = 'direct'

character (len=S2_STRING_LEN) :: method_str = METHOD_FFT_REAL

logical :: message = .true.

type (s2_sky) :: data

type (cswt_swav) :: swav_mother

type (cswt_tr) :: tr

logical :: time_it = .false.

logical :: isotropic_out = .false.

Description of Subroutines and Functions

subroutine parse_options ()
    ! Calls: cswt_error, getArgument
end subroutine parse_options

Version: 0.1 - November 2004