# ------------------ INPUTS TO MAIN PROGRAM ------------------- max_step = 50 stop_time = 0.2 # PROBLEM SIZE & GEOMETRY geometry.is_periodic = 0 0 0 geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical geometry.prob_lo = 0 0 0 geometry.prob_hi = 4.e4 2500 2500 amr.n_cell = 128 8 8 # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< # 0 = Interior 3 = Symmetry # 1 = Inflow 4 = SlipWall # 2 = Outflow 5 = NoSlipWall # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< castro.lo_bc = 3 4 4 castro.hi_bc = 2 4 4 # WHICH PHYSICS castro.do_hydro = 1 castro.do_react = 1 castro.ppm_type = 1 castro.disable_shock_burning = 1 # TIME STEP CONTROL castro.cfl = 0.25 # cfl number for hyperbolic system castro.init_shrink = 0.1 # 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.max_level = 2 # 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 = 64 amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est # CHECKPOINT FILES amr.check_file = det_x_chk # root name of checkpoint file amr.check_int = 100 # number of timesteps between checkpoints # PLOTFILES amr.plot_file = det_x_plt # root name of plotfile amr.plot_int = 100 amr.derive_plot_vars = density xmom ymom zmom eden Temp pressure # these variables appear in the plotfile # problem initialization problem.T_l = 4.e9 problem.T_r = 5.e7 problem.dens = 2.e8 problem.cfrac = 0.0 problem.smallx = 1.e-10 problem.idir = 1 problem.w_T = 5.e-4 problem.center_T = 0.3 # refinement amr.refinement_indicators = temperr tempgrad amr.refine.temperr.max_level = 5 amr.refine.temperr.value_greater = 4.e9 amr.refine.temperr.field_name = Temp amr.refine.tempgrad.max_level = 5 amr.refine.tempgrad.gradient = 1.e8 amr.refine.tempgrad.field_name = Temp # Microphysics integrator.call_eos_in_rhs = 1