inputs_3d - the "baseline" run
  -CASTRO initial model
  -burning_cutoff_density_hi = 5.e6 (disable burning for rho > 5.e6, which is r~2.2e8)
  -anelastic_cutoff_density = 1.e3
  -no velocity perturbation
  -HOEXTRAP boundary conditions
  -evolve_base_state=false
  -average_base_state=F (original temporal integrator)

inputs_3d_perturb
  -maestro.perturb_model = true (velocity perturbations in Si shell)

inputs_3d_MESA
  -maestro.model_file = "initial_models/15m_500_sec.hse.6400" (MESA initial model)

inputs_3d_cutoff
  -maestro.anelastic_cutoff_density = 1.e5

inputs_3d_heating
  -maestro.do_burning = false
  -maestro.do_heating = true
  -use_analytic_heating = T
  
inputs_3d_heatinner
  -maestro.do_burning = false
  -maestro.do_heating = true
  -use_analytic_heating = T
  -maestro.heating_cutoff_density_lo = 5.e5  # apply heating in the inner part of the star only

inputs_3d_noburn
  -maestro.do_burning = false

inputs_3d_burnrange
  -burning_cutoff_density_lo = 5.e5 (restrict burning to the Si shell)

inputs_3d_burntocore
  -maestro.init_shrink = 0.00001
  -no burning_cutoff_density_hi
