diff --git a/applications/utilities/thermophysical/adiabaticFlameT/adiabaticFlameT.C b/applications/utilities/thermophysical/adiabaticFlameT/adiabaticFlameT.C
index e00da5e9e5..58de0b6f77 100644
--- a/applications/utilities/thermophysical/adiabaticFlameT/adiabaticFlameT.C
+++ b/applications/utilities/thermophysical/adiabaticFlameT/adiabaticFlameT.C
@@ -69,6 +69,7 @@ int main(int argc, char *argv[])
     dictionary control(controlFile);
 
 
+    scalar P(readScalar(control.lookup("P")));
     scalar T0(readScalar(control.lookup("T0")));
     const word fuelName(control.lookup("fuel"));
     scalar n(readScalar(control.lookup("n")));
@@ -179,7 +180,7 @@ int main(int argc, char *argv[])
 
         Info<< "products " << (1/products.nMoles())*products << ';' << endl;
 
-        scalar Tad = products.THa(reactants.Ha(T0), 1000.0);
+        scalar Tad = products.THa(reactants.Ha(P, T0), P, 1000.0);
         Info<< "Tad = " << Tad << nl << endl;
     }
 
diff --git a/applications/utilities/thermophysical/equilibriumCO/equilibriumCO.C b/applications/utilities/thermophysical/equilibriumCO/equilibriumCO.C
index 00391cbd68..dde7e8678a 100644
--- a/applications/utilities/thermophysical/equilibriumCO/equilibriumCO.C
+++ b/applications/utilities/thermophysical/equilibriumCO/equilibriumCO.C
@@ -50,9 +50,8 @@ typedef specieThermo<janafThermo<perfectGas>, absoluteEnthalpy> thermo;
 
 int main(int argc, char *argv[])
 {
-
-#   include "setRootCase.H"
-#   include "createTime.H"
+    #include "setRootCase.H"
+    #include "createTime.H"
 
     Info<< nl << "Reading Burcat data IOdictionary" << endl;
 
@@ -71,6 +70,7 @@ int main(int argc, char *argv[])
 
 
 
+    scalar P = 1e5;
     scalar T = 3000.0;
 
     SLPtrList<thermo> EQreactions;
@@ -121,7 +121,7 @@ int main(int argc, char *argv[])
 
     forAllConstIter(SLPtrList<thermo>, EQreactions, iter)
     {
-        Info<< "Kc(EQreactions) = " << iter().Kc(T) << endl;
+        Info<< "Kc(EQreactions) = " << iter().Kc(P, T) << endl;
     }
 
     Info<< nl << "end" << endl;
@@ -131,4 +131,3 @@ int main(int argc, char *argv[])
 
 
 // ************************************************************************* //
-
diff --git a/applications/utilities/thermophysical/equilibriumFlameT/equilibriumFlameT.C b/applications/utilities/thermophysical/equilibriumFlameT/equilibriumFlameT.C
index cfb5d71d74..8cf2371898 100644
--- a/applications/utilities/thermophysical/equilibriumFlameT/equilibriumFlameT.C
+++ b/applications/utilities/thermophysical/equilibriumFlameT/equilibriumFlameT.C
@@ -242,7 +242,7 @@ int main(int argc, char *argv[])
 
 
             scalar equilibriumFlameTemperatureNew =
-                products.THa(reactants.Ha(T0), adiabaticFlameTemperature);
+                products.THa(reactants.Ha(P, T0), P, adiabaticFlameTemperature);
 
             if (j==0)
             {
diff --git a/applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C b/applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C
index 02d4787f28..728897c88a 100644
--- a/applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C
+++ b/applications/utilities/thermophysical/mixtureAdiabaticFlameT/mixtureAdiabaticFlameT.C
@@ -69,6 +69,7 @@ int main(int argc, char *argv[])
     dictionary control(controlFile);
 
 
+    scalar P(readScalar(control.lookup("P")));
     scalar T0(readScalar(control.lookup("T0")));
     mixture rMix(control.lookup("reactants"));
     mixture pMix(control.lookup("products"));
@@ -116,7 +117,7 @@ int main(int argc, char *argv[])
     }
 
     Info<< "Adiabatic flame temperature of mixture " << rMix.name() << " = "
-         << products.THa(reactants.Ha(T0), 1000.0) << " K" << endl;
+         << products.THa(reactants.Ha(P, T0), P, 1000.0) << " K" << endl;
 
     return 0;
 }