diff --git a/ProcessLib/HT/CreateHTProcess.cpp b/ProcessLib/HT/CreateHTProcess.cpp
index 7b412b43550a9c9c52d73e487407734e1c169086..192a112b5d1357473a9a27666da5a4cb009195c4 100644
--- a/ProcessLib/HT/CreateHTProcess.cpp
+++ b/ProcessLib/HT/CreateHTProcess.cpp
@@ -31,17 +31,22 @@ void checkMPLProperties(
     MaterialPropertyLib::MaterialSpatialDistributionMap const& media_map)
 {
     std::array const required_property_medium = {
+        MaterialPropertyLib::PropertyType::permeability,
         MaterialPropertyLib::PropertyType::porosity,
-        MaterialPropertyLib::PropertyType::permeability};
+        MaterialPropertyLib::PropertyType::thermal_conductivity,
+        MaterialPropertyLib::PropertyType::thermal_longitudinal_dispersivity,
+        MaterialPropertyLib::PropertyType::thermal_transversal_dispersivity};
 
     std::array const required_property_liquid_phase = {
         MaterialPropertyLib::PropertyType::viscosity,
         MaterialPropertyLib::PropertyType::density,
-        MaterialPropertyLib::PropertyType::specific_heat_capacity};
+        MaterialPropertyLib::PropertyType::specific_heat_capacity,
+        MaterialPropertyLib::PropertyType::thermal_conductivity};
 
     std::array const required_property_solid_phase = {
         MaterialPropertyLib::PropertyType::specific_heat_capacity,
         MaterialPropertyLib::PropertyType::density,
+        MaterialPropertyLib::PropertyType::thermal_conductivity,
         MaterialPropertyLib::PropertyType::storage};
 
     MaterialPropertyLib::checkMaterialSpatialDistributionMap(
diff --git a/ProcessLib/HT/MonolithicHTFEM.h b/ProcessLib/HT/MonolithicHTFEM.h
index d4a824c797900a41ccd5638b45f2ed3d7c9c3279..cb4ddee4a9e4d9c4f8da86a5416ca9e0400cff1b 100644
--- a/ProcessLib/HT/MonolithicHTFEM.h
+++ b/ProcessLib/HT/MonolithicHTFEM.h
@@ -134,7 +134,8 @@ public:
             vars[static_cast<int>(
                 MaterialPropertyLib::Variable::phase_pressure)] = p_int_pt;
 
-            vars[static_cast<int>(MaterialPropertyLib::Variable::liquid_saturation)] = 1.0;
+            vars[static_cast<int>(
+                MaterialPropertyLib::Variable::liquid_saturation)] = 1.0;
             // \todo the argument to getValue() has to be changed for non
             // constant storage model
             auto const specific_storage =
@@ -179,8 +180,8 @@ public:
             // matrix assembly
             GlobalDimMatrixType const thermal_conductivity_dispersivity =
                 this->getThermalConductivityDispersivity(
-                    vars, fluid_density, specific_heat_capacity_fluid,
-                    velocity, I, pos, t, dt);
+                    vars, fluid_density, specific_heat_capacity_fluid, velocity,
+                    I, pos, t, dt);
             KTT.noalias() +=
                 (dNdx.transpose() * thermal_conductivity_dispersivity * dNdx +
                  N.transpose() * velocity.transpose() * dNdx * fluid_density *