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

This code uses equation (3) and (8) of Vikhlinin et al., Astrophysical Journal, 640, 691 (2006)
to compute the electron density, X-ray temperature, and electron pressure profiles:

ne^2 = n0**2d0*(r/rc)**(-alpha)/(1d0+r**2d0/rc**2d0)**(3d0*beta-alpha/2d0) &
      /(1d0+r**gamma/rs**gamma)**(eps/gamma) + n02**2d0/(1d0+r**2d0/rc2**2d0)**(3d0*beta2) ! cm^-6

Te = 1.35d0*(t3/1.11d0/0.97d0)*((x/0.045d0)**1.9d0+0.45d0)/((x/0.045d0)**1.9d0+1d0) &
    /(1d0+(x/0.6d0)**2d0)**0.45d0 ! keV

Pe = ne*(1000*Te) [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_szprofiles"
- 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]