diff --git a/ProcessLib/HeatConduction/HeatConductionProcess.cpp b/ProcessLib/HeatConduction/HeatConductionProcess.cpp
index 7124426f8572959b4c29f2c64f115b487242108d..7affa6cac585719f49a9ae0b96a0930abd9194b8 100644
--- a/ProcessLib/HeatConduction/HeatConductionProcess.cpp
+++ b/ProcessLib/HeatConduction/HeatConductionProcess.cpp
@@ -132,7 +132,7 @@ void HeatConductionProcess::computeSecondaryVariableConcrete(
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &HeatConductionLocalAssemblerInterface::computeSecondaryVariable,
         _local_assemblers, pv.getActiveElementIDs(), getDOFTable(process_id), t,
-        dt, x, x_dot, _coupled_solutions);
+        dt, x, x_dot);
 }
 
 }  // namespace HeatConduction
diff --git a/ProcessLib/HeatTransportBHE/HeatTransportBHEProcess.cpp b/ProcessLib/HeatTransportBHE/HeatTransportBHEProcess.cpp
index a7f40e2c8bc8f1a307ef690af68234c128b69098..96a943db91055b877b343703f8d7b1a7b1c0bc82 100644
--- a/ProcessLib/HeatTransportBHE/HeatTransportBHEProcess.cpp
+++ b/ProcessLib/HeatTransportBHE/HeatTransportBHEProcess.cpp
@@ -196,7 +196,7 @@ void HeatTransportBHEProcess::computeSecondaryVariableConcrete(
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &HeatTransportBHELocalAssemblerInterface::computeSecondaryVariable,
         _local_assemblers, pv.getActiveElementIDs(), getDOFTable(process_id), t,
-        dt, x, x_dot, _coupled_solutions);
+        dt, x, x_dot);
 }
 
 #ifdef OGS_USE_PYTHON
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
index b994ada35a7098faa7fc7d3b937c341f39026e78..d3f93d024b4425d47224c8632f7dec083f7c9a7b 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
+++ b/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
@@ -524,12 +524,16 @@ void HydroMechanicsProcess<DisplacementDim>::computeSecondaryVariableConcrete(
     double const t, double const dt, GlobalVector const& x,
     GlobalVector const& x_dot, const int process_id)
 {
+    if (process_id != 0)
+    {
+        return;
+    }
+
     DBUG("Compute the secondary variables for HydroMechanicsProcess.");
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LocalAssemblerIF::computeSecondaryVariable, _local_assemblers,
-        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot,
-        _coupled_solutions);
+        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot);
 }
 
 template <int DisplacementDim>
diff --git a/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp b/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp
index 26deeeb36693e9bf686b377661a3fb3496343da5..9d07f5adb9fd330f9c0d3943c440c11929e8c4e9 100644
--- a/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp
+++ b/ProcessLib/LIE/SmallDeformation/SmallDeformationProcess.cpp
@@ -537,7 +537,7 @@ void SmallDeformationProcess<DisplacementDim>::computeSecondaryVariableConcrete(
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &SmallDeformationLocalAssemblerInterface::computeSecondaryVariable,
         _local_assemblers, pv.getActiveElementIDs(), getDOFTable(process_id), t,
-        dt, x, x_dot, _coupled_solutions);
+        dt, x, x_dot);
 }
 
 template <int DisplacementDim>
