diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index 2f2fe779e7af872abc27b6367bddc02da4e57cbc..7899b27f5bc2293c9d4e70e6e33b27a08a2f587d 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -1076,8 +1076,6 @@ void TH2MLocalAssembler<
         auto const s_L = current_state.S_L_data.S_L;
         auto const s_L_dot = (s_L - prev_state.S_L_data->S_L) / dt;
 
-        auto& alpha_B = ip_cv.biot_data();
-
         auto const& b = this->process_data_.specific_body_force;
 
         // ---------------------------------------------------------------------
@@ -1172,7 +1170,7 @@ void TH2MLocalAssembler<
         //  - displacement equation
         // ---------------------------------------------------------------------
 
-        KUpG.noalias() -= (BuT * alpha_B * m * Np) * w;
+        KUpG.noalias() -= (BuT * ip_cv.biot_data() * m * Np) * w;
 
         KUpC.noalias() +=
             BuT * Invariants::identity2 * Np * (ip_cv.fu_2_KupC.m * w);
@@ -1402,8 +1400,6 @@ void TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
         auto const& s_L = current_state.S_L_data.S_L;
         auto const s_L_dot = (s_L - prev_state.S_L_data->S_L) / dt;
 
-        auto const alpha_B = ip_cv.biot_data();
-
         auto const& b = this->process_data_.specific_body_force;
 
         // ---------------------------------------------------------------------
@@ -1747,7 +1743,7 @@ void TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
         //  - displacement equation
         // ---------------------------------------------------------------------
 
-        KUpG.noalias() -= (BuT * alpha_B * m * Np) * w;
+        KUpG.noalias() -= (BuT * ip_cv.biot_data() * m * Np) * w;
 
         // dfU_2/dp_GR = dKUpG/dp_GR * p_GR + KUpG. The former is zero, the
         // latter is handled below.