def sigma_eff_rad(t): z = nfwso.get_z(t) sm = stellar_mass(t) return 190.*(sm/1.e11)**0.2*(1.+z)**0.47
def get_eff_rad(t): z = nfwso.get_z(t) sm = stellar_mass(t) return 2500*(sm/1.e11)**0.73*(1.+z)**-0.98