diff options
| author | Aaron LI <aaronly.me@outlook.com> | 2017-01-07 16:52:28 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-06-01 16:33:38 +0800 |
| commit | 9d1fbf8d875d4ed754cbf6894109e3379c98b0bb (patch) | |
| tree | 0ac2a221d10ac22e340d6c5f0f5f9bb550dd7ce7 | |
| parent | a154d19d38fc0b693184ac1b2edbd0565fdcec55 (diff) | |
| download | fg21sim-9d1fbf8d875d4ed754cbf6894109e3379c98b0bb.tar.bz2 | |
clusters/halo: Add missing config "buffer_np"
| -rw-r--r-- | fg21sim/configs/20-extragalactic.conf.spec | 5 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 2 |
2 files changed, 7 insertions, 0 deletions
diff --git a/fg21sim/configs/20-extragalactic.conf.spec b/fg21sim/configs/20-extragalactic.conf.spec index 5c1d3aa..633ef78 100644 --- a/fg21sim/configs/20-extragalactic.conf.spec +++ b/fg21sim/configs/20-extragalactic.conf.spec @@ -77,6 +77,11 @@ # equation to calculate the electron spectrum. pgrid_num = integer(default=100, min=10) + # Number of grid points used as the buffer region near the lower + # boundary, and the value within this buffer region will be fixed to + # avoid unphysical pile-up of low-energy electrons. + buffer_np = integer(default=5, min=0) + # Time step for solving the Fokker-Planck equation (unit: Gyr) time_step = float(default=0.01, min=1e-5, max=1.0) diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 26ff6de..072cdff 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -98,6 +98,7 @@ class HaloSingle: self.pmin = self.configs.getn(comp+"/pmin") self.pmax = self.configs.getn(comp+"/pmax") self.pgrid_num = self.configs.getn(comp+"/pgrid_num") + self.buffer_np = self.configs.getn(comp+"/buffer_np") self.time_step = self.configs.getn(comp+"/time_step") self.injection_index = self.configs.getn(comp+"/injection_index") # Cosmology model @@ -121,6 +122,7 @@ class HaloSingle: fpsolver = FokkerPlanckSolver( xmin=self.pmin, xmax=self.pmax, grid_num=self.pgrid_num, + buffer_np=self.buffer_np, tstep=self.time_step, f_advection=self.fp_advection, f_diffusion=self.fp_diffusion, |