diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index efb73ecae8..9c424896a0 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -70,7 +70,7 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
     const CompositionModel<reactingCloudType>& composition =
         td.cloud().composition();
 
-    const scalar TMax = td.cloud().phaseChange().TMax(pc_, this->Tc_);
+    const scalar TMax = td.cloud().phaseChange().TMax(pc_);
     const scalar Tdash = min(T, TMax);
     const scalar Tsdash = min(Ts, TMax);
 
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
index 5d7616b158..a3fe5dff49 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
@@ -247,7 +247,7 @@ Foam::scalar Foam::LiquidEvaporation<CloudType>::TMax(const scalar pIn) const
     scalar T = -GREAT;
     forAll(liquids_, i)
     {
-        T = max(T, liquids_.properties()[i].pv(pIn));
+        T = max(T, liquids_.properties()[i].pvInvert(pIn));
     }
 
     return T;
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
index 468975d179..60dbfcbcce 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
@@ -349,7 +349,7 @@ Foam::scalar Foam::LiquidEvaporationBoil<CloudType>::TMax
     scalar T = -GREAT;
     forAll(liquids_, i)
     {
-        T = max(T, liquids_.properties()[i].pv(pIn));
+        T = max(T, liquids_.properties()[i].pvInvert(pIn));
     }
 
     return T;