STYLE: use zeroGradientType() instead of hard-coded or typeName

src/atmosphericModels/functionObjects/ObukhovLength/ObukhovLength.C

@@ -76,7 +76,7 @@ bool Foam::functionObjects::ObukhovLength::calcOL()
             ),
             mesh_,
             dimLength,
-            zeroGradientFvPatchField<scalar>::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         result1->store();
@@ -94,7 +94,7 @@ bool Foam::functionObjects::ObukhovLength::calcOL()
             ),
             mesh_,
             dimVelocity,
-            zeroGradientFvPatchField<scalar>::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         result2->store();

src/combustionModels/diffusionMulticomponent/diffusionMulticomponent.C

@@ -71,7 +71,7 @@ diffusionMulticomponent<ReactionThermo, ThermoType>::init()
                 ),
                 this->mesh_,
                 dimensionedScalar(dimMass/dimTime/dimVolume, Zero),
-                zeroGradientFvPatchScalarField::typeName
+                fvPatchFieldBase::zeroGradientType()
             )
         );
 
@@ -209,7 +209,7 @@ diffusionMulticomponent<ReactionThermo, ThermoType>::correct()
                     ),
                     this->mesh_,
                     dimensionedScalar(dimMass/dimTime/dimVolume, Zero),
-                    zeroGradientFvPatchScalarField::typeName
+                    fvPatchFieldBase::zeroGradientType()
                 )
             );
 
@@ -410,7 +410,7 @@ diffusionMulticomponent<ReactionThermo, ThermoType>::Qdot() const
             ),
             this->mesh(),
             dimensionedScalar(dimEnergy/dimTime/dimVolume, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         )
     );
 

src/combustionModels/eddyDissipationDiffusionModel/eddyDissipationDiffusionModel.C

@@ -91,7 +91,7 @@ eddyDissipationDiffusionModel<ReactionThermo, ThermoType>::rtDiff() const
             ),
             this->mesh(),
             dimensionedScalar(dimLength, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         )
     );
 

src/dynamicFaMesh/interfaceTrackingFvMesh/interfaceTrackingFvMesh.C

@@ -144,7 +144,7 @@ void Foam::interfaceTrackingFvMesh::makeUs() const
     wordList patchFieldTypes
     (
         aMesh().boundary().size(),
-        zeroGradientFaPatchVectorField::typeName
+        faPatchFieldBase::zeroGradientType()
     );
 
     forAll(aMesh().boundary(), patchI)
@@ -883,7 +883,7 @@ void Foam::interfaceTrackingFvMesh::updateSurfactantConcentration()
                 ),
                 aMesh(),
                 dimensionedScalar(dimMoles/dimVolume, Zero),
-                zeroGradientFaPatchScalarField::typeName
+                faPatchFieldBase::zeroGradientType()
             );
 
             Cb.ref().field() =

src/finiteArea/finiteArea/fac/facEdgeIntegrate.C

@@ -62,7 +62,7 @@ edgeIntegrate
             ),
             mesh,
             dimensioned<Type>(ssf.dimensions()/dimArea, Zero),
-            zeroGradientFaPatchField<Type>::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
     GeometricField<Type, faPatchField, areaMesh>& vf = tvf.ref();
@@ -134,7 +134,7 @@ edgeSum
             ),
             mesh,
             dimensioned<Type>(ssf.dimensions(), Zero),
-            zeroGradientFaPatchField<Type>::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
     GeometricField<Type, faPatchField, areaMesh>& vf = tvf.ref();

src/finiteArea/finiteArea/gradSchemes/leastSquaresFaGrad/leastSquaresFaGrad.C

@@ -77,7 +77,7 @@ leastSquaresFaGrad<Type>::grad
             ),
             mesh,
             dimensioned<GradType>(vsf.dimensions()/dimLength, Zero),
-            zeroGradientFaPatchField<GradType>::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
     GeometricField<GradType, faPatchField, areaMesh>& lsGrad = tlsGrad.ref();

src/finiteVolume/cfdTools/general/CorrectPhi/CorrectPhi.C

