diff --git a/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.cpp b/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.cpp
index 16d0507d1be48a1adb81bac8cd71eb0ab46758c5..34e56be520db6464e8777e7aaa77c8eb229d67a7 100644
--- a/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.cpp
@@ -18,6 +18,8 @@ namespace ProcessLib::ThermoRichardsMechanics
 template <int DisplacementDim>
 void FluidThermalExpansionModel<DisplacementDim>::eval(
     SpaceTimeData const& x_t, MediaData const& media_data,
+    CapillaryPressureData<DisplacementDim> const& p_cap_data,
+    TemperatureData<DisplacementDim> const& T_data,
     SolidThermalExpansionData<DisplacementDim> const& s_therm_exp_data,
     PorosityData const& poro_data, LiquidDensityData const& rho_L_data,
     BiotData const& biot_data, FluidThermalExpansionData& out) const
@@ -28,6 +30,9 @@ void FluidThermalExpansionModel<DisplacementDim>::eval(
 
     namespace MPL = MaterialPropertyLib;
     MPL::VariableArray variables;
+    variables[static_cast<int>(MPL::Variable::phase_pressure)] =
+        -p_cap_data.p_cap;
+    variables[static_cast<int>(MPL::Variable::temperature)] = -T_data.T;
 
     double const phi = poro_data.phi;
     double const alpha = biot_data.alpha;
diff --git a/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.h b/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.h
index b8bfe7525688b3e9ea5013d00a07d8aad752bb63..d960c79e41c7aab1195eacdd63c407703c8582b2 100644
--- a/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.h
+++ b/ProcessLib/ThermoRichardsMechanics/Constitutive/FluidThermalExpansion.h
@@ -27,6 +27,8 @@ struct FluidThermalExpansionModel
 {
     void eval(
         SpaceTimeData const& x_t, MediaData const& media_data,
+        CapillaryPressureData<DisplacementDim> const& p_cap_data,
+        TemperatureData<DisplacementDim> const& T_data,
         SolidThermalExpansionData<DisplacementDim> const& s_therm_exp_data,
         PorosityData const& poro_data, LiquidDensityData const& rho_L_data,
         BiotData const& biot_data, FluidThermalExpansionData& out) const;
diff --git a/ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.cpp b/ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.cpp
index 8d39624c77a9ceb2b298d75428062eb68f2e6c9c..c9bdb5cc718855f9e4720d04de7f8b9feb88280d 100644
--- a/ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.cpp
+++ b/ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.cpp
@@ -144,8 +144,9 @@ void ConstitutiveSetting<DisplacementDim>::eval(
                                       dS_L_data, poro_data, p_cap_data, T_data,
                                       vap_data);
 
-    models.f_therm_exp_model.eval(x_t, media_data, s_therm_exp_data, poro_data,
-                                  rho_L_data, biot_data, f_therm_exp_data);
+    models.f_therm_exp_model.eval(x_t, media_data, p_cap_data, T_data,
+                                  s_therm_exp_data, poro_data, rho_L_data,
+                                  biot_data, f_therm_exp_data);
 
     models.storage_model.eval(x_t, biot_data, poro_data, rho_L_data, S_L_data,
                               dS_L_data, prev_state.S_L_data, p_cap_data,