Source

bluewaters.support_files / enzo / 35Mpc_1024_2_nested.enzo

The default branch has multiple heads

#
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