Commits

Matthew Turk  committed 4c55289

Trying to get the ICs correct.

  • Participants
  • Parent commits 73462ae

Comments (0)

Files changed (4)

File 1200kpc_128_nested.conf

 levelmax		= 14
 padding			= 8
 overlap			= 16
-ref_offset		= 0.1953125,0.46875,0.59375
+ref_center		= 0.19635944, 0.4712808, 0.59101587
+#ref_offset		= 0.1953125,  0.46875  ,  0.5859375
 ref_extent		= 0.03125,0.03125,0.03125
 align_top		= yes
 baryons			= yes

File 1200kpc_512_unigrid.conf

 levelmax		= 9
 padding			= 8
 overlap			= 16
-ref_offset		= 0.4, 0.4, 0.4
+ref_offset		= 0.0, 0.0, 0.0
 ref_extent		= 0.2, 0.2, 0.2
 align_top		= yes
 baryons			= yes

File 128_nested_static.qsub

 #!/bin/bash
 #PBS -N nibblonia-128-nested
 #PBS -l nodes=8:ppn=32
-#PBS -l walltime=336:00:00
+#PBS -l walltime=12:00:00
 #PBS -o nibblonia-128-nested.log
 #PBS -e nibblonia-128-nested.log
 #PBS -V

File halo_find.py

 from mpi4py import MPI
 print MPI
-from yt.mods import *
+from yt.pmods 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")))
+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_full/%s.hop" % d)
-    halos.write_particle_lists("ic.enzo_full/%s.particles" % d)
+    halos.write_out("ic.enzo_512/%s.hop" % d)
+    halos.write_particle_lists("ic.enzo_512/%s.particles" % d)
     for halo in halos:
         print "TOTAL MASS: %0.3e, %s" % (
             halo.total_mass(), halo.center_of_mass())