@@ -60,7 +60,7 @@ void Foam::CorrectPhi
     wordList pcorrTypes
     (
         p.boundaryField().size(),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     // Set BCs of pcorr to fixed-value for patches at which p is fixed
@@ -140,7 +140,7 @@ void Foam::CorrectPhi
     wordList pcorrTypes
     (
         p.boundaryField().size(),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     // Set BCs of pcorr to fixed-value for patches at which p is fixed

src/finiteVolume/finiteVolume/fvc/fvcSmooth/fvcSmooth.C

@@ -345,7 +345,7 @@ void Foam::fvc::spreadSource
         ),
         mesh,
         dimensionedScalar(mDotOut.dimensions(), Zero),
-        zeroGradientFvPatchField<scalar>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     // Smearing of source term field

src/finiteVolume/fvMesh/fvMeshSubset/fvMeshSubsetProxyTemplates.C

@@ -51,7 +51,7 @@ Foam::fvMeshSubsetProxy::zeroGradientField
         io,
         df.mesh(),
         dimensioned<Type>(df.dimensions(), Zero),
-        zeroGradientFvPatchField<Type>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     tfield.ref().primitiveFieldRef() = df;
     tfield.ref().oriented() = df.oriented();

src/finiteVolume/fvMesh/wallDist/patchDistMethods/patchDistMethod/patchDistMethodTemplates.C

@@ -40,7 +40,7 @@ Foam::wordList Foam::patchDistMethod::patchTypes
     wordList yTypes
     (
         mesh.boundary().size(),
-        zeroGradientFvPatchField<Type>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     for (const label patchi : patchIDs)

src/functionObjects/field/CourantNo/CourantNo.C

@@ -100,7 +100,7 @@ bool Foam::functionObjects::CourantNo::calc()
                 ),
                 mesh_,
                 dimensionedScalar(dimless, Zero),
-                zeroGradientFvPatchScalarField::typeName
+                fvPatchFieldBase::zeroGradientType()
             );
             tCo.ref().ref() = Coi();
             tCo.ref().correctBoundaryConditions();

src/functionObjects/field/age/age.C

@@ -63,7 +63,7 @@ Foam::wordList Foam::functionObjects::age::patchTypes() const
     {
         if (isA<wallFvPatch>(mesh_.boundary()[patchi]))
         {
-            result[patchi] = zeroGradientFvPatchField<scalar>::typeName;
+            result[patchi] = fvPatchFieldBase::zeroGradientType();
         }
     }
 

src/functionObjects/field/blendingFactor/blendingFactor.C

@@ -96,7 +96,7 @@ Foam::functionObjects::blendingFactor::blendingFactor
         ),
         mesh_,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     store(resultName_, indicatorPtr);

src/functionObjects/field/resolutionIndex/resolutionIndexModels/resolutionIndexModel/resolutionIndexModel.C

@@ -58,7 +58,7 @@ Foam::tmp<Foam::volScalarField> Foam::resolutionIndexModel::V() const
         ),
         mesh_,
         dimVolume,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     tV.ref().primitiveFieldRef() = mesh_.V();
@@ -104,7 +104,7 @@ bool Foam::resolutionIndexModel::read(const dictionary& dict)
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         mesh_.objectRegistry::store(indexPtr);

src/functionObjects/field/stabilityBlendingFactor/stabilityBlendingFactor.C

@@ -71,7 +71,7 @@ Foam::functionObjects::stabilityBlendingFactor::indicator()
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         mesh_.objectRegistry::store(fldPtr);
@@ -332,7 +332,7 @@ bool Foam::functionObjects::stabilityBlendingFactor::init(bool first)
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
         auto& magGradCC = tmagGradCC.ref();
 
@@ -351,7 +351,7 @@ bool Foam::functionObjects::stabilityBlendingFactor::init(bool first)
                 ),
                 mesh_,
                 dimensionedScalar(dimLength, Zero),
-                zeroGradientFvPatchScalarField::typeName
+                fvPatchFieldBase::zeroGradientType()
             );
             cci = mesh_.C().component(i);
             cci.correctBoundaryConditions();
@@ -405,7 +405,7 @@ bool Foam::functionObjects::stabilityBlendingFactor::init(bool first)
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         auto& Co = CoPtr.ref();

src/functionObjects/field/zeroGradient/zeroGradientTemplates.C

@@ -96,7 +96,7 @@ int Foam::functionObjects::zeroGradient::apply
             ),
             mesh_,
             dimensioned<Type>(input.dimensions(), Zero),
-            zeroGradientFvPatchField<Type>::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         store(outputName, tzeroGrad);

src/functionObjects/utilities/solverInfo/solverInfo.C

@@ -217,7 +217,7 @@ bool Foam::functionObjects::solverInfo::write()
                 ),
                 mesh_,
                 dimensionedScalar(dimless, Zero),
-                zeroGradientFvPatchField<scalar>::typeName
+                fvPatchFieldBase::zeroGradientType()
             );
 
             residual.primitiveFieldRef() = *residualPtr;

src/fvMotionSolver/motionDiffusivity/inverseVolume/inverseVolumeDiffusivity.C

@@ -84,7 +84,7 @@ void Foam::inverseVolumeDiffusivity::correct()
         ),
         mesh(),
         dimless,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     V.primitiveFieldRef() = mesh().V();

src/fvMotionSolver/motionDiffusivity/motionDirectional/motionDirectionalDiffusivity.C

