aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rw-r--r--fg21sim/extragalactic/clusters/helper.py40
-rw-r--r--fg21sim/extragalactic/clusters/main.py24
2 files changed, 40 insertions, 24 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index f07f1c0..f20af5f 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -17,11 +17,16 @@ References
.. [cassano2007]
Cassano et al. 2007, MNRAS, 378, 1565;
http://adsabs.harvard.edu/abs/2007MNRAS.378.1565C
+
+.. [cassano2012]
+ Cassano et al. 2012, A&A, 548, A100
+ http://adsabs.harvard.edu/abs/2012A%26A...548A.100C
"""
import numpy as np
+from ...configs import configs
from ...utils import cosmo
from ...utils.units import (Constants as AC,
UnitConversions as AUC)
@@ -127,6 +132,11 @@ def density_number_thermal(mass, z=0.0):
"""
Calculate the number density of the ICM thermal plasma.
+ NOTE
+ ----
+ This number density is independent of cluster (virial) mass,
+ but (mostly) increases with redshifts.
+
Parameters
----------
mass : float
@@ -224,3 +234,33 @@ def time_crossing(M_main, M_sub, z=0.0):
uconv = AUC.kpc2km * AUC.s2Gyr
time = uconv * R_vir / vi # [Gyr]
return time
+
+
+def magnetic_field(mass):
+ """
+ Calculate the mean magnetic field strength according to the
+ scaling relation between magnetic field and cluster mass.
+
+ Parameters
+ ----------
+ mass : float
+ Cluster mass
+ Unit: [Msun]
+
+ Returns
+ -------
+ B : float
+ The mean magnetic field strength
+ Unit: [uG]
+
+ References
+ ----------
+ Ref.[cassano2012],Eq.(1)
+ """
+ comp = "extragalactic/clusters"
+ b_mean = configs.getn(comp+"/b_mean")
+ b_index = configs.getn(comp+"/b_index")
+
+ M_mean = 1.6e15 # [Msun]
+ B = b_mean * (mass/M_mean) ** b_index
+ return B
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 840acbe..db82e78 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -75,8 +75,6 @@ class GalaxyClusters:
self.merger_mass_min = self.configs.getn(comp+"/merger_mass_min")
self.ratio_major = self.configs.getn(comp+"/ratio_major")
self.tau_merger = self.configs.getn(comp+"/tau_merger")
- self.b_mean = self.configs.getn(comp+"/b_mean")
- self.b_index = self.configs.getn(comp+"/b_index")
self.filename_pattern = self.configs.getn("output/filename_pattern")
self.use_float = self.configs.getn("output/use_float")
@@ -241,30 +239,8 @@ class GalaxyClusters:
logger.info("%d (%.1f%%) clusters have recent major mergers." %
(num_major, 100*num_major/num))
- def _magnetic_field(self, mass):
"""
- Calculate the mean magnetic field strength according to the
- scaling relation between magnetic field and cluster mass.
-
- Parameters
- ----------
- mass : float
- Cluster mass
- Unit: [Msun]
-
- Returns
- -------
- B : float
- The mean magnetic field strength
- Unit: [uG]
-
- References
- ----------
- Ref.[cassano2012],Eq.(1)
"""
- M_mean = 1.6e15 # [Msun]
- B = self.b_mean * (mass/M_mean) ** self.b_index
- return B
def preprocess(self):
"""