summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2016-07-10 22:03:27 +0800
committerAaron LI <aaronly.me@outlook.com>2016-07-10 22:03:27 +0800
commit2bde3651a56296b3b2b9b5c7471244ab7d7478fb (patch)
treed753ec5492ee7c9de74374d9c0fd2636bc32c69e
parent2f89d7da79e4d7984575d12d88d8d42b44d45d05 (diff)
downloadcexcess-2bde3651a56296b3b2b9b5c7471244ab7d7478fb.tar.bz2
fit_tprofile.py: set default "rcut=3000"
-rwxr-xr-xfit_tprofile.py13
1 files changed, 9 insertions, 4 deletions
diff --git a/fit_tprofile.py b/fit_tprofile.py
index 40d888c..88f7fa7 100755
--- a/fit_tprofile.py
+++ b/fit_tprofile.py
@@ -1,7 +1,12 @@
#!/usr/bin/env python3
#
-# Weitian LI
-# 2016-07-04
+# Aaron LI
+# Created: 2016-07-04
+# Updated: 2016-07-04
+#
+# Change logs:
+# 2016-07-04:
+# * Set default "rcut=3000" for TemperatureProfile.extrapolate()
#
"""
@@ -25,7 +30,7 @@ redshift = <REDSHIFT>
t_profile_data = t_profile_data.txt
# cut radius to which stop the extrapolation (unit: kpc)
-rcut_extrap = <RCUT>
+rcut_extrap = 3000
# number of data points for the output temperature profile
num_dp = 1000
@@ -155,7 +160,7 @@ class TemperatureProfile:
xerr=self.r_err, yerr=self.t_err)
self.model.fit()
- def extrapolate(self, rcut, num=1000):
+ def extrapolate(self, rcut=3000, num=1000):
"""
Interpolate and extrapolate the fitted temperature profile.