diff options
| -rwxr-xr-x | mass_profile/extract_tcool.py | 42 |
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() |