aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2018-01-01 13:57:22 +0800
committerAaron LI <aly@aaronly.me>2018-01-01 13:57:22 +0800
commit21665d0966d456613e7846a74cc381b78943fd0e (patch)
treea2f49881a21ec12147d03c188ff647bf0b5a85d4
parent6a458c2c4676f622ac3151526bd0642f1a191100 (diff)
downloadfg21sim-21665d0966d456613e7846a74cc381b78943fd0e.tar.bz2
clusters/halo: adjust radio halo size estimation
-rw-r--r--fg21sim/extragalactic/clusters/helper.py6
1 files changed, 3 insertions, 3 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index f5b4f8c..29134f7 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -1,4 +1,4 @@
-# Copyright (c) 2017 Weitian LI <weitian@aaronly.me>
+# Copyright (c) 2017-2018 Weitian LI <weitian@aaronly.me>
# MIT license
"""
@@ -98,7 +98,7 @@ def radius_halo(M_main, M_sub, z=0.0, configs=CONFIGS):
The halo radius is estimated to be the same as the turbulence
injection scale, i.e.:
R_halo ≅ L ≅ R_vir / 3
- where R_vir the virial radius of the main cluster.
+ where R_vir the virial radius of the merged cluster.
Reference: [vazza2011],Sec.(3.6)
@@ -120,7 +120,7 @@ def radius_halo(M_main, M_sub, z=0.0, configs=CONFIGS):
# Turbulence injection scale factor
key = "extragalactic/halos/f_lturb"
f_lturb = configs.getn(key)
- R_halo = f_lturb * radius_virial(mass=M_main, z=z) # [kpc]
+ R_halo = f_lturb * radius_virial(mass=M_main+M_sub, z=z) # [kpc]
return R_halo