def get_rh(Mtot):
H = calculate_r200.get_H(z)
rs = calculate_r200.get_r200(Mtot, H)
concentration = 10.**(1.025-0.097*np.log(Mtot/((10.**12)/little_h)))
rh = rs*(0.6082 - 0.1843*np.log(concentration) - 0.1011*(np.log(concentration)**2.) + \
0.03918*(np.log(concentration)**3.))
return rh
def get_NFW_parameters(t, H, gal):
m200 = galaxy_mass(mg[gal], t)
r200 = get_r200(m200, H)
c = get_c(m200)
R_S = get_RS(r200, c)
rho_crit = get_rho_crit(H, m200, r200)
rho_0 = get_rho_0(rho_crit, c)
return (m200, r200, c, R_S, rho_crit, rho_0)