diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index 620a870625dc9337767709c897284f036f08b5c0..23d6388e81a7dc2c8694aab9b7973dfe9e95f42f 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -597,8 +597,10 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
 #endif
             ;
 
-        ip_cv.dadvection_C_dp_GR =
-            c.drho_C_GR_dp_GR * k_over_mu_G + c.drho_C_LR_dp_GR * k_over_mu_L;
+        ip_cv.dadvection_C_dp_GR = c.drho_C_GR_dp_GR * k_over_mu_G
+                                   // + rhoCGR * (dk_over_mu_G_dp_GR = 0)
+                                   // + rhoCLR * (dk_over_mu_L_dp_GR = 0)
+                                   + c.drho_C_LR_dp_GR * k_over_mu_L;
 
         ip_cv.dfC_4_LCpG_dT =
             c.drho_C_GR_dT * k_over_mu_G + c.drho_C_LR_dT * k_over_mu_L