*** Electron density, X-ray temperature, and electron pressure profiles of Coma cluster ***
November 20, 2011: E. Komatsu

This code uses the beta-model fit of Briel et al. (1992) and the generalized-NFW fit of Alexis Finoguenov (in prep) to compute the electron density, X-ray temperature, and electron pressure profiles:

ne = ne0*(1+(r/rc)^2)^(-3*beta/2) [cm^-3]

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

Te = Pe/ne/1e3 [keV] 

- 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=r/r500
  2nd column: r [h^-1 Mpc]
  3rd column: electron number density profile, ne [cm^-3]
  4th column: X-ray temperature profile, Te [keV]
  5th column: electron pressure profile, Pe [eV cm^-3]