*** Sunyaev-Zel'dovich profile (in Rayleigh-Jeans temperature) of Coma cluster ***
November 20, 2011: E. Komatsu

This code uses the generalized-NFW fit of Alexis Finoguenov (in prep) to compute the electron pressure profile of the Coma cluster,

Pe = Pe0*(r/r0)^(-gamma)*[1+(r/r0)^alpha]^{-(beta-gamma)/alpha} [eV cm^-3]

and then integrates the pressure along the line of sight to obtain the profile of the Sunyaev-Zel'dovich effect in units of the Rayleigh-Jeans temperature:

Tsz(theta) = -2*273[uK] * int dl Pe/(25eV cm^-3 Mpc)

The integration boundary is from -sqrt(rout^2-theta^2 Da^2) to +sqrt(rout^2-theta^2 Da^2).
[See equations (74) and (75) of E.Komatsu et al., Astrophysical Journal Supplement Series, 192, 18 (2011).]

- To compile and use the code, edit Makefile and simply "make"
- It will generate executables called "compute_coma_profiles"
- The output includes:
  1st column: x=theta/theta500, where theta500=r500/(angular diameter distance)
  2nd column: Sunyaev-Zel'dovich effect in RJ units, Tsz [uK]