diff options
| author | Aaron LI <aly@aaronly.me> | 2018-01-01 16:57:16 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2018-01-01 16:57:16 +0800 |
| commit | c9b0744c1056ebcaf425c4ff7d6eb98786b40699 (patch) | |
| tree | a87b025f3f8ac1a86de3afb4436ca59f6d166396 | |
| parent | 151dafecd3ae38c62d0fa2e10ff2160b7b17ff5c (diff) | |
| download | fg21sim-c9b0744c1056ebcaf425c4ff7d6eb98786b40699.tar.bz2 | |
clusters/halo: impose the maximum acceleration timescale
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 12 |
1 files changed, 9 insertions, 3 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index a8f218a..1da89b3 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -1,4 +1,4 @@ -# Copyright (c) 2017 Weitian LI <weitian@aaronly.me> +# Copyright (c) 2017-2018 Weitian LI <weitian@aaronly.me> # MIT license """ @@ -746,12 +746,18 @@ class RadioHalo: ---------- Ref.[donnert2013],Eq.(15) """ + # Maximum acceleration timescale when no turbulence acceleration + # NOTE: see the above WARNING! + tau_max = 10.0 # [Gyr] if (t < self.age_merger) or (t > self.age_merger+self.time_turbulence): - # NO acceleration (see also the above NOTE and WARNING!) - tau_acc = 10 # [Gyr] + # NO active turbulence acceleration + tau_acc = tau_max else: # Turbulence acceleration tau_acc = self.tau_acceleration # [Gyr] + # Impose the maximum acceleration timescale + if tau_acc > tau_max: + tau_acc = tau_max gamma = np.asarray(gamma) diffusion = gamma**2 / 4 / tau_acc |