WMAP Software for First Year Data Release:  Likelihood Code
*****************************************************************************
If you use this code in any publication, please reference Verde et al. (2003),
if you use the accompanying data files, please also reference 
Hinshaw et al. (2003) and Kogut et al. (2003). 
*****************************************************************************


Description: 
------------

FORTRAN subroutines are provided to compute the likelihood of a user-supplied
theoretical power spectrum against the measured WMAP TT and/or TE power 
spectrum.  These subroutines should compile with a FORTRAN90 compiler.

A sample FORTRAN90 wrapper program in included to illustrate the usage and
calling sequence for subroutines.


Contents:
---------

test_driver.f90           A sample wrapper program.  This program compares
                          a test power spectrum to the data.

compute_likelihood.F      File that contains the collected subroutines
                          that perform the likelihood computation.

tt_diag.dat,              ASCII data files containing diagonal and off-diagonal
te_diag.dat,              terms for Fisher matrix computation.
tt_offdiag.dat,
te_offdiag.dat,

testmodel.dat             ASCII data for a test model used to verify the code.	  


Building the code:
------------------

Commands may vary depending on platform.  A sample command for 
a unix system is:
> f90 -o test_driver test_driver.f90 compute_likelihood.F


User Input:
-----------

The user must provide a theoretical power spectrum, l*(l+1)/2pi * C_l, 
in units of uK^2, starting at l=2. 


Output:
-------
The log likelihood (ln L), or the Fisher matrix can be returned, 
depending on the subroutine called:

compute_mapt_likelihood  -- Returns ln L for TT power spectrum.
compute_mapt_fisher      -- Returns  TT Fisher matrix.
compute_mapte_likelihood -- Returns ln L for TE power spectrum.
compute_mapte_fisher     -- Returns TE Fisher matrix.

The chi-squared of the fit is roughly -2*(ln L).  There are 899 data points
in the WMAP TT power spectrum, and 449 in the TE power spectrum.