clusters/halo: Add option "kT_out" to tune cluster temperature

2 files changed, 7 insertions, 5 deletions

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index d3b5dd1..bf8e643 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -327,9 +327,7 @@ class RadioHalo:
         if not hasattr(self, "_tau_acceleration"):
             Mach = self.Mach_turbulence
             Rvir = helper.radius_virial(mass=self.M_main, z=self.z_merger)
-            kT = helper.kT_cluster(mass=self.M_main, z=self.z_merger,
-                                   radius=Rvir)  # [keV]
-            cs = helper.speed_sound(kT)  # [km/s]
+            cs = helper.speed_sound(self.kT_main)  # [km/s]
             # Turbulence injection scale
             L0 = self.f_lturb * Rvir  # [kpc]
             x = cs*AUC.km2cm / AC.c
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index d77c63e..86662e1 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -153,7 +153,9 @@ def kT_cluster(mass, z=0.0, radius=None):
     the cluster (near the virial radius) which can heat the gas,
     therefore the ICM has a higher temperature than the virial
     temperature, which can be estimated as:
-        T_icm = T_vir + 1.5*T_1 ~ T_vir + 0.8 [keV]
+        kT_icm ~ kT_vir + 1.5 * kT_out
+    where kT_out the temperature of the outer gas surround the cluster,
+    which may be ~0.5-1.0 keV.
 
     Reference: Ref.[fujita2003],Eq.(49)
 
@@ -163,8 +165,10 @@ def kT_cluster(mass, z=0.0, radius=None):
        The temperature of the cluster ICM.
        Unit: [keV]
     """
+    key = "extragalactic/clusters/kT_out"
+    kT_out = CONFIGS.getn(key)
     kT_vir = kT_virial(mass=mass, z=z, radius=radius)
-    kT_icm = kT_vir + 0.8
+    kT_icm = kT_vir + 1.5*kT_out
     return kT_icm