mod_ly/extract_tcool.py


1
2
3
4
5
6
7
8
9
10
11
12
13
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)