############################## CASTRO INPUTS ############################################### ############################################################################################ # Problem setup ############################################################################################ amr.probin_file = probin_test_wdmerger_3D # Name of the probin file max_step = 5 # Maximum coarse timestep geometry.is_periodic = 0 0 0 # Non-periodic boundary conditions geometry.coord_sys = 0 # Cartesian coordinate system geometry.prob_lo = -5.12e9 -5.12e9 -5.12e9 # Lower boundary limits in physical space geometry.prob_hi = 5.12e9 5.12e9 5.12e9 # Upper boundary limits in physical space castro.center = 0.0e0 0.0e0 0.0e0 # System center of mass castro.cfl = 0.5 # CFL number for hyperbolic system castro.init_shrink = 0.1 # Scale back initial timestep by this factor castro.change_max = 1.1 # Factor by which dt is allowed to change each timestep amr.subcycling_mode = None # Disable subcycling so the relaxation can proceed ############################################################################################ # Boundary conditions # 0 = Interior 3 = Symmetry # 1 = Inflow 4 = SlipWall # 2 = Outflow 5 = NoSlipWall ############################################################################################ castro.lo_bc = 2 2 2 # Boundary conditions on lo x, y, and z edges castro.hi_bc = 2 2 2 # Boundary conditions on hi x, y, and z edges castro.fill_ambient_bc = 1 castro.clamp_ambient_temp = 1 ############################################################################################ # Resolution, gridding and AMR ############################################################################################ amr.n_cell = 96 96 96 # Number of cells on the coarse grid amr.max_level = 1 # Maximum level number allowed amr.ref_ratio = 2 amr.refinement_indicators = temperature amr.refine.temperature.max_level = 4 amr.refine.temperature.value_greater = 2.e8 amr.refine.temperature.field_name = Temp amr.max_grid_size = 32 # Maximum grid size at each level amr.blocking_factor = 16 # Grid sizes must be a multiple of this amr.grid_eff = 0.9 # What constitutes an efficient grid ############################################################################################ # Physics to include ############################################################################################ castro.do_hydro = 1 # Whether or not to do hydrodynamics castro.do_grav = 1 # Whether or not to do gravity castro.do_react = 0 # Whether or not to do reactions castro.do_sponge = 1 # Whether or not to apply the sponge castro.add_ext_src = 1 # Whether or not to apply external source terms castro.ext_src_implicit = 1 castro.do_rotation = 0 # Whether or not to include the rotation source term castro.rotational_period = 100.0 # Rotational period of the rotating reference frame castro.implicit_rotation_update = 1 # Implicit rotation coupling ############################################################################################ # PPM options ############################################################################################ castro.ppm_type = 1 # Piecewise parabolic with the original limiters (0 is piecewise linear; 2 is new limiters) castro.ppm_temp_fix = 0 # Use the EOS in calculation of the edge states going into the Riemann solver castro.grav_source_type = 4 # How to include the gravity source term in the hydro equations castro.rot_source_type = 4 # How to include the rotation source term in the hydro equations # Explicitly limit fluxes to avoid hitting a negative density castro.limit_fluxes_on_small_dens = 0 ############################################################################################ # Thermodynamics ############################################################################################ castro.small_temp = 1.e5 # Minimum allowable temperature (K) castro.small_dens = 1.e-5 # Minimum allowable density (g / cm**3) castro.dual_energy_eta1 = 1.0e-3 # Threshold for when to use the internal energy in calculating pressure castro.dual_energy_eta2 = 1.0e-1 # Threshold for when to use (E - K) in updating internal energy eos.use_eos_coulomb = 0 eos.eos_input_is_constant = 1 castro.ambient_temp = 1.0e7 castro.ambient_density = 1.0e-4 ############################################################################################ # Gravity ############################################################################################ gravity.gravity_type = PoissonGrav # Full self-gravity with the Poisson equation gravity.max_multipole_order = 6 # Multipole expansion includes terms up to r**(-max_multipole_order) gravity.rel_tol = 1.e-10 # Relative tolerance for multigrid solver gravity.no_sync = 1 # Turn off sync solve for gravity after refluxing ############################################################################################ # Sponge ############################################################################################ castro.sponge_lower_radius = 3.840e9 castro.sponge_upper_radius = 4.352e9 castro.sponge_timescale = 0.01e0 ############################################################################################ # Diagnostics and I/O ############################################################################################ amr.plot_files_output = 1 # Whether or not to output plotfiles amr.checkpoint_files_output = 1 # Whether or not to output checkpoints amr.check_file = chk # Root name of checkpoint file amr.check_int = 10 # Number of timesteps between checkpoints amr.plot_file = plt # Root name of plot file amr.plot_int = 10 # Number of timesteps between plotfiles amr.v = 1 # Control verbosity in Amr.cpp castro.v = 1 # Control verbosity in Castro.cpp gravity.v = 1 # Control verbosity in Gravity.cpp mg.v = 2 # Control verbosity in the multigrid solver amr.data_log = grid_diag.out star_diag.out species_diag.out amr_diag.out primary_diag.out secondary_diag.out extrema_diag.out rotation_diag.out amr.derive_plot_vars = NONE