diff --git a/ProcessLib/TimeLoop.cpp b/ProcessLib/TimeLoop.cpp
index 048e5d88830a1238bd3d44403b0fbccbf3091038..284c8e17a0d4a1e5c78e7199180e66c371b83071 100644
--- a/ProcessLib/TimeLoop.cpp
+++ b/ProcessLib/TimeLoop.cpp
@@ -550,6 +550,8 @@ void TimeLoop::initialize()
         setCoupledSolutions();
     }
 
+    updateDeactivatedSubdomains(_per_process_data, _start_time);
+
     // Output initial conditions
     {
         const bool output_initial_condition = true;
@@ -564,8 +566,6 @@ void TimeLoop::initialize()
         computeTimeStepping(0.0, _current_time, _accepted_steps,
                             _rejected_steps, time_step_constraints);
 
-    updateDeactivatedSubdomains(_per_process_data, _start_time);
-
     calculateNonEquilibriumInitialResiduum(
         _per_process_data, _process_solutions, _process_solutions_prev);
 }