diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index 9cd09120aa1df03e546449f36369170c151e26b8..686f291c91d17a86ee350fdd6e39ad60d5c330f4 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -572,7 +572,7 @@ void TH2MLocalAssembler<
             k_over_mu_G * ip.rhoGR * b - k_over_mu_G * gradpGR;
 
         GlobalDimVectorType const w_LS = k_over_mu_L * gradpCap +
-                                         k_over_mu_L * ip.rhoGR * b -
+                                         k_over_mu_L * ip.rhoLR * b -
                                          k_over_mu_L * gradpGR;
 
         // ---------------------------------------------------------------------