diff options
author | Aaron LI <aly@aaronly.me> | 2019-01-18 00:01:25 +0800 |
---|---|---|
committer | Aaron LI <aly@aaronly.me> | 2019-01-18 13:00:52 +0800 |
commit | 4a5cbc5147ae4b0330e3154d5e4dda38601ecadd (patch) | |
tree | f59f0751cb0f823171146fe358bed97cbb9fc59e /fg21sim/extragalactic/clusters/helper.py | |
parent | 0abeded3cad5de40789c3d9bbe253aac6af8927f (diff) | |
download | fg21sim-4a5cbc5147ae4b0330e3154d5e4dda38601ecadd.tar.bz2 |
Minor tweaks
Diffstat (limited to 'fg21sim/extragalactic/clusters/helper.py')
-rw-r--r-- | fg21sim/extragalactic/clusters/helper.py | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py index 4d06ceb..e27d763 100644 --- a/fg21sim/extragalactic/clusters/helper.py +++ b/fg21sim/extragalactic/clusters/helper.py @@ -137,7 +137,7 @@ def radius_virial(mass, z=0.0): return R_vir * AUC.cm2kpc # [kpc] -def radius_stripping(M_main, M_sub, z, configs=CONFIGS): +def radius_stripping(M_main, M_sub, z, f_rc=0.1): """ Calculate the stripping radius of the in-falling sub-cluster, which is determined by the equipartition between the static and ram pressure. @@ -162,7 +162,7 @@ def radius_stripping(M_main, M_sub, z, configs=CONFIGS): rs = optimize.brentq(lambda r: f_rho_sub(r) - rhs, a=0.1*r_vir, b=r_vir, xtol=1e-1) except ValueError: - rs = 0.1 * r_vir + rs = 2*f_rc * r_vir return rs # [kpc] |