diff options
| author | Aaron LI <aly@aaronly.me> | 2019-02-16 11:23:23 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2019-02-16 11:23:23 +0800 |
| commit | 0da3d12a40ead2c2553e425043cfc903f19de0fb (patch) | |
| tree | 091836c28b480cace9c0e33756ec01cd5e98071e /fg21sim/extragalactic | |
| parent | 04a5cae1c09b46e485f30f148e60a556ebb11748 (diff) | |
| download | fg21sim-0da3d12a40ead2c2553e425043cfc903f19de0fb.tar.bz2 | |
clusters/halo: Clean up the magnetic field calculation a bit
Diffstat (limited to 'fg21sim/extragalactic')
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 7 |
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): """ |