# ------------------ INPUTS TO MAIN PROGRAM ------------------- max_step = 9900000 stop_time = 6.e-5 # PROBLEM SIZE & GEOMETRY geometry.is_periodic = 0 0 geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical geometry.prob_lo = 0 0 geometry.prob_hi = 2.56e3 amr.n_cell = 128 # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< # 0 = Interior 3 = Symmetry # 1 = Inflow 4 = SlipWall # 2 = Outflow 5 = NoSlipWall # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< castro.lo_bc = 1 castro.hi_bc = 2 castro.xl_ext_bc_type = 1 castro.hse_interp_temp = 1 castro.hse_reflect_vels = 1 castro.fill_ambient_bc = 1 castro.ambient_fill_dir = 0 castro.ambient_outflow_vel = 1 # WHICH PHYSICS castro.do_hydro = 1 castro.do_react = 0 castro.do_rotation = 0 castro.do_grav = 1 castro.do_sponge = 0 castro.small_temp = 1.e6 castro.small_dens = 1.e-5 castro.ppm_type = 1 castro.grav_source_type = 2 gravity.gravity_type = ConstantGrav gravity.const_grav = -1.5e14 castro.diffuse_temp = 0 # TIME STEP CONTROL castro.cfl = 0.8 # cfl number for hyperbolic system castro.init_shrink = 1.0 # scale back initial timestep castro.change_max = 1.1 # max time step growth castro.use_retry = 1 castro.max_subcycles = 16 castro.fixed_dt = 3.75e-8 # 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 = 0 # maximum level number allowed amr.ref_ratio = 4 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 = 128 amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est # CHECKPOINT FILES amr.check_file = flame_wave_chk # root name of checkpoint file amr.check_int = 250000 # number of timesteps between checkpoints # PLOTFILES amr.plot_file = flame_wave_128_plt # root name of plotfile amr.plot_per = 1.e-5 # number of seconds between plotfiles amr.derive_plot_vars = ALL # PROBLEM PARAMETERS problem.dens_base = 3.43e6 problem.T_star = 1.e8 problem.T_hi = 2.e8 problem.T_lo = 8.e6 problem.H_star = 1700.e0 problem.atm_delta = 50.0 problem.fuel1_name = "helium-4" problem.fuel1_frac = 1.0e0 problem.ash1_name = "nickel-56" problem.ash1_frac = 1.0e0 problem.low_density_cutoff = 1.e-4 problem.cutoff_density = 2.5e4 problem.X_min = 1.e-2