############################## 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.dt_cutoff = 1.e-8                           # Level 0 timestep below which we halt
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

############################################################################################ 
# 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.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.rotational_dPdt = -0.0                      # Time rate of change of the rotational period
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

############################################################################################
# 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

############################################################################################
# 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 bndy_diag.out primary_diag.out secondary_diag.out extrema_diag.out rotation_diag.out

amr.derive_plot_vars = NONE