@@ -87,7 +87,8 @@ void Foam::motionDirectionalDiffusivity::correct()
             diffusivityVector_.y()*vector::one
           + (diffusivityVector_.x() - diffusivityVector_.y())*cellMotionU
            /(mag(cellMotionU) + dimensionedScalar("small", dimVelocity, SMALL)),
-            zeroGradientFvPatchVectorField::typeName
+
+            fvPatchFieldBase::zeroGradientType()
         );
         D.correctBoundaryConditions();
 

src/fvMotionSolver/motionInterpolation/patchCorrected/patchCorrectedInterpolationTemplates.C

@@ -65,7 +65,7 @@ void Foam::patchCorrectedInterpolation::interpolateType
             wordList
             (
                 pointUncorrectedDisplacement.boundaryField().size(),
-                zeroGradientPointPatchField<Type>::typeName
+                pointPatchFieldBase::zeroGradientType()
             )
         );
 
@@ -107,7 +107,7 @@ void Foam::patchCorrectedInterpolation::interpolateDataFromPatchGroups
         ),
         data.mesh(),
         dimensionedScalar(dimless, Zero),
-        zeroGradientPointPatchField<scalar>::typeName
+        pointPatchFieldBase::zeroGradientType()
     );
     data = dimensioned<Type>(data.dimensions(), Zero);
 
@@ -124,7 +124,7 @@ void Foam::patchCorrectedInterpolation::interpolateDataFromPatchGroups
             ),
             data.mesh(),
             dimensionedScalar(data.dimensions(), Zero),
-            zeroGradientPointPatchField<scalar>::typeName
+            pointPatchFieldBase::zeroGradientType()
         );
         GeometricField<Type, pointPatchField, pointMesh> patchData(data);
 

src/fvOptions/sources/derived/acousticDampingSource/acousticDampingSource.C

@@ -98,8 +98,9 @@ Foam::fv::acousticDampingSource::acousticDampingSource
                 IOobject::NO_WRITE
             ),
             mesh_,
-            dimensionedScalar("blend0", dimless, 1.0),
+            scalar(1),
-            zeroGradientFvPatchField<vector>::typeName
+            dimless,
+            fvPatchFieldBase::zeroGradientType()
         )
     ),
     frequency_("frequency", dimless/dimTime, 0),

src/fvOptions/sources/derived/jouleHeatingSource/jouleHeatingSource.C

@@ -68,7 +68,7 @@ Foam::fv::jouleHeatingSource::transformSigma
         ),
         mesh_,
         dimensionedSymmTensor(sigmaLocal.dimensions(), Zero),
-        zeroGradientFvPatchField<symmTensor>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     auto& sigma = tsigma.ref();
 

src/fvOptions/sources/derived/solidificationMeltingSource/solidificationMeltingSource.C

@@ -192,7 +192,7 @@ Foam::fv::solidificationMeltingSource::solidificationMeltingSource
         ),
         mesh,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     curTimeIndex_(-1),
     deltaT_(cells_.size(), 0)

src/fvOptions/sources/interRegion/interRegionHeatTransfer/interRegionHeatTransferModel/interRegionHeatTransferModel.C

@@ -139,7 +139,7 @@ Foam::fv::interRegionHeatTransferModel::interRegionHeatTransferModel
         ),
         mesh,
         dimensionedScalar(dimEnergy/dimTime/dimTemperature/dimVolume, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     TName_(coeffs_.getOrDefault<word>("T", "T")),
     TNbrName_(coeffs_.getOrDefault<word>("TNbr", "T"))

src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C

@@ -717,7 +717,7 @@ Foam::DSMCCloud<ParcelType>::DSMCCloud
         ),
         mesh_,
         dimensionedScalar(dimensionSet(0, 3, -1, 0, 0), Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     collisionSelectionRemainder_
     (

src/lagrangian/intermediate/submodels/MPPIC/AveragingMethods/Basic/Basic.C

@@ -80,7 +80,7 @@ void Foam::AveragingMethods::Basic<Type>::updateGrad()
         ),
         this->mesh_,
         dimensioned<Type>(dimless, Zero),
-        zeroGradientFvPatchField<Type>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     tempData.primitiveFieldRef() = data_;
     tempData.correctBoundaryConditions();

src/lagrangian/intermediate/submodels/MPPIC/PackingModels/Implicit/Implicit.C

@@ -56,7 +56,7 @@ Foam::PackingModels::Implicit<CloudType>::Implicit
         ),
         this->owner().mesh(),
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     phiCorrect_(nullptr),
     uCorrect_(nullptr),
@@ -150,7 +150,7 @@ void Foam::PackingModels::Implicit<CloudType>::cacheFields(const bool store)
             ),
             mesh,
             dimensionedScalar(dimDensity, Zero),
