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

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 0    0
geometry.coord_sys   = 1       # 1 => RZ

# since we are modeling only one quadrant, we can make the domain
# smaller -- but make sure that the resolution agrees with inputs.1d.sph
# and inputs.3d.sph
geometry.prob_lo     = 0    0
geometry.prob_hi     = 0.5  0.5
amr.n_cell           = 16   16

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

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

castro.ppm_temp_fix = 0

castro.use_flattening = 1

castro.riemann_solver = 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

# 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       = 4       # 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   = 32

amr.refinement_indicators = denerr dengrad presserr pressgrad

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

amr.refine.dengrad.max_level = 3
amr.refine.dengrad.gradient = 0.01
amr.refine.dengrad.field_name = density

amr.refine.presserr.max_level = 3
amr.refine.presserr.value_greater = 3
amr.refine.presserr.field_name = pressure

amr.refine.pressgrad.max_level = 3
amr.refine.pressgrad.gradient = 0.01
amr.refine.pressgrad.field_name = pressure

# CHECKPOINT FILES
amr.check_file      = sedov_2d_sph_in_cyl_chk     # root name of checkpoint file
amr.check_int       = 1000     # number of timesteps between checkpoints

# PLOTFILES
amr.plot_file       = sedov_2d_sph_in_cyl_plt
amr.plot_int        = 1000
amr.derive_plot_vars=ALL

# PROBLEM PARAMETERS
problem.r_init = 0.01
problem.p_ambient = 1.e-5
problem.exp_energy = 1.0
problem.dens_ambient = 1.0
problem.nsub = 10

# EOS
eos.eos_assume_neutral = 1
eos.eos_gamma = 1.4