clusters/halo: Implement radius_turb_eff() method

1 files changed, 27 insertions, 0 deletions

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 94cb1c0..c5da9eb 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -778,6 +778,33 @@ class RadioHaloAM(RadioHalo1M):
         """
         return max([self.radius_turb(tm) for tm in self.t_merger])
 
+    def radius_turb_eff(self, t, use_last=True):
+        """
+        Get the effective turbulence radius, i.e., the largest one if
+        multiple mergers are active at the given time.
+
+        Parameters
+        ----------
+        use_last : bool
+            If ``True``, return the turbulence radius of the last merger
+            event when there is no active turbulence at the given time.
+            Otherwise, return 0.
+
+        Unit: [kpc]
+        """
+        mergers = [(t, t+self.duration_turb(t), self.radius_turb(t))
+                   for t in self.t_merger]  # time decreasing
+        try:
+            r_eff = max([r for t1, t2, r in mergers if t >= t1 and t < t2])
+        except ValueError:
+            # No active turbulence at this time
+            if use_last:
+                r_eff = next(r for __, t2, r in mergers if t >= t2)
+            else:
+                r_eff = 0
+
+        return r_eff
+
     @property
     def t_begin(self):
         """