Commits

Matthew Turk committed 73462ae

Adding more redshift outputs and correcting a bug in the MUSIC parameter files.

Comments (0)

Files changed (5)

1200kpc_128_nested.conf

 transfer		= eisenstein
 
 [random]
-seed[7]		= 122140
-seed[8]		= 480428
+seed[7]		= 480428
+seed[8]		= 122140
 seed[9]		= 800208
 seed[10]	= 190184
 seed[11]	= 650517

128_nested_static.qsub

 #!/bin/bash
 #PBS -N nibblonia-128-nested
-#PBS -l nodes=16:ppn=16
-#PBS -l walltime=72:00:00
+#PBS -l nodes=8:ppn=32
+#PBS -l walltime=336:00:00
 #PBS -o nibblonia-128-nested.log
 #PBS -e nibblonia-128-nested.log
 #PBS -V
 #PBS -A col117
 #PBS -m abe
 
-export NP=512
+export NP=256
 export MX_RCACHE=0
 export MPI_DIR=/home/diag/opt/mvapich2/1.5.1p1/intel/
 cd $PBS_O_WORKDIR/ic.enzo_full
 from yt.mods import *
 import glob
 
-fns = list(sorted(glob.glob("ic.enzo_512/*/*.hierarchy")))
+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_512/%s.hop" % 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())
 #$MPI_DIR/bin/mpirun -hostfile $PBS_NODEFILE -np $NP ./ring.exe pv ParticlePositions.3 ParticleVelocities.3 |& tee -a estd.out
 #$MPI_DIR/bin/mpirun -hostfile $PBS_NODEFILE -np $NP ./enzo.exe -d FirstStar_WMAP7_Real1_AMR.enzo |& tee -a estd.out
 #./MUSIC 1200kpc_512_unigrid.conf |& tee -a MUSIC.out
-setenv OMP_NUM_THREADS 16
+setenv OMP_NUM_THREADS 32
 ./MUSIC 1200kpc_128_nested.conf |& tee -a MUSIC_full.out
 

parameter_file_128_nested_static.enzo

 #
 # region allowed for further refinement
 #
-RefineRegionAutoAdjust                   = 1
+RefineRegionAutoAdjust                   = 0
 RefineRegionLeftEdge                     =         0.484375        0.484375        0.484375
 RefineRegionRightEdge                     =         0.515625         0.515625         0.515625
 
 RedshiftDumpName = RD
 DataDumpName = DD
+#dtDataDump = 1.0
+MultiSpecies = 0
 
-CosmologyOutputRedshift[00] = 30
-CosmologyOutputRedshift[01] = 25
-CosmologyOutputRedshift[02] = 24
-CosmologyOutputRedshift[03] = 23
-CosmologyOutputRedshift[04] = 22
-CosmologyOutputRedshift[05] = 21
-CosmologyOutputRedshift[06] = 20
-CosmologyOutputRedshift[07] = 19
-CosmologyOutputRedshift[08] = 18
-CosmologyOutputRedshift[09] = 17
-CosmologyOutputRedshift[10] = 16
-CosmologyOutputRedshift[11] = 15
-CosmologyOutputRedshift[12] = 14
-CosmologyOutputRedshift[13] = 13
+CosmologyOutputRedshift[00] = 99
+CosmologyOutputRedshift[01] = 90
+CosmologyOutputRedshift[02] = 80
+CosmologyOutputRedshift[03] = 70
+CosmologyOutputRedshift[04] = 60
+CosmologyOutputRedshift[05] = 50
+CosmologyOutputRedshift[06] = 40
+CosmologyOutputRedshift[07] = 30
+CosmologyOutputRedshift[08] = 28
+CosmologyOutputRedshift[09] = 27
+CosmologyOutputRedshift[10] = 26
+CosmologyOutputRedshift[11] = 25
+CosmologyOutputRedshift[12] = 24
+CosmologyOutputRedshift[13] = 23
+CosmologyOutputRedshift[14] = 22
+CosmologyOutputRedshift[15] = 21
+CosmologyOutputRedshift[16] = 20
+CosmologyOutputRedshift[17] = 19
+CosmologyOutputRedshift[18] = 18
+CosmologyOutputRedshift[19] = 17
+CosmologyOutputRedshift[20] = 16
+CosmologyOutputRedshift[21] = 15
+CosmologyOutputRedshift[22] = 14
+CosmologyOutputRedshift[23] = 13
 
 CellFlaggingMethod = 6
+HydroMethod = 0
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.