max_step = 100000000
stop_time = 4.0

geometry.is_periodic = 0
geometry.coord_sys   = 0
geometry.prob_lo     = -1.6384e9
geometry.prob_hi     =  1.6384e9
amr.n_cell           = 256

castro.use_retry = 1
castro.max_subcycles = 128

castro.add_ext_src = 1

castro.fill_ambient_bc = 1

castro.lo_bc = 2
castro.hi_bc = 2

castro.do_hydro = 1
castro.do_react = 1

castro.cfl            = 0.5
castro.init_shrink    = 0.1
castro.change_max     = 1.25
castro.dtnuc_e        = 1.e200
castro.dtnuc_X        = 1.e200

castro.small_temp     = 1.e5

castro.sum_interval   = 0
castro.v              = 1
amr.v                 = 1

amr.max_level       = 1
amr.ref_ratio       = 4
amr.regrid_int      = 1
amr.blocking_factor = 64
amr.max_grid_size   = 64

amr.n_error_buf     = 2

amr.refinement_indicators = tempgrad_rel dengrad_rel dxnucerr

amr.refine.tempgrad_rel.max_level = 0
amr.refine.tempgrad_rel.relative_gradient = 0.5
amr.refine.tempgrad_rel.field_name = Temp

amr.refine.dengrad_rel.max_level = 0
amr.refine.dengrad_rel.relative_gradient = 0.5
amr.refine.dengrad_rel.field_name = density

amr.refine.dxnucerr.max_level = 15
amr.refine.dxnucerr.value_greater = 1.e-16
amr.refine.dxnucerr.field_name = t_sound_t_enuc

castro.use_post_step_regrid = 1

amr.check_file       = det_x_chk
amr.check_per        = 0.1
amr.plot_file        = det_x_plt
amr.plot_per         = 0.1
amr.plot_vars        = ALL
amr.derive_plot_vars = ALL
amr.small_plot_per   = 0.01
amr.small_plot_vars  = density Temp rho_e rho_c12 rho_o16 rho_si28 rho_ni56 enuc
amr.derive_small_plot_vars = pressure soundspeed x_velocity t_sound_t_enuc

amr.plotfile_on_restart = 0
amr.checkpoint_on_restart = 0
amr.write_plotfile_with_checkpoint = 0
castro.plot_per_is_exact = 1
castro.small_plot_per_is_exact = 1
castro.output_at_completion = 1

castro.stopping_criterion_field = Temp
castro.stopping_criterion_value = 4.e9

# problem initialization

problem.T_l = 1.e7
problem.T_r = 1.e7

problem.dens = 5.e6
problem.cfrac = 0.5
problem.ofrac = 0.45

problem.smallx = 1.e-12

problem.idir = 1

problem.w_T = 1.e-12
problem.center_T = 0.5

problem.vel = 2.0e8

problem.grav_acceleration = 1.1e8

# Microphysics

network.small_x = 1.e-12

eos.use_eos_coulomb = 1
eos.eos_input_is_constant = 1

integrator.call_eos_in_rhs = 1

integrator.rtol_spec = 1.e-6
integrator.atol_spec = 1.e-6

integrator.rtol_enuc = 1.e-6
integrator.atol_enuc = 1.e-6

integrator.renormalize_abundances = 1

integrator.MAX_TEMP = 1.e10

network.use_tables = 1