diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 05f99990569327a6e7782a5dcd90cfd193217189..4d1b32c3ad6d182bbcbfd9b824b2d7389cc92a3d 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -777,6 +777,9 @@ std::vector<double> const& RichardsMechanicsLocalAssembler<
                 .template value<double>(variables, x_position, t, dt);
         variables[static_cast<int>(MPL::Variable::temperature)] = temperature;
 
+        variables[static_cast<int>(MPL::Variable::porosity)] =
+            _ip_data[ip].porosity;
+
         auto const mu = liquid_phase.property(MPL::PropertyType::viscosity)
                             .template value<double>(variables, x_position, t, dt);
         auto const rho_LR =