Commits

Matthew Turk committed 6e07409

Adding a few simple scripts

  • Participants
  • Parent commits 6a886d1

Comments (0)

Files changed (2)

File halo_find.py

+from mpi4py import MPI
+print MPI
+from yt.mods import *
+import glob
+
+fns = list(sorted(glob.glob("ic.enzo_512/*/*.hierarchy")))
+ds = TimeSeriesData.from_filenames(fns)
+for d in ds:
+    halos = amods.halo_finding.parallelHF(d)
+    halos.write_out("ic.enzo_512/%s.hop" % d)

File make_plots.py

+from mpi4py import MPI
+print MPI
+from yt.mods import *
+import glob
+
+fns = list(sorted(glob.glob("ic.enzo_512/*/*.hierarchy")))
+ds = TimeSeriesData.from_filenames(fns)
+for d in ds:
+    proj = d.h.proj(0, "Density")
+    frb = proj.to_frb(1.0, 2048)
+    write_image(na.log10(frb["Density"]), "frames/%s_Density.png" % d)