From abfe30454a973322dc08cb70c38ce3d2c6300ede Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
Date: Tue, 23 Jan 2024 15:26:04 +0000
Subject: [PATCH] ENH: finiteArea: consolidate various IOobject calls

---
 .../faMesh/faMeshTools/faMeshTools.C          |   2 +-
 .../EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C   |  48 ++++----
 .../EulerFaDdtScheme/EulerFaDdtScheme.C       |  96 ++++++++--------
 .../backwardFaDdtScheme/backwardFaDdtScheme.C |  96 ++++++++--------
 .../boundedBackwardFaDdtScheme.C              |  96 ++++++++--------
 .../steadyStateFaDdtScheme.C                  | 104 +++++++-----------
 6 files changed, 209 insertions(+), 233 deletions(-)

diff --git a/src/finiteArea/faMesh/faMeshTools/faMeshTools.C b/src/finiteArea/faMesh/faMeshTools/faMeshTools.C
index d9a346a52d1..0154e992980 100644
--- a/src/finiteArea/faMesh/faMeshTools/faMeshTools.C
+++ b/src/finiteArea/faMesh/faMeshTools/faMeshTools.C
@@ -144,7 +144,7 @@ Foam::faMeshTools::newMesh
     IOobject cmptIO(meshIO, "faceLabels", meshSubDir);
     cmptIO.readOpt(IOobject::MUST_READ);
     cmptIO.writeOpt(IOobject::NO_WRITE);
-    cmptIO.registerObject(false);
+    cmptIO.registerObject(IOobject::NO_REGISTER);
 
 
     // Check who has a mesh
diff --git a/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C b/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
index 0788c39e60d..e80016d9360 100644
--- a/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
+++ b/src/finiteArea/finiteArea/d2dt2Schemes/EulerFaD2dt2Scheme/EulerFaD2dt2Scheme.C
@@ -71,13 +71,13 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     scalar coefft   = (deltaT + deltaT0)/(2*deltaT);
     scalar coefft00 = (deltaT + deltaT0)/(2*deltaT0);
 
