diff --git a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.C b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.C
index 4704bb39..e2b99d08 100644
--- a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.C
+++ b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.C
@@ -53,6 +53,16 @@ Foam::PengRobinsonGas<Specie>::PengRobinsonGas
 {}
 
 
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+
+template<class Specie>
+void Foam::PengRobinsonGas<Specie>::write(Ostream& os) const
+{
+    Specie::write(os);
+}
+
+
 // * * * * * * * * * * * * * * * Ostream Operator  * * * * * * * * * * * * * //
 
 template<class Specie>
diff --git a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.H b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.H
index fd141117..7d1df112 100644
--- a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.H
+++ b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGas.H
@@ -220,6 +220,10 @@ public:
 };
 
 
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 #include "PengRobinsonGasI.H"
diff --git a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGasI.H b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGasI.H
index 861f2e09..aae0b92f 100644
--- a/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGasI.H
+++ b/src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGasI.H
@@ -218,6 +218,7 @@ inline void Foam::PengRobinsonGas<Specie>::operator+=
     omega_ = molr1*omega_ + molr2*pg.omega_;
 }
 
+
 template<class Specie>
 inline void Foam::PengRobinsonGas<Specie>::operator-=
 (
@@ -236,6 +237,7 @@ inline void Foam::PengRobinsonGas<Specie>::operator-=
     omega_ = molr1*omega_ - molr2*pg.omega_;
 }
 
+
 template<class Specie>
 inline void Foam::PengRobinsonGas<Specie>::operator*=(const scalar s)
 {