#
# AMR PROBLEM DEFINITION FILE: Cosmology Simulation (amr version)
#
#  define problem
#
ProblemType                = 30      // cosmology simulation
TopGridRank                = 3
TopGridDimensions          = 512 512 512   
SelfGravity                = 1       // gravity on
TopGridGravityBoundary     = 0       // Periodic BC for gravity
LeftFaceBoundaryCondition  = 3 3 3   // same for fluid
RightFaceBoundaryCondition = 3 3 3
#
#  problem parameters
#
CosmologySimulationOmegaBaryonNow       = 0.044
CosmologySimulationOmegaCDMNow          = 0.226 
CosmologyOmegaMatterNow                 = 0.27 
CosmologyOmegaLambdaNow                 = 0.73  
#CosmologySimulationInitialTemperature   = 100.0         // K
CosmologySimulationDensityName          = GridDensity
CosmologySimulationVelocity1Name        = GridVelocities
CosmologySimulationVelocity2Name        = GridVelocities
CosmologySimulationVelocity3Name        = GridVelocities
CosmologySimulationParticlePositionName = ParticlePositions
CosmologySimulationParticleVelocityName = ParticleVelocities
CosmologySimulationNumberOfInitialGrids = 1
##CosmologySimulationSubgridsAreStatic    = 1
#
#  define cosmology parameters
#
ComovingCoordinates        = 1       // Expansion ON
CosmologyHubbleConstantNow = 0.71    // in km/s/Mpc
CosmologyComovingBoxSize   = 13.632  // in Mpc/h
CosmologyMaxExpansionRate  = 0.015   // maximum allowed delta(a)/a
CosmologyInitialRedshift   = 60      // 
CosmologyFinalRedshift     = 1.9     //
GravitationalConstant      = 1       // this must be true for cosmology
#
#  set I/O and stop/start parameters
#
#StopCycle              = 40          // stop after this many cycles
##dtDataDump             = 0
#CycleSkipDataDump      = 20
#DataDumpName           = 
CosmologyOutputRedshift[0] = 25.0 
CosmologyOutputRedshift[1] = 10.0
CosmologyOutputRedshift[2] = 4.0
CosmologyOutputRedshift[3] = 3.0
CosmologyOutputRedshift[4] = 2.1
CosmologyOutputRedshift[5] = 2.0
CosmologyOutputRedshift[6] = 1.9
#
#  set hydro parameters
#
Gamma                  = 1.6667
PPMDiffusionParameter  = 0       // diffusion off
DualEnergyFormalism    = 1       // use total & internal energy
InterpolationMethod    = 1     // SecondOrderA
CourantSafetyNumber    = 0.5
ParticleCourantSafetyNumber = 0.8
FluxCorrection         = 1
ConservativeInterpolation = 0
HydroMethod            = 0
#
#  set cooling parameters
#
RadiativeCooling       = 1
MultiSpecies           = 1
RadiationFieldType     = 12
StarParticleCreation   = 0
StarParticleFeedback   = 0
#
#  set grid refinement parameters
#
StaticHierarchy           = 1    // Not dynamic hierarchy
ParallelParticleIO	  = 1    // for use with...
ParallelRootGridIO	  = 1    // ... poe 
##MaximumRefinementLevel    = 8    // use up to this many levels
##RefineBy                  = 2    // refinement factor
##CellFlaggingMethod        = 2 4  // use baryon mass for refinement
##MinimumEfficiency         = 0.4  // fraction efficiency
##MinimumOverDensityForRefinement = 4.5 9.0 // times the initial density
##MinimumMassForRefinementLevelExponent = 0.0 0.0
##MinimumEnergyRatioForRefinement = 0.4 // min Egas/Etot for shock refinement
##RefineRegionLeftEdge            = 0.152 0.164 0.303
##RefineRegionRightEdge           = 0.396 0.398 0.466
#
#  set some global parameters
#
GreensFunctionMaxNumber   = 100   // # of greens function at any one time
##UseMinimumPressureSupport = 0