-    IOobject d2dt2IOobject
+    const IOobject d2dt2IOobject
     (
-        "d2dt2("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "d2dt2("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -130,13 +130,13 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     dimensionedScalar rDeltaT2 =
         4.0/sqr(mesh().time().deltaT() + mesh().time().deltaT0());
 
-    IOobject d2dt2IOobject
+    const IOobject d2dt2IOobject
     (
-        "d2dt2("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "d2dt2("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -208,13 +208,13 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     dimensionedScalar rDeltaT2 =
         4.0/sqr(mesh().time().deltaT() + mesh().time().deltaT0());
 
-    IOobject d2dt2IOobject
+    const IOobject d2dt2IOobject
     (
-        "d2dt2("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "d2dt2("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -293,13 +293,13 @@ EulerFaD2dt2Scheme<Type>::facD2dt2
     dimensionedScalar rDeltaT2 =
         4.0/sqr(mesh().time().deltaT() + mesh().time().deltaT0());
 
-    IOobject d2dt2IOobject
+    const IOobject d2dt2IOobject
     (
-        "d2dt2("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "d2dt2("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
diff --git a/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
index 33452b4ccb5..3359014fe4b 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/EulerFaDdtScheme/EulerFaDdtScheme.C
@@ -49,13 +49,13 @@ EulerFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -96,13 +96,13 @@ EulerFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     tmp<GeometricField<Type, faPatchField, areaMesh>> tdtdt0
@@ -134,13 +134,13 @@ EulerFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -187,13 +187,13 @@ EulerFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -234,13 +234,13 @@ EulerFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -287,13 +287,13 @@ EulerFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -336,13 +336,13 @@ EulerFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
@@ -393,13 +393,13 @@ EulerFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     if (mesh().moving())
diff --git a/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
index 98d1889a46c..457e744df21 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/backwardFaDdtScheme/backwardFaDdtScheme.C
@@ -80,13 +80,13 @@ backwardFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -135,13 +135,13 @@ backwardFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -182,13 +182,13 @@ backwardFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -255,13 +255,13 @@ backwardFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -326,13 +326,13 @@ backwardFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -399,13 +399,13 @@ backwardFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -470,13 +470,13 @@ backwardFaDdtScheme<Type>::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -547,13 +547,13 @@ backwardFaDdtScheme<Type>::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
diff --git a/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
index 32aeac7c158..2e41f528971 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/boundedBackwardFaDdtScheme/boundedBackwardFaDdtScheme.C
@@ -63,13 +63,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
 
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -118,13 +118,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
 
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+dt.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+dt.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -163,13 +163,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -260,13 +260,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -355,13 +355,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -452,13 +452,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -547,13 +547,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
@@ -648,13 +648,13 @@ tmp<areaScalarField> boundedBackwardFaDdtScheme::facDdt0
 {
     dimensionedScalar rDeltaT = 1.0/mesh().time().deltaT();
 
-    IOobject ddtIOobject
+    const IOobject ddtIOobject
     (
-        "ddt0("+rho.name()+','+vf.name()+')',
-        mesh().time().timeName(),
-        mesh().thisDb(),
-        IOobject::NO_READ,
-        IOobject::NO_WRITE
+        mesh().thisDb().newIOobject
+        (
+            "ddt0("+rho.name()+','+vf.name()+')',
+            { IOobject::REGISTER }
+        )
     );
 
     scalar deltaT = deltaT_();
diff --git a/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C b/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
index d430680d880..c562a9b2401 100644
--- a/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
+++ b/src/finiteArea/finiteArea/ddtSchemes/steadyStateFaDdtScheme/steadyStateFaDdtScheme.C
@@ -48,16 +48,13 @@ steadyStateFaDdtScheme<Type>::facDdt
     const dimensioned<Type> dt
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt("+dt.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt("+dt.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(dt.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        dt.dimensions()/dimTime
     );
 }
 
@@ -69,16 +66,13 @@ steadyStateFaDdtScheme<Type>::facDdt0
     const dimensioned<Type> dt
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt("+dt.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt("+dt.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(dt.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        dt.dimensions()/dimTime
     );
 }
 
@@ -90,16 +84,13 @@ steadyStateFaDdtScheme<Type>::facDdt
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt("+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt("+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        vf.dimensions()/dimTime
     );
 }
 
@@ -111,16 +102,13 @@ steadyStateFaDdtScheme<Type>::facDdt0
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt0("+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt0("+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        vf.dimensions()/dimTime
     );
 }
 
@@ -133,16 +121,13 @@ steadyStateFaDdtScheme<Type>::facDdt
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt("+rho.name()+','+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt("+rho.name()+','+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        rho.dimensions()*vf.dimensions()/dimTime
     );
 }
 
@@ -154,16 +139,13 @@ steadyStateFaDdtScheme<Type>::facDdt0
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt0("+rho.name()+','+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt0("+rho.name()+','+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        rho.dimensions()*vf.dimensions()/dimTime
     );
 }
 
@@ -176,16 +158,13 @@ steadyStateFaDdtScheme<Type>::facDdt
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt("+rho.name()+','+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt("+rho.name()+','+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        rho.dimensions()*vf.dimensions()/dimTime
     );
 }
 
@@ -198,16 +177,13 @@ steadyStateFaDdtScheme<Type>::facDdt0
     const GeometricField<Type, faPatchField, areaMesh>& vf
 )
 {
-    return tmp<GeometricField<Type, faPatchField, areaMesh>>::New
+    return GeometricField<Type, faPatchField, areaMesh>::New
     (
-        IOobject
-        (
-            "ddt0("+rho.name()+','+vf.name()+')',
-            mesh().time().timeName(),
-            mesh().thisDb()
-        ),
+        "ddt0("+rho.name()+','+vf.name()+')',
+        IOobject::REGISTER,
         mesh(),
-        dimensioned<Type>(rho.dimensions()*vf.dimensions()/dimTime, Zero)
+        Foam::zero{},
+        rho.dimensions()*vf.dimensions()/dimTime
     );
 }
 
-- 
GitLab