ProblemTypeName            = AgoraRestart      // Agora Isolated Galaxy Simulation
TopGridRank                = 3
TopGridDimensions          = 128 128 128
SelfGravity                = 0       // gravity off
TopGridGravityBoundary     = 1       // isolated gravity BCs
LeftFaceBoundaryCondition  = 3 3 3   // periodic
RightFaceBoundaryCondition = 3 3 3
DomainLeftEdge             = 0 0 0
DomainRightEdge            = 1 1 1
CycleSkipDataDump          = 1
StopCycle                  = 2

#
# problem parameters
#
# All parameters are in code units except as noted.
AgoraRestartCenterPosition           = 0.5 0.5 0.5 // code units
# The following physical parameters are completely non-physical -- just used for
# getting a test case which reproduces the metallicity issue
AgoraRestartScaleLength              = .3288817978
AgoraRestartScaleHeight              = .0325288817978
AgoraRestartDiskMass                 = 1900.79142
AgoraRestartGasFraction              = 0.9
AgoraRestartDiskTemperature          = 1.0e4     // Kelvin
AgoraRestartHaloMass                 = 77.8844181 // 10^-6 cm^-3
AgoraRestartHaloTemperature          = 1.0e6
AgoraRestartRefineAtStart            = 1
AgoraRestartUseMetallicityField      = 1
AgoraRestartCutOffZ                  = 0.114440899
AgoraRestartCutOffRadius             = 0.9440899
AgoraRestartMetalFractionByMass      = 1.0
AgoraRestartDiskMetallicity          = 1.0
AgoraRestartHaloMetallicity          = 1.0

#
#  define cosmology parameters
#
ComovingCoordinates        = 0       // Expansion OFF

#
# units
#
MassUnits          = 1.9891e42     // 10^9 msun in grams
LengthUnits        = 4.04446e24    // 1.31072 mpc (chosen to get dx = 80 pc on level 8) ???
TimeUnits          = 3.15569e16    // 1 Gyr in s
GravitationalConstant      = 0.0000251055  // 4*pi*G*DensityUnits*TimeUnits^2

#
#  set I/O and stop/start parameters
#
dtDataDump             = 0.3
DataDumpName           = DD
ParallelRootGridIO     = 0

Initialdt = 1.0e-6

#
#  set hydro parameters
#
Gamma                       = 1.6667
PPMDiffusionParameter       = 0       // diffusion off
DualEnergyFormalism         = 1
DualEnergyFormalismEta1     = .001 // default
DualEnergyFormalismEta2     = .4 // default
CourantSafetyNumber         = 0.1
RadiativeCooling            = 0
MultiSpecies                = 0
FluxCorrection              = 1
ConservativeInterpolation   = 1
RiemannSolver               = 5  // two shock
CorrectParentBoundaryFlux   = 1
HydroMethod                 = 0  // PPM
Mu              = 1.4
OutputTemperature       = 1

#
#  set grid refinement parameters
#
StaticHierarchy           = 0      // dynamic hierarchy
MaximumRefinementLevel    = 0
RefineBy                  = 2      // refinement factor
CellFlaggingMethod        = 2 4 6     // gas mass, particle mass, jeans length, metallicity, must-refine particles, slope
RefineByJeansLengthSafetyFactor = 4
MinimumMassForRefinement =  .0003 .0003 0 0 0
MinimumMassForRefinementLevelExponent  =  -.937403 -0.444916 0 0 0
MustRefineRegionLeftEdge = 0.475 0.475 0.475
MustRefineRegionRightEdge = 0.525 0.525 0.525
MustRefineRegionMinRefinementLevel = 5
SecondDerivativeFlaggingFields = 20  // density, metallicity. --- see typedefs.h
SecondDerivativeForRefinement = 1.99
SecondDerivativeEpsilon = 0.1
SlopeFlaggingFields = 2
MinimumSlopeForRefinement = 500.0
MetallicityRefinementMinMetallicity = 0.5
MetallicityRefinementMinLevel = 3
#MetallicityRefinementMinDensity =


#
#  set some global parameters
#
PotentialIterations = 10



MaximumGravityRefinementLevel  = 0
MaximumParticleRefinementLevel = 0
SubgridSizeAutoAdjust          = 1
OptimalSubgridsPerProcessor    = 8
MinimumSubgridEdge             = 40
MaximumSubgridSize             = 50000
MinimumEfficiency              = 0.4
NumberOfBufferZones            = 4
LoadBalancing                  = 4 // space-filling curve


RadiativeCoolingModel          = 1
GadgetEquilibriumCooling       = 0
CIECooling                     = 0
H2OpticalDepthApproximation    = 1
ThreeBodyRate                  = 0
H2FormationOnDust              = 0
CloudyElectronFractionFactor   = 0.009154
MetalCooling                   = 1
MetalCoolingTable              = metal_cool_VS2010.dat
HydrogenFractionByMass         = 0.760000
DeuteriumToHydrogenRatio       = 0.000068
SolarMetalFractionByMass       = 0.020410
PhotoelectricHeating           = 1
PhotoelectricHeatingRate       = 1e-27
OutputCoolingTime              = 1
OutputTemperature              = 1
JeansRefinementColdTemperature  = -1.000000
RefineByResistiveLengthSafetyFactor  = 2.000000
MustRefineParticlesRefineToLevel = 0
MustRefineParticlesRefineToLevelAutoAdjust = 0
MustRefineParticlesMinimumMass = 0.000000
ParticleTypeInFile               = 1
WriteGhostZones                  = 0
ReadGhostZones                   = 0
OutputParticleTypeGrouping       = 0
MoveParticlesBetweenSiblings     = 1
ComovingCoordinates                   = 0
StarParticleCreation                  = 0
BigStarFormation                      = 0
BigStarFormationDone                  = 0
BigStarSeparation                     = 0.250000
SimpleQ                               = 1e+50
SimpleRampTime                        = 0.100000
StarFormationOncePerRootGridTimeStep  = 0
StarParticleFeedback                  = 0
NumberOfParticleAttributes            = 3
AddParticleAttributes                 = 0
SinkMergeDistance                     = 10000000272564224.000000
SinkMergeMass                         = 0.100000
StellarWindFeedback                   = 0
StellarWindTurnOnMass                 = 0.100000
MSStellarWindTurnOnMass                 = 10.000000
StarMakerOverDensityThreshold         = 6.46753e10 //1.9454e-22
StarMakerSHDensityThreshold           = 7e-26
StarMakerMassEfficiency               = 0.5
StarMakerMinimumMass                  = 200
StarMakerMinimumDynamicalTime         = 3000000
StarMassEjectionFraction              = 0.8
StarMetalYield                        = 0.02
StarEnergyToThermalFeedback           = 5.93884481e-06
StarEnergyToStellarUV                 = 1e-06
StarEnergyToQuasarUV                  = 1.67e-06
StarFeedbackDistRadius                = 0
StarFeedbackDistCellStep              = 0
StarMakerTypeIaSNe                    = 0
StarMakerPlanetaryNebulae             = 0
MultiMetals                           = 0

CoolingCutOffDensity1      = 0
CoolingCutOffDensity2      = 1e+10
CoolingCutOffTemperature   = 0
CoolingPowerCutOffDensity1 = 0
CoolingPowerCutOffDensity2 = 1e+10
UseGasDrag                       = 0
GasDragCoefficient               = 0.000000
AngularVelocity              = 0.001000
VelocityGradient             = 1.000000
ShearingVelocityDirection    = -1
ShearingBoxProblemType    = 0