-            zeroGradientFvPatchField<scalar>::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
         rho.primitiveFieldRef() = max(rhoAverage.primitiveField(), rhoMin_);
         rho.correctBoundaryConditions();
@@ -171,7 +171,7 @@ void Foam::PackingModels::Implicit<CloudType>::cacheFields(const bool store)
             ),
             mesh,
             dimensionedScalar(dimPressure, Zero),
-            zeroGradientFvPatchField<scalar>::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         tauPrime.primitiveFieldRef() =
@@ -256,7 +256,7 @@ void Foam::PackingModels::Implicit<CloudType>::cacheFields(const bool store)
                 ),
                 mesh,
                 dimensionedVector(dimVelocity, Zero),
-                fixedValueFvPatchField<vector>::typeName
+                fvPatchFieldBase::zeroGradientType()
             );
             U.primitiveFieldRef() = uAverage.primitiveField();
             U.correctBoundaryConditions();

src/lagrangian/molecularDynamics/molecule/mdTools/createMDFields.H

@@ -57,7 +57,7 @@ forAll(allSpeciesRhoN, rN)
             ),
             mesh,
             dimless/dimVolume,
-            "zeroGradient"
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     allSpeciesRhoN[rN].primitiveFieldRef() = 0;
@@ -78,7 +78,7 @@ volScalarField totalRhoN
     ),
     mesh,
     dimless/dimVolume,
-    "zeroGradient"
+    fvPatchFieldBase::zeroGradientType()
 );
 totalRhoN.primitiveFieldRef() = 0;
 totalRhoN.correctBoundaryConditions();
@@ -112,7 +112,7 @@ forAll(allSpeciesRhoM, rM)
             ),
             mesh,
             dimMass/dimVolume,
-            "zeroGradient"
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     allSpeciesRhoM[rM].primitiveFieldRef() = 0;
@@ -133,7 +133,7 @@ volScalarField totalRhoM
     ),
     mesh,
     dimMass/dimVolume,
-    "zeroGradient"
+    fvPatchFieldBase::zeroGradientType()
 );
 totalRhoM.primitiveFieldRef() = 0;
 totalRhoM.correctBoundaryConditions();
@@ -167,7 +167,7 @@ forAll(allSpeciesVelocity, v)
             ),
             mesh,
             dimVelocity,
-            "zeroGradient"
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     allSpeciesVelocity[v].primitiveFieldRef() = Zero;
@@ -188,7 +188,7 @@ Info << "    Creating total velocity field" << endl;
 //     ),
 //     mesh,
 //     dimVelocity,
-//     "zeroGradient"
+//     fvPatchFieldBase::zeroGradientType()
 // );
 // totalVelocity.primitiveFieldRef() = Zero;
 // totalVelocity.correctBoundaryConditions();
@@ -237,7 +237,7 @@ forAll(allSpeciesTemperature, t)
             ),
             mesh,
             dimTemperature,
-            "zeroGradient"
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     allSpeciesTemperature[t].primitiveFieldRef() = 0;
@@ -258,7 +258,7 @@ volScalarField totalTemperature
     ),
     mesh,
     dimTemperature,
-    "zeroGradient"
+    fvPatchFieldBase::zeroGradientType()
 );
 totalTemperature.primitiveFieldRef() = 0;
 totalTemperature.correctBoundaryConditions();
@@ -293,7 +293,7 @@ forAll(allSpeciesMeanKE, mKE)
             ),
             mesh,
             dimensionSet(1, 2, -2, 0, 0, 0, 0),
-            "zeroGradient"
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     allSpeciesMeanKE[mKE].primitiveFieldRef() = 0;
@@ -314,7 +314,7 @@ volScalarField totalMeanKE
     ),
     mesh,
     dimensionSet(1, 2, -2, 0, 0, 0, 0),
-    "zeroGradient"
+    fvPatchFieldBase::zeroGradientType()
 );
 totalMeanKE.primitiveFieldRef() = 0;
 totalMeanKE.correctBoundaryConditions();

src/mesh/snappyHexMesh/meshRefinement/meshRefinement.C

@@ -2406,27 +2406,27 @@ Foam::label Foam::meshRefinement::appendPatch
     addPatchFields<volScalarField>
     (
         mesh,
-        calculatedFvPatchField<scalar>::typeName
+        fvPatchFieldBase::calculatedType()
     );
     addPatchFields<volVectorField>
     (
         mesh,
-        calculatedFvPatchField<vector>::typeName
+        fvPatchFieldBase::calculatedType()
     );
     addPatchFields<volSphericalTensorField>
     (
         mesh,
-        calculatedFvPatchField<sphericalTensor>::typeName
+        fvPatchFieldBase::calculatedType()
     );
     addPatchFields<volSymmTensorField>
     (
         mesh,
-        calculatedFvPatchField<symmTensor>::typeName
+        fvPatchFieldBase::calculatedType()
     );
     addPatchFields<volTensorField>
     (
         mesh,
-        calculatedFvPatchField<tensor>::typeName
+        fvPatchFieldBase::calculatedType()
     );
 
     // Surface fields
