From 157028160de3a7baa640cbe87bf5f1a3c1b90ad7 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Thu, 21 Feb 2019 20:52:14 +0800 Subject: clusters/halo: Several minor improvements --- fg21sim/extragalactic/clusters/halo.py | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'fg21sim/extragalactic') diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index bcda7bb..eb873ec 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -210,11 +210,10 @@ class RadioHalo1M: @lru_cache() def t_merger_end(self): """ - The time when the merger/turbulence ends. + The time when the merger (i.e., turbulence) ends. Unit: [Gyr] """ - tau_turb = self.duration_turb(self.t_merger) - return self.t_merger + tau_turb + return self.t_merger + self.duration_turb(self.t_merger) @property def radius(self): @@ -847,7 +846,7 @@ class RadioHaloAM(RadioHalo1M): Return the most recent merger event happend before the given time, i.e., the merger event that the given time locates in. """ - idx = (self.t_merger > t).sum() + idx = (self.t_merger > t).sum() # 't_merger' in decreasing order return { "idx": idx, "M_main": self.M_main[idx], @@ -908,6 +907,8 @@ class RadioHaloAM(RadioHalo1M): (i.e., smallest ``tau_acc``) is chosen and its beginning time is returned. Otherwise, the most recent merger event happened before the given time is chosen. + + Unit: [Gyr] """ mergers = [(tm, tm+self.duration_turb(tm), self.tau_acceleration(tm)) for tm in self.t_merger] -- cgit v1.2.2