nibblonia / halo_find.py

from mpi4py import MPI
print MPI
from yt.mods import *
import glob

fns = list(sorted(glob.glob("ic.enzo_full/RD001[6-9]/RD001[6-9].hierarchy")
                + glob.glob("ic.enzo_full/RD002[0-4]/RD002[0-4].hierarchy")))
ds = TimeSeriesData.from_filenames(fns)
for d in ds:
    halos = amods.halo_finding.parallelHF(d)
    halos.write_out("ic.enzo_full/%s.hop" % d)
    halos.write_particle_lists("ic.enzo_full/%s.particles" % d)
    for halo in halos:
        print "TOTAL MASS: %0.3e, %s" % (
            halo.total_mass(), halo.center_of_mass())
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.