diff options
| author | Aaron LI <aly@aaronly.me> | 2017-06-24 21:52:41 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-06-24 21:52:41 +0800 |
| commit | f9f7a6e672367639cf1c3c93513f1a5e7e01063b (patch) | |
| tree | edbd1e2ed29d0fbe77307574e02f839345f71a93 /fg21sim | |
| parent | 374795fc794e7306cbd4260942f6a3051489e656 (diff) | |
| download | fg21sim-f9f7a6e672367639cf1c3c93513f1a5e7e01063b.tar.bz2 | |
halo.py: Add property "zgrid" as well as config option "zbinsize"
Diffstat (limited to 'fg21sim')
| -rw-r--r-- | fg21sim/configs/20-extragalactic.conf.spec | 4 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 8 |
2 files changed, 11 insertions, 1 deletions
diff --git a/fg21sim/configs/20-extragalactic.conf.spec b/fg21sim/configs/20-extragalactic.conf.spec index 4e84d85..7678a6a 100644 --- a/fg21sim/configs/20-extragalactic.conf.spec +++ b/fg21sim/configs/20-extragalactic.conf.spec @@ -48,6 +48,8 @@ # Maximum redshift until where to tracing the cluster merging history # (e.g., when calculating the electron spectrum) zmax = float(default=3.0, min=0.0) + # Redshift bin size for, e.g., calculating acceleration coefficients + zbinsize = float(default=0.01, min=0.0, max=0.1) # Mass threshold of the sub-cluster to be regarded as a significant # merger. (unit: Msun) @@ -59,7 +61,7 @@ # Radius of the giant radio halo in clusters (unit: kpc) # XXX: currently only support a constant radius of halos - radius = float(default=500, min=100) + radius = float(default=500.0, min=100.0) # Magnetic field scaling relation for clusters # Reference: Cassano et al. 2012, A&A, 548, A100, Eq.(1) diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 3fbd50c..8895b3a 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -103,6 +103,14 @@ class HaloSingle: logger.info("Loaded and set up configurations") @property + def zgrid(self): + """ + Redshift grid between 0 and ``self.zmax``. + """ + return np.arange(start=0.0, stop=self.zmax+self.zbinsize, + step=self.zbinsize) + + @property def magnetic_field(self): """ Calculate the mean magnetic field of this cluster from the |