@@ -2434,27 +2434,27 @@ Foam::label Foam::meshRefinement::appendPatch
     addPatchFields<surfaceScalarField>
     (
         mesh,
-        calculatedFvPatchField<scalar>::typeName
+        fvsPatchFieldBase::calculatedType()
     );
     addPatchFields<surfaceVectorField>
     (
         mesh,
-        calculatedFvPatchField<vector>::typeName
+        fvsPatchFieldBase::calculatedType()
     );
     addPatchFields<surfaceSphericalTensorField>
     (
         mesh,
-        calculatedFvPatchField<sphericalTensor>::typeName
+        fvsPatchFieldBase::calculatedType()
     );
     addPatchFields<surfaceSymmTensorField>
     (
         mesh,
-        calculatedFvPatchField<symmTensor>::typeName
+        fvsPatchFieldBase::calculatedType()
     );
     addPatchFields<surfaceTensorField>
     (
         mesh,
-        calculatedFvPatchField<tensor>::typeName
+        fvsPatchFieldBase::calculatedType()
     );
     return patchi;
 }
@@ -2574,7 +2574,7 @@ Foam::label Foam::meshRefinement::addMeshedPatch
 //            mesh_,
 //            ppPtr(),
 //            dictionary(),       // optional field values
-//            calculatedFvPatchField<scalar>::typeName,
+//            fvPatchFieldBase::calculatedType(),
 //            true
 //        );
 
@@ -3595,7 +3595,7 @@ void Foam::meshRefinement::dumpRefinementLevel() const
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         const labelList& cellLevel = meshCutter_.cellLevel();

src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C

@@ -3576,7 +3576,7 @@ void Foam::meshRefinement::zonify
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         forAll(cellToZone, cellI)

src/optimisation/adjointOptimisation/adjoint/ATCModel/ATCModel/ATCModel.C

@@ -101,8 +101,9 @@ ATCModel::ATCModel
             IOobject::NO_WRITE
         ),
         mesh_,
-        dimensionedScalar("limiter", dimless, 1.0),
+        scalar(1),
-        zeroGradientFvPatchField<scalar>::typeName
+        dimless,
+        fvPatchFieldBase::zeroGradientType()
     ),
     ATC_
     (
@@ -231,17 +232,18 @@ tmp<volScalarField> ATCModel::createLimiter
     (
         new volScalarField
         (
-           IOobject
+            IOobject
-           (
+            (
-               "limiter",
+                "limiter",
-               mesh.time().timeName(),
+                mesh.time().timeName(),
-               mesh,
+                mesh,
-               IOobject::NO_READ,
+                IOobject::NO_READ,
-               IOobject::NO_WRITE
+                IOobject::NO_WRITE
-           ),
+            ),
-           mesh,
+            mesh,
-           dimensionedScalar("limiter", dimless, 1.0),
+            scalar(1),
-           zeroGradientFvPatchField<scalar>::typeName
+            dimless,
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     volScalarField& limiter = tlimiter.ref();

src/optimisation/adjointOptimisation/adjoint/ATCModel/ATCModel/ATCModelTemplates.C

@@ -50,8 +50,9 @@ void Foam::ATCModel::smoothFieldBasedOnCells
             IOobject::NO_WRITE
         ),
         mesh_,
-        dimensionedScalar("limiter", dimless, 1.0),
+        scalar(1),
-        zeroGradientFvPatchField<scalar>::typeName
+        dimless,
+        fvPatchFieldBase::zeroGradientType()
     );
 
     computeLimiter(limiter, cells, nSmooth_);

src/optimisation/adjointOptimisation/adjoint/dynamicMesh/motionSolver/elasticityMotionSolver/elasticityMotionSolver.C

@@ -153,7 +153,7 @@ Foam::elasticityMotionSolver::elasticityMotionSolver
         ),
         fvMesh_,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     exponent_(this->coeffDict().get<scalar>("exponent")),
     nSteps_(this->coeffDict().get<label>("steps")),

src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointEikonalSolver/adjointEikonalSolverIncompressible.C

@@ -55,7 +55,7 @@ wordList adjointEikonalSolver::patchTypes() const
 
     for (const label patchi : wallPatchIDs_)
     {
-        daTypes[patchi] = zeroGradientFvPatchScalarField::typeName;
+        daTypes[patchi] = fvPatchFieldBase::zeroGradientType();
     }
 
     return daTypes;

