diff --git a/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C b/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
index ff371e230393c12ff75670e18d2902f3f6475a1c..dca99a5a61afae4429aea198cb964b39f2f32686 100644
--- a/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
+++ b/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
@@ -111,7 +111,7 @@ Foam::Ostream& Foam::operator<<
 {
     os  << static_cast<const EquationOfState&>(pt) << tab
         << pt.Hf_/pt.W() << tab
-        << pt.Sf_ << tab
+        << pt.Sf_/pt.W() << tab
         << "cpPolynomial" << tab << pt.cpPolynomial_/pt.W();
 
     os.check