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
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