aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py18
1 files changed, 12 insertions, 6 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index a18aa17..55107e5 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -241,17 +241,23 @@ class RadioHalo:
The estimated radius of the simulated radio halo.
Unit: [kpc]
"""
- return self.radius_turbulence()
+ return self.radius_turbulence(self.age_merger)
@lru_cache()
- def radius_turbulence(self, t=None):
+ def radius_turbulence(self, t):
"""
- The radius of the turbulence injection region.
+ The radius of the turbulence injection region, which is calculated
+ as the mean of the stripping radius of the sub-cluster and the core
+ radius of the main cluster.
+
Unit: [kpc]
"""
- rs = helper.radius_stripping(self.M_main, self.M_sub, self.z_merger,
- configs=self.configs) # [kpc]
- return self.f_lturb * rs
+ z = COSMO.redshift(t)
+ M_main = self.mass_main(t)
+ M_sub = self.mass_sub(t)
+ rs = helper.radius_stripping(M_main, M_sub, z, configs=self.configs)
+ R_vir = helper.radius_virial(M_main, z)
+ return (R_vir + rs) / 2
@property
def angular_radius(self):