# ------------------  INPUTS TO MAIN PROGRAM  -------------------
max_step = 50000
stop_time = 0.06

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic =  1 1
geometry.coord_sys   =  0       # 0 => cart
geometry.prob_lo     =  0    0
geometry.prob_hi     =  1.e8    1.e8
amr.n_cell           = 128  128

# >>>>>>>>>>>>>  BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior           3 = Symmetry
# 1 = Inflow             4 = SlipWall
# 2 = Outflow            5 = NoSlipWall
# >>>>>>>>>>>>>  BC FLAGS <<<<<<<<<<<<<<<<
castro.lo_bc       =  0   0
castro.hi_bc       =  0   0

# WHICH PHYSICS
castro.do_hydro = 1
castro.do_react = 1

# TIME STEP CONTROL

castro.cfl            = 0.5     # cfl number for hyperbolic system
castro.init_shrink    = 0.01    # scale back initial timestep
castro.change_max     = 1.1     # maximum increase in dt over successive steps
castro.fixed_dt = 1.e-4

castro.sdc_order=4
castro.time_integration_method=2
castro.limit_fourth_order=1
castro.use_reconstructed_gamma1=1

# DIAGNOSTICS & VERBOSITY
castro.sum_interval   = 1       # timesteps between computing mass
castro.v              = 1       # verbosity in Castro.cpp
amr.v                 = 1       # verbosity in Amr.cpp
#amr.grid_log         = grdlog  # name of grid logging file

# REFINEMENT / REGRIDDING
amr.max_level       = 0       # maximum level number allowed
amr.ref_ratio       = 2 2 2 2 # refinement ratio
amr.regrid_int      = 2       # how often to regrid
amr.blocking_factor = 8       # block factor in grid generation
amr.max_grid_size   = 256

# CHECKPOINT FILES
amr.checkpoint_files_output = 0
amr.check_file      = react_converge_128_chk   # root name of checkpoint file
amr.check_int       = 100       # number of timesteps between checkpoints

# PLOTFILES
amr.plot_file       = react_converge_128_plt
amr.plot_int        = -1
amr.plot_per        = 0.24
amr.derive_plot_vars=ALL


# Problem initialization

problem.rho0 = 5.e5
problem.T0 = 3.e8
problem.dp_fact = 2.0
problem.L_pert = 2.e7

# Refinement

amr.refinement_indicators = denerr

amr.refine.denerr.max_level = 3
amr.refine.denerr.value_greater = 9.e5
amr.refine.denerr.field_name = density

# Microphysics

network.small_x = 1.e-10

integrator.rtol_spec = 1.e-6
integrator.atol_spec = 1.e-10
integrator.rtol_enuc = 1.e-6
integrator.atol_enuc = 1.e-10