Cleanup

Docs/sphinx_documentation/source/InputsTimeStepping.rst

@@ -28,6 +28,10 @@ The inputs below must be preceded by "incflo."
 +----------------------+-----------------------------------------------------------------------+-------------+--------------+
 | cfl                  | CFL constraint (dt < cfl * dx / u) if fixed_dt not > 0                |    Real     |   0.5        |
 +----------------------+-----------------------------------------------------------------------+-------------+--------------+
+| dt_change_max        | maximum percentage dt can increase                                    |    Real     |   1.1        |
++----------------------+-----------------------------------------------------------------------+-------------+--------------+
+| init_shrink          | reduce the initial time step by this safety factor                    |    Real     |   0.1        |
++----------------------+-----------------------------------------------------------------------+-------------+--------------+
 
 Setting the Time Step
 ---------------------

src/incflo_apply_predictor.cpp

@@ -63,7 +63,7 @@ using namespace amrex;
 //
 //     vel = u** - dt * grad p / rho^nph
 //
-// It is assumed that the ghost cels of the old data have been filled and
+// It is assumed that the ghost cells of the old data have been filled and
 // the old and new data are the same in valid region.
 //
 void incflo::ApplyPredictor (bool incremental_projection)
@@ -139,6 +139,7 @@ void incflo::ApplyPredictor (bool incremental_projection)
 
     // *************************************************************************************
     // Compute explicit viscous term
+    // Note that for !advect_momentum, this actually computes divtau / rho
     // *************************************************************************************
     if (need_divtau() || use_tensor_correction )
     {
@@ -160,7 +161,7 @@ void incflo::ApplyPredictor (bool incremental_projection)
     // *************************************************************************************
     bool include_pressure_gradient = !(m_use_mac_phi_in_godunov);
     compute_vel_forces(GetVecOfPtrs(vel_forces), get_velocity_old_const(),
-                       get_density_old_const(), get_tracer_old_const(), get_tracer_new_const(),
+                       get_density_old_const(), get_tracer_old_const(), get_tracer_old_const(),
                        include_pressure_gradient);
     compute_MAC_projected_velocities(get_velocity_old_const(), get_density_old_const(),
                                      AMREX_D_DECL(GetVecOfPtrs(u_mac), GetVecOfPtrs(v_mac),

src/incflo_compute_dt.cpp

@@ -229,7 +229,7 @@ void incflo::ComputeDt (int initialization, bool explicit_diffusion)
         dt_new = Real(0.5) * m_dt;
     }
 
-    // Don't let the timestep grow by more than 10% per step
+    // Don't let the timestep grow by more than m_dt_change_max per step
     // unless the previous time step was unduly shrunk to match m_plot_per_exact
     Real allowed_change_factor = m_dt_change_max;
     if( (m_dt > Real(0.0)) && !(m_plot_per_exact > 0 && m_last_plt == m_nstep && m_nstep > 0) )