Program s2_sky2cl

program s2_sky2cl

        ! Uses
    use s2_types_mod
    use s2_error_mod
    use s2_sky_mod
    use s2_pl_mod

        ! Calls
    call parse_options
    call s2_error
    call s2_pl_conv
    call s2_pl_free
    call s2_pl_io_fits_write
    call s2_sky_compute_alm
    call s2_sky_dilate
    call s2_sky_free

        ! Variables
    character (len=*), parameter :: MAP_FILE = 'map'
    character (len=*), parameter :: SKY_FILE = 'sky'
    character (len=S2_STRING_LEN) :: filename_in
    character (len=S2_STRING_LEN) :: filename_out
    character (len=S2_STRING_LEN) :: filename_beam
    character (len=S2_STRING_LEN) :: file_type = SKY_FILE
    type (s2_sky) :: sky
    type (s2_pl) :: cl
    type (s2_pl) :: beam
    integer :: lmax = 0
    integer :: mmax = 0
    logical :: new_lmax = .false.
    logical :: new_mmax = .false.
    logical :: dil1_status = .false.
    logical :: dil2_status = .false.
    logical :: norm_pres = .false.
    logical :: beam_status = .false.
    real (kind=s2_sp) :: dil1
    real (kind=s2_sp) :: dil2

        ! Subroutines and functions
    subroutine parse_options ()

end program s2_sky2cl

Read a map from a fits file, compute the alms of the map and then the cl spectrum. Write the computed cl spectrum values to the standard output.

Usage:

[-help]: Display usage information.
[-inp filename_in]: Name of sky fits file containing map to compute cls of.
[-out filename_out]: Name of ouput pl file.
[-beam filename_beam]: Name of beam file to convolve with only if present.
[-type file_type (''map'' or ''sky'')]: String to specify type of input file.
[-lmax lmax]: lmax of alm representation of map.
[-mmax mmax]: mmax of alm representation of map.

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

Version: 0.1 - March 2005

Description of Variables

MAP_FILE

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

SKY_FILE

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

filename_in

character (len=S2_STRING_LEN) :: filename_in

filename_out

character (len=S2_STRING_LEN) :: filename_out

filename_beam

character (len=S2_STRING_LEN) :: filename_beam

file_type

character (len=S2_STRING_LEN) :: file_type = SKY_FILE

sky

type (s2_sky) :: sky

cl

type (s2_pl) :: cl

beam

type (s2_pl) :: beam

lmax

integer :: lmax = 0

mmax

integer :: mmax = 0

new_lmax

logical :: new_lmax = .false.

new_mmax

logical :: new_mmax = .false.

dil1_status

logical :: dil1_status = .false.

dil2_status

logical :: dil2_status = .false.

norm_pres

logical :: norm_pres = .false.

beam_status

logical :: beam_status = .false.

dil1

real (kind=s2_sp) :: dil1

dil2

real (kind=s2_sp) :: dil2

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