aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile/extract_tcool.py
diff options
context:
space:
mode:
Diffstat (limited to 'mass_profile/extract_tcool.py')
-rwxr-xr-xmass_profile/extract_tcool.py42
1 files changed, 0 insertions, 42 deletions
diff --git a/mass_profile/extract_tcool.py b/mass_profile/extract_tcool.py
deleted file mode 100755
index f2b7774..0000000
--- a/mass_profile/extract_tcool.py
+++ /dev/null
@@ -1,42 +0,0 @@
-#!/usr/bin/env python3
-#
-# Extract the cooling time corresponding to the cooling radius.
-#
-# Junhua GU
-# 2012-12-20
-# Weitian LI
-# 2016-06-07
-#
-
-import argparse
-import numpy as np
-
-
-def get_tcool(data, rcool):
- """
- Get the cooling time *at* the specified cooling radius.
-
- XXX: whether to interpolate first?
- """
- radius = data[:, 0]
- ctime = data[:, 1]
- tcool = np.min(ctime[radius > rcool])
- return tcool
-
-
-def main():
- parser = argparse.ArgumentParser(
- description="Extract cooling time w.r.t the given cooling radius")
- parser.add_argument("infile", help="input cooling time data file")
- parser.add_argument("rcool", type=float, help="cooling radius (kpc)")
- args = parser.parse_args()
-
- data = np.loadtxt(args.infile)
- tcool = get_tcool(data, rcool=args.rcool)
- print("cooling time at %f kpc=%f Gyr" % (args.rcool, tcool))
- print("cooling_time= %f Gyr" % tcool)
- print("cooling_radius= %f kpc" % args.rcool)
-
-
-if __name__ == "__main__":
- main()