# ------------------  INPUTS TO MAIN PROGRAM  -------------------
max_step = 20

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 0 0
geometry.coord_sys   = 1
geometry.prob_lo     = 0.   0.
geometry.prob_hi     = 7.5e8 7.5e8

amr.n_cell           =  256  256

# >>>>>>>>>>>>>  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.add_ext_src   = 0
castro.do_sponge     = 0
castro.ppm_type      = 1
castro.do_grav       = 1
gravity.gravity_type = PoissonGrav
gravity.drdxfac      = 4
gravity.drdxfac      = 1

gravity.rel_tol = 1.e-9

# RELEVANT FOR EOS
castro.small_temp    = 1.e-3
castro.small_dens    = 1.e-6

# TIME STEP CONTROL
castro.cfl            = 0.5     # cfl number for hyperbolic system
castro.init_shrink    = 0.1     # scale back initial timestep
castro.change_max     = 1.05    # scale back initial timestep

# SPONGE
castro.sponge_lower_density = 1.0e-3
castro.sponge_upper_density = 1.0e-3
castro.sponge_timescale = 1.0e-3

# DIAGNOSTICS & VERBOSITY
castro.sum_interval   = 1       # timesteps between computing mass
castro.v              = 0       # 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       = 2         # maximum level number allowed
castro.dual_energy_eta2 = 1.e-4 # this is the default

# Note that if you run with max_level = 3, you need to relax the eta2 value.
#amr.max_level       = 3       # maximum level number allowed
#castro.dual_energy_eta2 = 1.e-3

amr.ref_ratio       = 2 2 2 2 # refinement ratio
amr.regrid_int      = 2 2 2 2 # how often to regrid
amr.blocking_factor = 16      # block factor in grid generation
amr.max_grid_size   = 64
amr.n_error_buf     = 2 2 2 2 # number of buffer cells in error est

amr.refinement_indicators = dengrad

amr.refine.dengrad.gradient = 1.0e7
amr.refine.dengrad.field_name = density
amr.refine.dengrad.max_level = 3

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

# PLOTFILES
amr.plot_file       = plt_2d_  # root name of plotfile
amr.plot_int        = 10000         # number of timesteps between plotfiles

amr.plot_vars        = density xmom ymom rho_E rho_e Temp grav_x grav_y
amr.derive_plot_vars = NONE

# PROBLEM PARAMETERS
problem.rho_0 = 1.e9
problem.r_0   = 6.5e8
problem.p_0   = 1.e15
problem.rho_ambient = 1.0e-5
problem.smooth_delta = 4.e6

problem.center_x = 0.e0
problem.center_y = 0.e0
problem.center_z = 0.e0

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