src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurface/sensitivitySurfaceIncompressible.C

@@ -320,7 +320,7 @@ void sensitivitySurface::smoothSensitivities()
         ),
         aMesh,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFaPatchField<scalar>::typeName
+        faPatchFieldBase::zeroGradientType()
     );
 
     // Copy sensitivities to area field

src/optimisation/adjointOptimisation/adjoint/turbulenceModels/incompressibleAdjoint/adjointRAS/adjointkOmegaSST/adjointkOmegaSST.C

@@ -2293,7 +2293,7 @@ tmp<volTensorField> adjointkOmegaSST::FISensitivityTerm()
             ),
             mesh_,
             dimensionedTensor(sqr(dimLength)/pow3(dimTime), Zero),
-            zeroGradientFvPatchTensorField::typeName
+            fvPatchFieldBase::zeroGradientType()
         )
     );
     volTensorField& FISens = tFISens.ref();

src/overset/cellCellStencil/cellCellStencil/cellCellStencilTemplates.C

@@ -171,7 +171,7 @@ Foam::cellCellStencil::createField
         ),
         mesh,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     volScalarField& fld = tfld.ref();
 

src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C

@@ -534,7 +534,7 @@ Foam::cellCellStencils::cellVolumeWeight::cellVolumeWeight
         ),
         mesh_,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     allowInterpolatedDonors_
     (

src/overset/cellCellStencil/inverseDistance/inverseDistanceCellCellStencil.C

@@ -1736,7 +1736,7 @@ Foam::cellCellStencils::inverseDistance::inverseDistance
         ),
         mesh_,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     )
 {
     // Protect local fields from interpolation

src/overset/include/createCellMask.H

@@ -46,8 +46,9 @@ volScalarField cellMask
         IOobject::NO_WRITE
     ),
     mesh,
-    dimensionedScalar("one", dimless, 1.0),
+    scalar(1),
-    zeroGradientFvPatchScalarField::typeName
+    dimless,
+    fvPatchFieldBase::zeroGradientType()
 );
 
 surfaceScalarField faceMask

src/overset/include/createInterpolatedCells.H

@@ -46,8 +46,9 @@ volScalarField interpolatedCells
         IOobject::NO_WRITE
     ),
     mesh,
-    dimensionedScalar("one", dimless, 1.0),
+    scalar(1),
-    zeroGradientFvPatchScalarField::typeName
+    dimless,
+    fvPatchFieldBase::zeroGradientType()
 );
 
 #include "setInterpolatedCells.H"

src/overset/oversetFvMesh/oversetFvMeshBase.C

@@ -656,7 +656,7 @@ bool Foam::oversetFvMeshBase::writeObject
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         forAll(volTypes.internalField(), cellI)
@@ -680,7 +680,7 @@ bool Foam::oversetFvMeshBase::writeObject
             ),
             mesh_,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
 
         const cellCellStencilObject& overlap = Stencil::New(mesh_);
@@ -719,8 +719,9 @@ bool Foam::oversetFvMeshBase::writeObject
                 IOobject::NO_REGISTER
             ),
             mesh_,
-            dimensionedScalar("minOne", dimless, scalar(-1)),
+            scalar(-1),
-            zeroGradientFvPatchScalarField::typeName
+            dimless,
+            fvPatchFieldBase::zeroGradientType()
         );
 
         forAll(cellStencil, cellI)

src/phaseSystemModels/multiphaseInter/phasesSystem/InterfaceCompositionModel/InterfaceCompositionModel.C

@@ -90,7 +90,7 @@ Foam::InterfaceCompositionModel<Thermo, OtherThermo>::getSpecieMassFraction
         ),
         mesh,
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     auto& Ys = tY.ref();
@@ -123,8 +123,9 @@ Foam::InterfaceCompositionModel<Thermo, OtherThermo>::getSpecieMassFraction
             IOobject::NO_WRITE
         ),
         mesh,
-        dimensionedScalar("one", dimless, scalar(1)),
+        scalar(1),
-        zeroGradientFvPatchScalarField::typeName
+        dimless,
+        fvPatchFieldBase::zeroGradientType()
     );
 }
 
@@ -156,7 +157,7 @@ Foam::InterfaceCompositionModel<Thermo, OtherThermo>::MwMixture
             dimMass/dimMoles,
             1e-3*mixture.cellMixture(0).W()
         ),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 }
 
@@ -328,7 +329,7 @@ Foam::InterfaceCompositionModel<Thermo, OtherThermo>::L
         ),
         p.mesh(),
         dimensionedScalar(dimEnergy/dimMass, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     auto& L = tmpL.ref();

src/phaseSystemModels/reactingEuler/multiphaseSystem/interfacialModels/dragModels/segregated/segregated.C

@@ -104,7 +104,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
         ),
         mesh,
         dimensionedScalar(dimLength),
