diff options
| author | Aaron LI <aly@aaronly.me> | 2018-11-08 15:52:37 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2018-11-08 15:52:37 +0800 |
| commit | 00979a4da68b1fc3bf5c1693975febc7d066763c (patch) | |
| tree | 4fefb234a7517c0c4438930778b587a95547c3e5 | |
| parent | 450767336f1d37e38f3e71c40695ac4a3cd20162 (diff) | |
| download | fg21sim-00979a4da68b1fc3bf5c1693975febc7d066763c.tar.bz2 | |
clusters/halo: Add and use 'f_rc' and 'beta' options
| -rw-r--r-- | fg21sim/configs/config.spec | 7 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 5 |
2 files changed, 11 insertions, 1 deletions
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): |