# ------------------  INPUTS TO MAIN PROGRAM  -------------------
max_step = 100000
stop_time =  1.0

# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 0 
geometry.coord_sys   = 2 
geometry.prob_lo     = 0 
geometry.prob_hi     = 5.e8

amr.n_cell           = 640

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

# WHICH PHYSICS
castro.do_hydro      = 1
castro.do_react      = 0
castro.add_ext_src   = 1
castro.do_grav       = 1
gravity.gravity_type = MonopoleGrav
castro.do_sponge     = 0

# THERM
castro.small_temp    = 5.e6

# TIME STEP CONTROL
castro.cfl            = 0.9     # 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_upper_density = 1.e4
castro.sponge_lower_density = 1.e3
castro.sponge_timescale     = 1.e-3

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

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

# PLOTFILES
amr.plot_file       = hse_adjust_plt  # root name of plotfile
amr.plot_int        = 1000            # 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.model_name = "spherical.hse.640"

problem.heating_time = 0.5e0
problem.heating_rad = 0.0e0
problem.heating_peak = 1.e16
problem.heating_sigma = 1.e7
problem.prob_type = 1