# ------------------  INPUTS TO MAIN PROGRAM  -------------------
max_step = 10000
stop_time =  0.15

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 0
geometry.coord_sys   = 0

geometry.prob_lo     = 0.0
geometry.prob_hi     = 1.0

amr.n_cell           = 128
amr.max_level        = 0

#amr.n_cell           = 2048
#amr.max_level        = 0

#amr.n_cell           = 512
#amr.max_level        = 2
#amr.ref_ratio         = 2 2

#amr.n_cell            = 128
#amr.max_level         = 2
#amr.ref_ratio         = 4 4

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

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

# TIME STEP CONTROL

castro.cfl            = 0.9     # cfl number for hyperbolic system
castro.init_shrink    = 1.0     # scale back initial timestep
castro.change_max     = 1.05    # scale back initial timestep

# 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.regrid_int      = 2       # how often to regrid
amr.blocking_factor = 16       # block factor in grid generation
amr.max_grid_size   = 2048
amr.n_error_buf     = 2 2 2 2 # number of buffer cells in error est

amr.refinement_indicators = denerr dengrad presserr pressgrad velgrad_x velgrad_y velgrad_z

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

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

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

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

amr.refine.velgrad_x.gradient = 0.01
amr.refine.velgrad_x.field_name = x_velocity
amr.refine.velgrad_x.max_level = 5

amr.refine.velgrad_y.gradient = 0.01
amr.refine.velgrad_y.field_name = y_velocity
amr.refine.velgrad_y.max_level = 5

amr.refine.velgrad_z.gradient = 0.01
amr.refine.velgrad_z.field_name = z_velocity
amr.refine.velgrad_z.max_level = 5

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

# PLOTFILES
amr.plot_files_output = 0
amr.plot_file       = test2_1d_plt     # root name of plotfile
amr.plot_int        = 10000           # number of timesteps between plotfiles
amr.derive_plot_vars= density xmom ymom zmom eden Temp pressure # these variables appear in the plotfile

# PROBLEM PARAMETERS
problem.p_l = 0.4
problem.u_l = -2.0
problem.rho_l = 1.0

problem.p_r = 0.4
problem.u_r = 2.0
problem.rho_r = 1.0

problem.idir = 1
problem.frac = 0.5

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