-        zeroGradientFvPatchField<scalar>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     L.primitiveFieldRef() = cbrt(mesh.V());
     L.correctBoundaryConditions();

src/phaseSystemModels/twoPhaseEuler/twoPhaseSystem/interfacialModels/dragModels/segregated/segregated.C

@@ -104,7 +104,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::segregated::K() const
         ),
         mesh,
         dimensionedScalar(dimLength, Zero),
-        zeroGradientFvPatchField<scalar>::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
     L.primitiveFieldRef() = cbrt(mesh.V());
     L.correctBoundaryConditions();

src/regionFaModels/KirchhoffShell/KirchhoffShell.C

@@ -297,7 +297,7 @@ const tmp<areaScalarField> KirchhoffShell::rho() const
             ),
             regionMesh(),
             dimensionedScalar("rho", dimDensity, solid().rho()),
-            zeroGradientFaPatchScalarField::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
 }

src/regionFaModels/thermalShell/thermalShell.C

@@ -199,7 +199,7 @@ const tmp<areaScalarField> thermalShell::Cp() const
             ),
             regionMesh(),
             dimensionedScalar(dimEnergy/dimTemperature/dimMass, thermo_.Cp()),
-            zeroGradientFaPatchScalarField::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
 }
@@ -222,7 +222,7 @@ const tmp<areaScalarField> thermalShell::rho() const
             ),
             regionMesh(),
             dimensionedScalar(dimDensity, thermo_.rho()),
-            zeroGradientFaPatchScalarField::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
 }
@@ -249,7 +249,7 @@ const tmp<areaScalarField> thermalShell::kappa() const
                 dimPower/dimLength/dimTemperature,
                 thermo_.kappa()
             ),
-            zeroGradientFaPatchScalarField::typeName
+            faPatchFieldBase::zeroGradientType()
         )
     );
 }

src/regionModels/regionModel/singleLayerRegion/singleLayerRegion.C

@@ -56,11 +56,11 @@ void Foam::regionModels::singleLayerRegion::constructMeshObjects()
                 time_.timeName(),
                 regionMesh(),
                 IOobject::READ_IF_PRESENT,
-                NO_WRITE
+                IOobject::NO_WRITE
             ),
             regionMesh(),
             dimensionedVector(dimless, Zero),
-            zeroGradientFvPatchField<vector>::typeName
+            fvPatchFieldBase::zeroGradientType()
         )
     );
 
@@ -75,11 +75,11 @@ void Foam::regionModels::singleLayerRegion::constructMeshObjects()
                 time_.timeName(),
                 regionMesh(),
                 IOobject::READ_IF_PRESENT,
-                NO_WRITE
+                IOobject::NO_WRITE
             ),
             regionMesh(),
             dimensionedScalar(dimArea, Zero),
-            zeroGradientFvPatchField<scalar>::typeName
+            fvPatchFieldBase::zeroGradientType()
         )
     );
 }

src/regionModels/regionModel/singleLayerRegion/singleLayerRegionTemplates.C

