# ------------------ INPUTS TO MAIN PROGRAM ------------------- max_step = 250 stop_time = 0.1 # PROBLEM SIZE & GEOMETRY geometry.is_periodic = 0 geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical geometry.prob_lo = 0 geometry.prob_hi = 102.4 amr.n_cell = 1024 # >>>>>>>>>>>>> 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 = 1 castro.diffuse_temp = 1 castro.diffuse_cutoff_density = 1.e-2 castro.ppm_type = 1 # TIME STEP CONTROL castro.cfl = 0.6 # cfl number for hyperbolic system castro.init_shrink = 0.1 # scale back initial timestep castro.change_max = 1.1 # max time step growth # SPONGE castro.sponge_upper_density = 5.0e-8 castro.sponge_lower_density = 1.0e-8 castro.sponge_timescale = 1.0e-6 # DIAGNOSTICS & VERBOSITY castro.sum_interval = 1 # timesteps between computing mass castro.v = 1 # verbosity in Castro.cpp amr.v = 1 # verbosity in Amr.cpp # REFINEMENT / REGRIDDING amr.max_level = 1 # 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 = 64 # 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.checkpoint_files_output = 0 amr.check_file = chk # root name of checkpoint file amr.check_int = 1000 # number of timesteps between checkpoints # PLOTFILES amr.plot_file = plt # root name of plotfile amr.plot_int = 1000 amr.derive_plot_vars = ALL castro.store_omegadot = 1 # problem initialization problem.rho_fuel = 1.0e0 problem.T_fuel = 1.0e0 problem.pert_frac = 0.2e0 problem.pert_delta = 0.05 # refinement amr.refinement_indicators = denerr dengrad amr.refine.denerr.max_level = 5 amr.refine.denerr.value_greater = 1.e-7 amr.refine.denerr.field_name = density amr.refine.dengrad.max_level = 5 amr.refine.dengrad.gradient = 0.01 amr.refine.dengrad.field_name = density # Microphysics eos.eos_assume_neutral = 1 network.rtilde = 40.0 network.nu = 4.0 network.specific_q_burn = 1.24e8 network.T_burn_ref = 1.0 network.rho_burn_ref = 1.0 network.f_act = 1.2 conductivity.const_conductivity = 5.e9