aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py7
1 files changed, 4 insertions, 3 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 0ffca1b..6ecdffa 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -154,6 +154,7 @@ class RadioHalo1M:
self.eta_turb = configs.getn(sec+"/eta_turb")
self.eta_e = configs.getn(sec+"/eta_e")
self.x_cr = configs.getn(sec+"/x_cr")
+ self.eta_b = self.x_cr # Equipartition between magnetic field and CR
self.mass_index = configs.getn(sec+"/mass_index")
self.gamma_min = configs.getn(sec+"/gamma_min")
self.gamma_max = configs.getn(sec+"/gamma_max")
@@ -170,6 +171,8 @@ class RadioHalo1M:
self.bolo_freq_max = configs.getn(sec+"/bolo_freq_max")
self.bolo_freq_num = configs.getn(sec+"/bolo_freq_num")
+ self.kT_out = configs.getn("extragalactic/clusters/kT_out")
+
def _set_solver(self):
self.fpsolver = FokkerPlanckSolver(
xmin=self.gamma_min,
@@ -681,12 +684,10 @@ class RadioHalo1M:
Unit: [uG]
"""
- eta_b = self.x_cr # Equipartition between magnetic field and CR
- kT_out = self.configs.getn("extragalactic/clusters/kT_out")
z = COSMO.redshift(t)
mass = self.mass_main(t) # [Msun]
return helper.magnetic_field(mass=mass, z=z,
- eta_b=eta_b, kT_out=kT_out)
+ eta_b=self.eta_b, kT_out=self.kT_out)
def _is_turb_active(self, t):
"""