*** Electron density, X-ray temperature, and electron pressure profiles ***
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]

- To compile and use the code, edit Makefile and simply "make"
- It will generate executables called "compute_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]