Commits

Britton Smith committed d40a325

Adding enzo parameter file for the plan B production run (1024^3 root grid
with 3 nested level and 1536^3 most refined region.

  • Participants
  • Parent commits e6e903c

Comments (0)

Files changed (1)

enzo/35Mpc_1024_3_nested_planB.enzo

+#
+ProblemType                              = 30      // cosmology simulation
+TopGridRank                              = 3
+TopGridDimensions                        = 1024 1024 1024
+SelfGravity                              = 1       // gravity on
+TopGridGravityBoundary                   = 0       // Periodic BC for gravity
+LeftFaceBoundaryCondition                = 3 3 3   // same for fluid
+RightFaceBoundaryCondition               = 3 3 3
+
+#
+# problem specific parameters
+#
+CosmologySimulationOmegaBaryonNow        = 0.0455
+CosmologySimulationOmegaCDMNow           = 0.2265
+CosmologySimulationDensityName           = GridDensity
+CosmologySimulationVelocity1Name         = GridVelocities_x
+CosmologySimulationVelocity2Name         = GridVelocities_y
+CosmologySimulationVelocity3Name         = GridVelocities_z
+CosmologySimulationCalculatePositions    = 1
+CosmologySimulationParticleVelocity1Name = ParticleVelocities_x
+CosmologySimulationParticleVelocity2Name = ParticleVelocities_y
+CosmologySimulationParticleVelocity3Name = ParticleVelocities_z
+CosmologySimulationParticleDisplacement1Name = ParticleDisplacements_x
+CosmologySimulationParticleDisplacement2Name = ParticleDisplacements_y
+CosmologySimulationParticleDisplacement3Name = ParticleDisplacements_z
+CosmologySimulationNumberOfInitialGrids  = 4
+CosmologySimulationGridDimension[1]      =              412              412              412
+CosmologySimulationGridLeftEdge[1]       =      0.400390625      0.400390625      0.400390625
+CosmologySimulationGridRightEdge[1]      =        0.6015625        0.6015625        0.6015625
+CosmologySimulationGridLevel[1]          = 1
+CosmologySimulationGridDimension[2]      =              788              788              788
+CosmologySimulationGridLeftEdge[2]       =      0.404296875      0.404296875      0.404296875
+CosmologySimulationGridRightEdge[2]      =     0.5966796875     0.5966796875     0.5966796875
+CosmologySimulationGridLevel[2]          = 2
+CosmologySimulationGridDimension[3]      =             1536             1536             1536
+CosmologySimulationGridLeftEdge[3]       =          0.40625          0.40625          0.40625
+CosmologySimulationGridRightEdge[3]      =          0.59375          0.59375          0.59375
+CosmologySimulationGridLevel[3]          = 3
+CosmologySimulationUseMetallicityField   = 1
+
+#
+# cosmology parameters
+#
+ComovingCoordinates                      = 1       // Expansion ON
+CosmologyOmegaMatterNow                  = 0.272
+CosmologyOmegaLambdaNow                  = 0.728
+CosmologyHubbleConstantNow               = 0.704   // in 100 km/s/Mpc
+CosmologyComovingBoxSize                 = 35.0      // in Mpc/h
+CosmologyMaxExpansionRate                = 0.015   // maximum allowed delta(a)/a
+CosmologyInitialRedshift                 = 99.0      //
+CosmologyFinalRedshift                   = 8.9     // go to here first for cooling
+#CosmologyFinalRedshift                  = 6      // restart with this
+GravitationalConstant                    = 1.0       // this must be true for cosmology
+
+#
+# I/O parameters
+#
+dtDataDump                               = 0.2299  // 1.0 in code units = 21.75 MYrs
+DataDumpName                             = DD
+DataDumpDir                              = DD
+RedshiftDumpName                         = RD
+RedshiftDumpDir                          = RD
+OutputCoolingTime                        = 1
+OutputTemperature                        = 1
+
+CosmologyOutputRedshift[0]               = 99.0
+CosmologyOutputRedshift[1]               = 50.0
+CosmologyOutputRedshift[2]               = 40.0
+CosmologyOutputRedshift[3]               = 35.0
+CosmologyOutputRedshift[4]               = 30.0
+CosmologyOutputRedshift[5]               = 28.0
+CosmologyOutputRedshift[6]               = 26.0
+CosmologyOutputRedshift[7]               = 24.0
+CosmologyOutputRedshift[8]               = 22.0
+CosmologyOutputRedshift[9]               = 20.0
+CosmologyOutputRedshift[10]               = 19.0
+CosmologyOutputRedshift[11]               = 18.0
+CosmologyOutputRedshift[12]               = 17.0
+CosmologyOutputRedshift[13]               = 16.0
+CosmologyOutputRedshift[14]               = 15.0
+CosmologyOutputRedshift[15]               = 14.0
+CosmologyOutputRedshift[16]               = 13.0
+CosmologyOutputRedshift[17]               = 12.0
+CosmologyOutputRedshift[18]               = 11.0
+CosmologyOutputRedshift[19]               = 10.0
+CosmologyOutputRedshift[20]               = 9.0
+CosmologyOutputRedshift[21]               = 8.9
+CosmologyOutputRedshift[22]               = 8.0
+CosmologyOutputRedshift[23]               = 7.0
+CosmologyOutputRedshift[24]               = 6.0
+
+#
+# set hydro parameters
+#
+HydroMethod                              = 0 
+Gamma                                    = 1.6667
+DualEnergyFormalism                      = 1
+ReconstructionMethod                     = 1 // PPM
+RiemannSolver                            = 4 // HLLC
+RiemannSolverFallback                    = 1 // HLL if there's a dnu<=0
+InterpolationMethod                      = 1
+CourantSafetyNumber                      = 0.4
+ParticleCourantSafetyNumber              = 0.8
+UseMinimumPressureSupport                = 1
+MinimumPressureSupportParameter          = 64.0
+FluxCorrection                           = 1
+ConservativeInterpolation                = 0
+
+#
+# refinement parameters
+#
+StaticHierarchy                          = 0
+RefineBy                                 = 2 
+MaximumRefinementLevel                   = 6
+CellFlaggingMethod                       = 2 4
+MinimumEfficiency                        = 0.4
+MinimumOverDensityForRefinement          = 0.0078125 0.0078125
+MinimumMassForRefinementLevelExponent    = -0.1 -0.1
+NumberOfBufferZones                      = 1 
+RefineRegionAutoAdjust                   = 0
+RefineRegionLeftEdge                     =          0.40625         0.40625         0.40625
+RefineRegionRightEdge                    =          0.59375          0.59375          0.59375
+SubgridSizeAutoAdjust = 1  // gives best performance
+
+# LB = 4 hilbert curve
+LoadBalancing = 4
+
+#
+#  chemistry, cooling, radiation background
+#
+RadiativeCooling                         = 1
+MultiSpecies                             = 2
+MetalCooling                             = 3 // Cloudy cooling
+CloudyCoolingGridFile                    = /u/sciteam/brittons/cloudy_cooling_data/solar_2009_4D_metals.h5 // until z = 8.9, then 
+                                                               // solar_2009_5D_metals.h5
+IncludeCloudyHeating                     = 0 // until z = 8.9, then 1.
+CMBTemperatureFloor                      = 1
+RadiationFieldType                       = 13
+AdjustUVBackground                       = 0
+RadiationRedshiftOn                      = 15
+RadiationRedshiftFullOn                  = 10
+
+#
+# star formation and feedback parameters
+#
+StarParticleCreation                     = 1
+StarParticleFeedback                     = 1
+StarMakerUsePhysicalDensityThreshold     = 1
+StarMakerOverDensityThreshold            = 10.0
+StarMakerMinimumMass                     = 3.0e+5  # recalculate this
+StarMakerMinimumDynamicalTime            = 1e6
+StarMakerMassEfficiency                  = 0.01
+StarMassEjectionFraction                 = 0.25
+StarMetalYield                           = 0.025
+StarEnergyToThermalFeedback              = 1e-5
+StarFeedbackDistRadius                   = 1
+StarFeedbackDistCellStep                 = 3
+
+#
+# root grid i/o
+#
+Unigrid                                  = 1
+ParallelRootGridIO                       = 1
+ParallelParticleIO                       = 1
+PartitionNestedGrids                     = 1
+
+ReadGhostZones                           = 1
+WriteGhostZones                          = 1
+ParticleTypeInFile                       = 1 # particle type field
+ParticleTypeGrouping                     = 1 # H5S for different particle types
+HierarchyFileOutputFormat                = 2