From 00979a4da68b1fc3bf5c1693975febc7d066763c Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Thu, 8 Nov 2018 15:52:37 +0800 Subject: clusters/halo: Add and use 'f_rc' and 'beta' options --- fg21sim/configs/config.spec | 7 +++++++ fg21sim/extragalactic/clusters/halo.py | 5 ++++- 2 files changed, 11 insertions(+), 1 deletion(-) diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec index fefb8b8..85606eb 100644 --- a/fg21sim/configs/config.spec +++ b/fg21sim/configs/config.spec @@ -447,3 +447,10 @@ stream = option("stderr", "stdout", "", default="stderr") # Unit: [MHz] fiducial_freq = float(default=150.0) fiducial_factor = float(default=3.0, min=1) + + # Parameters of the beta-model that is used to describe the gas density + # profile of the cluster. + # The fraction of the core radius to cluster's virial radius. + f_rc = float(default=0.1) + # The slope parameter (i.e., beta). + beta = float(default=0.8) diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 113e434..a18aa17 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -167,6 +167,8 @@ class RadioHalo: self.injection_index = configs.getn(comp+"/injection_index") self.fiducial_freq = configs.getn(comp+"/fiducial_freq") self.fiducial_factor = configs.getn(comp+"/fiducial_factor") + self.f_rc = configs.getn(comp+"/f_rc") + self.beta = configs.getn(comp+"/beta") def _set_solver(self): self.fpsolver = FokkerPlanckSolver( @@ -673,7 +675,8 @@ class RadioHalo: z = COSMO.redshift(t) M_main = self.mass_main(t) M_sub = self.mass_sub(t) - return helper.calc_gas_density_profile(mass=M_main+M_sub, z=z) + return helper.calc_gas_density_profile(mass=M_main+M_sub, z=z, + f_rc=self.f_rc, beta=self.beta) @lru_cache() def _velocity_turb(self, t): -- cgit v1.2.2