Program s2_sky2plplot

program s2_sky2plplot

        ! Uses
    use s2_types_mod
    use s2_pl_mod
    use s2_sky_mod

        ! Calls
    call parse_options
    call s2_pl_free
    call s2_pl_plot
    call s2_sky_compute_alm
    call s2_sky_dilate
    call s2_sky_free

        ! Variables
    character (len=S2_STRING_LEN) :: filename_sky
    character (len=S2_STRING_LEN) :: filename_out
    logical :: scale = .true.
    character (len=S2_STRING_LEN) :: title = ''
    type (s2_sky) :: sky
    type (s2_pl) :: cl
    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.
    real (kind=s2_sp) :: dil1
    real (kind=s2_sp) :: dil2

        ! Subroutines and functions
    subroutine parse_options ()

end program s2_sky2plplot

Read a full fits s2_sky file and compute and plot the cl spectrum (dilate sky also if required, see notes below).

Notes:

When constructing optimal filters, template sky is saved for dilation of one, when actual template used may be at different scale. Thus functionality incorporated here to dilate sky before computing and plotting cl spectrum.

Usage:

[-help]: Display usage information.
[-inp filename_sky]: Name of full s2_sky fits file containing sky to compute and plot pls of.
[-out filename_out]: Name of output postscript to print.
[-scale scale]: Logical to specify whether to scale y axis values to plot dls rather than cls.
[-title title]: Optional title to print on plot.
[-lmax lmax]: Spherical harmonic lmax if computing alms.
[-mmax mmax]: Spherical harmonic mmax if computing alms.
[-dil1 dil1]: First dilation value if dilating sky before computing alms.
[-dil2 dil2]: Second dilation value if dilating sky before computing alms.
[-npres norm_pres]: Logical to specify whether norm preserving dilation is to be performed, if indeed performing a dilation.

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

Version: 0.1 - April 2005

Description of Variables

filename_sky

character (len=S2_STRING_LEN) :: filename_sky

filename_out

character (len=S2_STRING_LEN) :: filename_out

scale

logical :: scale = .true.

title

character (len=S2_STRING_LEN) :: title = ''

sky

type (s2_sky) :: sky

cl

type (s2_pl) :: cl

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.

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