1 files changed, 14 insertions, 0 deletions
diff --git a/mod_ly/extract_tcool.py b/mod_ly/extract_tcool.py
new file mode 100755
index 0000000..05c836b
--- /dev/null
+++ b/mod_ly/extract_tcool.py
@@ -0,0 +1,14 @@
+#!/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))
+        print("rcool= %f kpc" % rcool)
+        print("cooling_time= %f Gyr @ %f kpc" % (t,rcool))
+        sys.exit(0)