aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/halo.py50
-rw-r--r--fg21sim/extragalactic/clusters/main.py4
2 files changed, 0 insertions, 54 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index ed73a47..458c1d0 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -852,56 +852,6 @@ class RadioHaloAM(RadioHalo1M):
rate = (mass0 - mass1) / (t0 - t1)
return (mass1 + rate * (t - t1))
- @property
- def duration_turb_avg(self):
- """
- Calculate the time-averaged turbulence acceleration active time
- within the period from ``age_begin`` to ``age_obs``.
-
- Unit: [Gyr]
- """
- dt = self.time_step
- xt = np.arange(self.age_begin, self.age_obs+dt/2, step=dt)
- t_turb = np.array([self.duration_turb(t) for t in xt])
- return np.sum(t_turb * dt) / (len(xt) * dt)
-
- @property
- def mach_turbulence_avg(self):
- """
- Calculate the time-averaged turbulence Mach number within the
- period from ``age_begin`` to ``age_obs``.
- """
- dt = self.time_step
- xt = np.arange(self.age_begin, self.age_obs+dt/2, step=dt)
- mach = np.array([self.mach_turbulence(t) for t in xt])
- return np.sum(mach * dt) / (len(xt) * dt)
-
- @property
- def tau_acceleration_avg(self):
- """
- Calculate the time-averaged turbulence acceleration timescale
- (i.e., efficiency) within the period from ``age_begin`` to
- ``age_obs``.
-
- Unit: [Gyr]
- """
- dt = self.time_step
- xt = np.arange(self.age_begin, self.age_obs+dt/2, step=dt)
- tau = np.array([self.tau_acceleration(t) for t in xt])
- return np.sum(tau * dt) / (len(xt) * dt)
-
- @property
- def time_acceleration_fraction(self):
- """
- Calculate the fraction of time within the period from
- ``age_begin`` to ``age_obs`` that the turbulence acceleration
- is active.
- """
- dt = self.fpsolver.tstep
- xt = np.arange(self.age_begin, self.age_obs+dt/2, step=dt)
- active = np.array([self._is_turb_active(t) for t in xt], dtype=int)
- return active.mean()
-
class RadioHalo:
"""
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 061b8ef..e169ab9 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -304,10 +304,6 @@ class GalaxyClusters:
n_e=n_e, # [cm^-3]
gamma=halo.gamma, # Lorentz factors
Ke=halo.injection_rate, # [cm^-3 Gyr^-1]
- Mach_turb=halo.mach_turbulence_avg, # Mach number
- tau_turb=halo.duration_turb_avg, # [Gyr]
- tau_acc=halo.tau_acceleration_avg, # [Gyr]
- tfrac_acc=halo.time_acceleration_fraction,
)
return haloinfo