#
# Method paper parameter study parameter file
#
ProblemType                = 30
TopGridRank                = 3
SelfGravity                = 1
TopGridGravityBoundary     = 0
LeftFaceBoundaryCondition  = 3 3 3
RightFaceBoundaryCondition = 3 3 3
BaryonSelfGravityApproximation = 1
TopGridDimensions = 1024 1024 1024
#
#  problem parameters
#
CosmologySimulationOmegaBaryonNow       = 0.0478
CosmologySimulationOmegaCDMNow          = 0.2602
CosmologySimulationNumberOfInitialGrids = 1
CosmologySimulationDensityName          = GridDensity
CosmologySimulationVelocity1Name        = GridVelocities_x
CosmologySimulationVelocity2Name        = GridVelocities_y
CosmologySimulationVelocity3Name        = GridVelocities_z
CosmologySimulationParticleVelocity1Name = ParticleVelocities_x
CosmologySimulationParticleVelocity2Name = ParticleVelocities_y
CosmologySimulationParticleVelocity3Name = ParticleVelocities_z
CosmologySimulationCalculatePositions = 1

#  define cosmology parameters
#
ComovingCoordinates        = 1
CosmologyOmegaMatterNow    = 0.308
CosmologyOmegaLambdaNow    = 0.692
CosmologyHubbleConstantNow = 0.678
CosmologyComovingBoxSize   = 57.83
CosmologyMaxExpansionRate  = 0.02
CosmologyInitialRedshift   = 40
CosmologyFinalRedshift     = 7.9
GravitationalConstant      = 1

#
#  set I/O and stop/start parameters
#
CycleSkipDataDump = 15
dtDataDump             = 0.0
RedshiftDumpName = csf_B0RD
RedshiftDumpDir = csf_B0RD
DataDumpName = csf_B0DD
DataDumpDir = csf_B0DD


CosmologyOutputRedshift[0] = 8
CosmologyOutputRedshift[1] = 6
CosmologyOutputRedshift[2] = 4
CosmologyOutputRedshift[3] = 3
CosmologyOutputRedshift[4] = 2
CosmologyOutputRedshift[5] = 1
CosmologyOutputRedshift[6] = 0.5
CosmologyOutputRedshift[7] = 0.02

#
#  set hydro parameters
#
Gamma                  = 1.6667
PPMDiffusionParameter  = 0
DualEnergyFormalism    = 1
ParticleCourantSafetyNumber    = 0.6
UseMinimumPressureSupport = 0

MultiSpecies         = 1
MultiMetals = 1
RadiativeCooling     = 1
RadiationFieldType        = 1
RiemannSolver =1
RiemannSolverFallback = 1
ReconstructionMethod = 0
ParticleSubgridDepositMode     = 0
MetalCooling = 0
CloudyCoolingGridFile     = solar_2008_3D_metals.h5

# star particle parameters^M
#^M
StarParticleCreation          = 0
StarParticleFeedback          = 0
OutputGriddedStarParticle = 0
StarMakerOverDensityThreshold = 0.001
StarMakerSHDensityThreshold           = 7e-26
StarMakerMinimumDynamicalTime = 1.5e+09
StarMakerMinimumMass          = 1e7
StarMakerMassEfficiency       = 0.01
StarMassEjectionFraction      = 0.25
StarMetalYield                = 0.025
StarEnergyToThermalFeedback   = 1e-8
StarFeedbackDistRadius        = 0
StarFeedbackDistCellStep      = 0

MBHInsertLocationFilename  =  /scratch/daint/vazzaf/testC/BH_location_z4.in

MBHAccretion = 0
MBHAccretionRadius = 50
MBHAccretingMassRatio = 1
MBHMinDynamicalTime = 1e3
MBHMinimumMass = 1
MBHParticleIO  = 1


UseDivergenceCleaning = 0
DualEnergyFormalism = 1
HydroMethod            = 4
CourantSafetyNumber    = 0.4
Theta_Limiter=1.0
InterpolationMethod = 1
CosmologySimulationInitialUniformBField = 1e-6 1e-6 1e-6

ParticleCourantSafetyNumber    = 0.6
SetUVBAmplitude                = 1
SetHeIIHeatingScale            = 1.8
RadiationFieldLevelRecompute   = 0
RadiationSpectrumNormalization = 1e-21
RadiationSpectrumSlope         = 1.5
FluxCorrection            = 1
ConservativeInterpolation = 0

agn= 0
ijet= 4
agn_dens= 100
k0= 0.01
agn_B= 0.01
sf=1
sf_B= 0.00
sf_dens=0.00
sf_z=8.00
agn_z= 10

#
#  set grid refinement parameters
#
StaticHierarchy           = 1
MaximumRefinementLevel    = 0
MaximumGravityRefinementLevel = 0
RefineBy                  = 2
CellFlaggingMethod        = 2 3
MinimumEfficiency         = 0.2
MinimumOverDensityForRefinement = 2
MinimumPressureJumpForRefinement = 5
RefineRegionLeftEdge =  0.3 0.4 0.3125
RefineRegionRightEdge = 0.55 0.8125 0.4125
MinimumMassForRefinementLevelExponent = -0.3
MinimumEnergyRatioForRefinement = 0.4
# some global parameters
#

GreensFunctionMaxNumber   = 40
ParallelRootGridIO        = 1
ParallelParticleIO        = 1
Unigrid                   = 1
UseCUDA = 1