Actually return rho_b instead of log10(rho_b)