diff --git a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
index ae54d245ff1d829935e1652332735943f28e08a1..51c7b989f77218da344e4ec24c0f811cb69d7d66 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
+++ b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
@@ -120,7 +120,7 @@ void LiquidFlowProcess::computeSecondaryVariableConcrete(
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LiquidFlowLocalAssemblerInterface::computeSecondaryVariable,
         _local_assemblers, pv.getActiveElementIDs(), getDOFTable(process_id), t,
-        dt, x, x_dot, _coupled_solutions);
+        dt, x, x_dot);
 }
 
 Eigen::Vector3d LiquidFlowProcess::getFlux(
diff --git a/ProcessLib/LocalAssemblerInterface.cpp b/ProcessLib/LocalAssemblerInterface.cpp
index 40d0880dc770a66eee9af54a4c9474d9f7a98d86..a1e30460f7265e846fd76ca1f0c4f12c76e51b38 100644
--- a/ProcessLib/LocalAssemblerInterface.cpp
+++ b/ProcessLib/LocalAssemblerInterface.cpp
@@ -72,16 +72,9 @@ void LocalAssemblerInterface::assembleWithJacobianForStaggeredScheme(
 void LocalAssemblerInterface::computeSecondaryVariable(
     std::size_t const mesh_item_id,
     NumLib::LocalToGlobalIndexMap const& dof_table, double const t,
-    double const dt, GlobalVector const& x, GlobalVector const& x_dot,
-    CoupledSolutionsForStaggeredScheme const* coupled_xs)
+    double const dt, GlobalVector const& x, GlobalVector const& x_dot)
 {
     auto const indices = NumLib::getIndices(mesh_item_id, dof_table);
-
-    if (coupled_xs != nullptr)
-    {
-        return;
-    }
-
     auto const local_x = x.get(indices);
     auto const local_x_dot = x_dot.get(indices);
     computeSecondaryVariableConcrete(t, dt, local_x, local_x_dot);
diff --git a/ProcessLib/LocalAssemblerInterface.h b/ProcessLib/LocalAssemblerInterface.h
index ca90d44fbff9749d055ef49132bcf63caa5e64fa..9a42e20b46d663a6ae00f6bf7fe55d81d46b570d 100644
--- a/ProcessLib/LocalAssemblerInterface.h
+++ b/ProcessLib/LocalAssemblerInterface.h
@@ -81,8 +81,7 @@ public:
         std::size_t const mesh_item_id,
         NumLib::LocalToGlobalIndexMap const& dof_table, double const t,
         double const dt, GlobalVector const& local_x,
-        GlobalVector const& local_x_dot,
-        CoupledSolutionsForStaggeredScheme const* coupled_xs);
+        GlobalVector const& local_x_dot);
 
     virtual void preTimestep(std::size_t const mesh_item_id,
                              NumLib::LocalToGlobalIndexMap const& dof_table,
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp b/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
index c7e64883de94864c307c8b0b2aed628e1eb17363..0f117c1f91266abe383ea5f5018aff5afb9818cd 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsProcess.cpp
@@ -586,8 +586,7 @@ void RichardsMechanicsProcess<DisplacementDim>::
 
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LocalAssemblerIF::computeSecondaryVariable, _local_assemblers,
-        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot,
-        _coupled_solutions);
+        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot);
 }
 
 template <int DisplacementDim>
diff --git a/ProcessLib/SmallDeformation/SmallDeformationProcess.cpp b/ProcessLib/SmallDeformation/SmallDeformationProcess.cpp
index b2b04be822b964b98d532c8bf03507230029b523..f66e58fbf6779b8fc214ffa93bc50c46d542a28b 100644
--- a/ProcessLib/SmallDeformation/SmallDeformationProcess.cpp
+++ b/ProcessLib/SmallDeformation/SmallDeformationProcess.cpp
@@ -270,8 +270,7 @@ void SmallDeformationProcess<DisplacementDim>::computeSecondaryVariableConcrete(
     ProcessLib::ProcessVariable const& pv = getProcessVariables(process_id)[0];
     GlobalExecutor::executeSelectedMemberOnDereferenced(
         &LocalAssemblerInterface::computeSecondaryVariable, _local_assemblers,
-        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot,
-        _coupled_solutions);
+        pv.getActiveElementIDs(), getDOFTable(process_id), t, dt, x, x_dot);
 }
 template class SmallDeformationProcess<2>;
 template class SmallDeformationProcess<3>;
diff --git a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsProcess.cpp b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsProcess.cpp
index 33d4a079f9c92573a85f5de92bb9e94db1ff908f..e99e71ba8d8c06cb9be97f6e4e42115dfc432f9c 100644
--- a/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsProcess.cpp
+++ b/ProcessLib/ThermoHydroMechanics/ThermoHydroMechanicsProcess.cpp
@@ -385,7 +385,7 @@ void ThermoHydroMechanicsProcess<DisplacementDim>::
     DBUG("Compute the secondary variables for ThermoHydroMechanicsProcess.");
     GlobalExecutor::executeMemberOnDereferenced(
         &LocalAssemblerInterface::computeSecondaryVariable, _local_assemblers,
-        getDOFTable(process_id), t, dt, x, x_dot, _coupled_solutions);
+        getDOFTable(process_id), t, dt, x, x_dot);
 }
 
 template <int DisplacementDim>