Add mass_profile tools

* These tools are mainly use to calculate the total gravitational mass
  profile, as well as the intermediate products (e.g., surface
  brightness profile fitting, gas density profile, NFW fitting, etc.)
* There are additional tools for calculating the luminosity and flux.
* These tools mainly developed by Junhua GU, and contributed by
  Weitian (Aaron) LI, and Zhenghao ZHU.

1 files changed, 12 insertions, 0 deletions

diff --git a/mass_profile/extract_tcool.py b/mass_profile/extract_tcool.py
new file mode 100755
index 0000000..393af27
--- /dev/null
+++ b/mass_profile/extract_tcool.py
@@ -0,0 +1,12 @@
+#!/usr/bin/env python
+
+import sys
+rcool=float(sys.argv[1])
+
+for l in open('cooling_time.dat'):
+    r,t=l.split()
+    r=float(r)
+    t=float(t)
+    if r>rcool:
+        print("cooling time at %f kpc=%f Gyr"%(rcool,t))
+        sys.exit(0)