Minor fix for initial dt. We need to reset dt at the end of (#83)

InitialIterations in order to get the specifed dt. Otherwise, we can end up with dt = initial_dt*change_max.
AMReX-Fluids · Nov 2, 2023 · 58de1d9 · 58de1d9
1 parent afa2a5c
commit 58de1d9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/setup/init.cpp b/src/setup/init.cpp
@@ -311,6 +311,10 @@ void incflo::InitialIterations ()
         copy_from_old_to_new_density();
         copy_from_old_to_new_tracer();
     }
+
+    // Reset dt to get initial step as specified, otherwise we can see increase to dt
+    m_prev_dt = Real(-1.0);
+    m_dt = Real(-1.0);
 }
 
 // Project velocity field to make sure initial velocity is divergence-free