diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
index ffef9994d4ab10fd5ed0be678b6d25514d623cea..3ac77ae34c758e7c999fe6875d3007b4445086fa 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/EulerDdtScheme/EulerDdtScheme.C
@@ -122,12 +122,10 @@ EulerDdtScheme<Type>::fvcDdt
             new GeometricField<Type, fvPatchField, volMesh>
             (
                 ddtIOobject,
-                mesh(),
-                rDeltaT.dimensions()*vf.dimensions(),
-                rDeltaT.value()*
+                rDeltaT*
                 (
-                    vf.primitiveField()
-                  - vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
+                    vf()
+                  - vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
                 ),
                 rDeltaT.value()*
                 (
@@ -174,12 +172,10 @@ EulerDdtScheme<Type>::fvcDdt
             new GeometricField<Type, fvPatchField, volMesh>
             (
                 ddtIOobject,
-                mesh(),
-                rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
-                rDeltaT.value()*rho.value()*
+                rDeltaT*rho*
                 (
-                    vf.primitiveField()
-                  - vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
+                    vf()
+                  - vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
                 ),
                 rDeltaT.value()*rho.value()*
                 (
@@ -226,13 +222,11 @@ EulerDdtScheme<Type>::fvcDdt
             new GeometricField<Type, fvPatchField, volMesh>
             (
                 ddtIOobject,
-                mesh(),
-                rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
-                rDeltaT.value()*
+                rDeltaT*
                 (
-                    rho.primitiveField()*vf.primitiveField()
-                  - rho.oldTime().primitiveField()
-                   *vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
+                    rho()*vf()
+                  - rho.oldTime()()
+                   *vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
                 ),
                 rDeltaT.value()*
                 (
@@ -282,18 +276,15 @@ EulerDdtScheme<Type>::fvcDdt
             new GeometricField<Type, fvPatchField, volMesh>
             (
                 ddtIOobject,
-                mesh(),
-                rDeltaT.dimensions()
-               *alpha.dimensions()*rho.dimensions()*vf.dimensions(),
                 rDeltaT.value()*
                 (
-                    alpha.primitiveField()
-                   *rho.primitiveField()
-                   *vf.primitiveField()
+                    alpha()
+                   *rho()
+                   *vf()
 
-                  - alpha.oldTime().primitiveField()
-                   *rho.oldTime().primitiveField()
-                   *vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
+                  - alpha.oldTime()()
+                   *rho.oldTime()()
+                   *vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
                 ),
                 rDeltaT.value()*
                 (