@@ -36,13 +36,14 @@ template<class Type>
 Foam::wordList
 Foam::regionModels::singleLayerRegion::mappedFieldAndInternalPatchTypes() const
 {
-    wordList bTypes(regionMesh().boundaryMesh().size());
+    wordList bTypes
+    (
+        regionMesh().boundaryMesh().size(),
+        fvPatchFieldBase::zeroGradientType()
+    );
 
-    bTypes = zeroGradientFvPatchField<Type>::typeName;
+    for (const label patchi : intCoupledPatchIDs_)
-
-    forAll(intCoupledPatchIDs_, i)
     {
-        const label patchi = intCoupledPatchIDs_[i];
         bTypes[patchi] = mappedFixedInternalValueFvPatchField<Type>::typeName;
     }
 
@@ -54,13 +55,14 @@ template<class Type>
 Foam::wordList
 Foam::regionModels::singleLayerRegion::mappedPushedFieldPatchTypes() const
 {
-    wordList bTypes(regionMesh().boundaryMesh().size());
+    wordList bTypes
+    (
+        regionMesh().boundaryMesh().size(),
+        fvPatchFieldBase::zeroGradientType()
+    );
 
-    bTypes = zeroGradientFvPatchField<Type>::typeName;
+    for (const label patchi : intCoupledPatchIDs_)
-
-    forAll(intCoupledPatchIDs_, i)
     {
-        const label patchi = intCoupledPatchIDs_[i];
         bTypes[patchi] =
             mappedFixedPushedInternalValueFvPatchField<Type>::typeName;
     }

src/regionModels/surfaceFilmModels/functionObjects/filmFlux/filmFlux.C

@@ -129,7 +129,7 @@ bool Foam::functionObjects::filmFlux::execute()
         ),
         filmMesh,
         dimensionedScalar(dimLength, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     );
 
     auto& heightc = height.ref();

src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C

@@ -468,7 +468,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimDensity, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     mu_
     (
@@ -482,7 +482,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimPressure*dimTime, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     sigma_
     (
@@ -496,7 +496,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimMass/sqr(dimTime), Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
 
     delta_
@@ -523,7 +523,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimless, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     U_
     (
@@ -548,7 +548,7 @@ kinematicSingleLayer::kinematicSingleLayer
             IOobject::NO_WRITE
         ),
         U_,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     Uw_
     (
@@ -561,7 +561,7 @@ kinematicSingleLayer::kinematicSingleLayer
             IOobject::NO_WRITE
         ),
         U_,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     deltaRho_
     (
@@ -575,7 +575,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(delta_.dimensions()*rho_.dimensions(), Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
 
     phi_
@@ -604,7 +604,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimMass, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     cloudMassTrans_
     (
@@ -618,7 +618,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimMass, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     cloudDiameterTrans_
     (
@@ -632,7 +632,7 @@ kinematicSingleLayer::kinematicSingleLayer
         ),
         regionMesh(),
         dimensionedScalar("minus1", dimLength, -1.0),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
 
     USp_

src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/curvatureSeparation/curvatureSeparation.C

@@ -204,11 +204,13 @@ tmp<scalarField> curvatureSeparation::calcCosAngle
                 "cosAngle",
                 mesh.time().timeName(),
                 mesh,
-                IOobject::NO_READ
+                IOobject::NO_READ,
+                IOobject::NO_WRITE,
+                IOobject::NO_REGISTER
             ),
             mesh,
             dimensionedScalar(dimless, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
         volCosAngle.primitiveFieldRef() = cosAngle;
         volCosAngle.correctBoundaryConditions();
@@ -346,11 +348,13 @@ void curvatureSeparation::correct
                 "Fnet",
                 mesh.time().timeName(),
                 mesh,
-                IOobject::NO_READ
+                IOobject::NO_READ,
+                IOobject::NO_WRITE,
+                IOobject::NO_REGISTER
             ),
             mesh,
             dimensionedScalar(dimForce, Zero),
-            zeroGradientFvPatchScalarField::typeName
+            fvPatchFieldBase::zeroGradientType()
         );
         volFnet.primitiveFieldRef() = Fnet;
         volFnet.correctBoundaryConditions();

src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C

@@ -85,7 +85,7 @@ standardRadiation::standardRadiation
         ),
         film.regionMesh(),
         dimensionedScalar(dimMass/pow3(dimTime), Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     beta_(coeffDict_.get<scalar>("beta")),
     kappaBar_(coeffDict_.get<scalar>("kappaBar"))

src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/waxSolventViscosity/waxSolventViscosity.C

@@ -126,7 +126,7 @@ waxSolventViscosity::waxSolventViscosity
         ),
         film.regionMesh(),
         dimensionedScalar(dimDynamicViscosity, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     muWaxModel_
     (
@@ -149,7 +149,7 @@ waxSolventViscosity::waxSolventViscosity
         ),
         film.regionMesh(),
         dimensionedScalar(dimDynamicViscosity, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     muSolventModel_
     (

src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/solidification/solidification.C

@@ -83,7 +83,7 @@ solidification::solidification
         ),
         film.regionMesh(),
         dimensionedScalar(dimMass, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     thickness_
     (
@@ -97,7 +97,7 @@ solidification::solidification
         ),
         film.regionMesh(),
         dimensionedScalar(dimLength, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     )
 {}
 

src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C

@@ -331,7 +331,7 @@ thermoSingleLayer::thermoSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimEnergy/dimMass/dimTemperature, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     kappa_
     (
@@ -345,7 +345,7 @@ thermoSingleLayer::thermoSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimEnergy/dimTime/dimLength/dimTemperature, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
 
     T_
@@ -371,7 +371,7 @@ thermoSingleLayer::thermoSingleLayer
             IOobject::NO_WRITE
         ),
         T_,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     Tw_
     (
@@ -384,7 +384,7 @@ thermoSingleLayer::thermoSingleLayer
             IOobject::NO_WRITE
         ),
         T_,
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
     hs_
     (
@@ -413,7 +413,7 @@ thermoSingleLayer::thermoSingleLayer
         ),
         regionMesh(),
         dimensionedScalar(dimEnergy, Zero),
-        zeroGradientFvPatchScalarField::typeName
+        fvPatchFieldBase::zeroGradientType()
     ),
 
     deltaWet_(coeffs_.get<scalar>("deltaWet")),