diff --git a/applications/solvers/combustion/dieselEngineFoam/createFields.H b/applications/solvers/combustion/dieselEngineFoam/createFields.H
index 6987608006..a75c6a7d07 100644
--- a/applications/solvers/combustion/dieselEngineFoam/createFields.H
+++ b/applications/solvers/combustion/dieselEngineFoam/createFields.H
@@ -88,7 +88,7 @@ volScalarField DpDt =
 
 multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-forAll (Y, i)
+forAll(Y, i)
 {
     fields.add(Y[i]);
 }
diff --git a/applications/solvers/combustion/reactingFoam/createFields.H b/applications/solvers/combustion/reactingFoam/createFields.H
index d92766db54..0b2238031c 100644
--- a/applications/solvers/combustion/reactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/createFields.H
@@ -77,7 +77,7 @@ volScalarField DpDt =
 
 multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-forAll (Y, i)
+forAll(Y, i)
 {
     fields.add(Y[i]);
 }
diff --git a/applications/solvers/combustion/rhoReactingFoam/createFields.H b/applications/solvers/combustion/rhoReactingFoam/createFields.H
index d44b4b9171..d2f6f5e4d4 100644
--- a/applications/solvers/combustion/rhoReactingFoam/createFields.H
+++ b/applications/solvers/combustion/rhoReactingFoam/createFields.H
@@ -77,7 +77,7 @@ volScalarField DpDt =
 
 multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-forAll (Y, i)
+forAll(Y, i)
 {
     fields.add(Y[i]);
 }
diff --git a/applications/solvers/lagrangian/coalChemistryFoam/createFields.H b/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
index 14d6f6f229..f79e48c17b 100644
--- a/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
@@ -28,7 +28,7 @@
 
     multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-    forAll (Y, i)
+    forAll(Y, i)
     {
         fields.add(Y[i]);
     }
diff --git a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/createFields.H b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/createFields.H
index 3d6c5500ea..68c0f8aee8 100644
--- a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/createFields.H
@@ -84,7 +84,7 @@
     Info<< "Creating multi-variate interpolation scheme\n" << endl;
     multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-    forAll (Y, i)
+    forAll(Y, i)
     {
         fields.add(Y[i]);
     }
diff --git a/applications/solvers/lagrangian/reactingParcelFoam/createFields.H b/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
index 359599a61f..e825c2cffb 100644
--- a/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
@@ -90,7 +90,7 @@
 
     multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
 
-    forAll (Y, i)
+    forAll(Y, i)
     {
         fields.add(Y[i]);
     }
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
index 7d6e41bf4c..95b35e9031 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/GidaspowErgunWenYu/GidaspowErgunWenYu.C
@@ -89,7 +89,7 @@ Foam::tmp<Foam::volScalarField> Foam::GidaspowErgunWenYu::K
     volScalarField& KWenYu = tKWenYu();
 
     // Ergun
-    forAll (beta, cellj)
+    forAll(beta, cellj)
     {
         if (beta[cellj] <= 0.8)
         {
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/SyamlalOBrien/SyamlalOBrien.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/SyamlalOBrien/SyamlalOBrien.C
index 52d847c854..497612462d 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/SyamlalOBrien/SyamlalOBrien.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/SyamlalOBrien/SyamlalOBrien.C
@@ -73,7 +73,7 @@ Foam::tmp<Foam::volScalarField> Foam::SyamlalOBrien::K
     volScalarField A = pow(beta, 4.14);
     volScalarField B = 0.8*pow(beta, 1.28);
 
-    forAll (beta, celli)
+    forAll(beta, celli)
     {
         if (beta[celli] > 0.85)
         {
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
index 3abba4d80b..3c992a4cbe 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/Schaeffer/SchaefferFrictionalStress.C
@@ -125,7 +125,7 @@ Foam::tmp<Foam::volScalarField> Foam::SchaefferFrictionalStress::muf
 
     volScalarField& muff = tmuf();
 
-    forAll (D, celli)
+    forAll(D, celli)
     {
         if (alpha[celli] > alphaMax.value()-5e-2)
         {
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/packingLimiter.H b/applications/solvers/multiphase/twoPhaseEulerFoam/packingLimiter.H
index 457e476b28..60c3bc8fa0 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/packingLimiter.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/packingLimiter.H
@@ -7,20 +7,20 @@
         labelListList neighbour = mesh.cellCells();
         scalarField cellVolumes = mesh.cellVolumes();
 
-        forAll (alphaEx, celli)
+        forAll(alphaEx, celli)
         {
             // Finding the labels of the neighbouring cells
             labelList neighbourCell = neighbour[celli];
 
             // Initializing neighbouring cells contribution
             scalar neighboursEx = 0.0;
-           
-            forAll (neighbourCell, cellj)
+
+            forAll(neighbourCell, cellj)
             {
                 labelList neighboursNeighbour = neighbour[neighbourCell[cellj]];
                 scalar neighboursNeighbourCellVolumes = 0.0;
-        
-                forAll (neighboursNeighbour, cellk)
+
+                forAll(neighboursNeighbour, cellk)
                 {
                     neighboursNeighbourCellVolumes +=
                         cellVolumes[neighboursNeighbour[cellk]];
diff --git a/applications/test/nearWallDist-wave/testWallDist2.C b/applications/test/nearWallDist-wave/testWallDist2.C
index 724702810c..a6a59d1780 100644
--- a/applications/test/nearWallDist-wave/testWallDist2.C
+++ b/applications/test/nearWallDist-wave/testWallDist2.C
@@ -73,7 +73,7 @@ int main(int argc, char *argv[])
 
     // Count walls
     label nWalls = 0;
-    forAll (mesh.boundary(), patchI)
+    forAll(mesh.boundary(), patchI)
     {
         const fvPatch& patch = mesh.boundary()[patchI];
 
@@ -87,13 +87,13 @@ int main(int argc, char *argv[])
     labelList changedFaces(nWalls);
 
     label nChangedFaces = 0;
-    forAll (mesh.boundary(), patchI)
+    forAll(mesh.boundary(), patchI)
     {
         const fvPatch& patch = mesh.boundary()[patchI];
 
         if (isA<wallFvPatch>(patch))
         {
-            forAll (patch.Cf(), patchFaceI)
+            forAll(patch.Cf(), patchFaceI)
             {
                 const polyPatch& polyPatch = mesh.boundaryMesh()[patchI];
 
@@ -170,7 +170,7 @@ int main(int argc, char *argv[])
         }
 
         // Copy boundary values
-        forAll (wallDistUncorrected.boundaryField(), patchI)
+        forAll(wallDistUncorrected.boundaryField(), patchI)
         {
             fvPatchScalarField& patchField =
                 wallDistUncorrected.boundaryField()[patchI];
diff --git a/applications/test/readCHEMKINIII/readCHEMKINIII.C b/applications/test/readCHEMKINIII/readCHEMKINIII.C
index 5449297127..01653ae705 100644
--- a/applications/test/readCHEMKINIII/readCHEMKINIII.C
+++ b/applications/test/readCHEMKINIII/readCHEMKINIII.C
@@ -72,7 +72,7 @@ int main(int argc, char *argv[])
 
         PtrList<gasReaction> testReactions(nReactions);
 
-        forAll (testReactions, i)
+        forAll(testReactions, i)
         {
             testReactions.set
             (
diff --git a/applications/test/router/processorRouter.C b/applications/test/router/processorRouter.C
index 5f68822cf3..3965d5b019 100644
--- a/applications/test/router/processorRouter.C
+++ b/applications/test/router/processorRouter.C
@@ -53,7 +53,7 @@ labelList procNeighbours(const polyMesh& mesh)
 
     label nNeighbours = 0;
 
-    forAll (mesh.boundaryMesh(), patchI)
+    forAll(mesh.boundaryMesh(), patchI)
     {
         if (isA<processorPolyPatch>(mesh.boundaryMesh()[patchI]))
         {
@@ -65,7 +65,7 @@ labelList procNeighbours(const polyMesh& mesh)
 
     nNeighbours = 0;
 
-    forAll (mesh.boundaryMesh(), patchI)
+    forAll(mesh.boundaryMesh(), patchI)
     {
         if (isA<processorPolyPatch>(mesh.boundaryMesh()[patchI]))
         {
diff --git a/applications/utilities/mesh/advanced/refinementLevel/refinementLevel.C b/applications/utilities/mesh/advanced/refinementLevel/refinementLevel.C
index 38842e7394..d96afd1619 100644
--- a/applications/utilities/mesh/advanced/refinementLevel/refinementLevel.C
+++ b/applications/utilities/mesh/advanced/refinementLevel/refinementLevel.C
@@ -216,7 +216,7 @@ int main(int argc, char *argv[])
 
     List<polyPatch*> p(patches.size());
 
-    forAll (p, patchI)
+    forAll(p, patchI)
     {
         p[patchI] = patches[patchI].clone(fMesh.boundaryMesh()).ptr();
     }
diff --git a/applications/utilities/mesh/conversion/cfx4ToFoam/cfx4ToFoam.C b/applications/utilities/mesh/conversion/cfx4ToFoam/cfx4ToFoam.C
index e7a5a5ba56..7ddfeb5fa2 100644
--- a/applications/utilities/mesh/conversion/cfx4ToFoam/cfx4ToFoam.C
+++ b/applications/utilities/mesh/conversion/cfx4ToFoam/cfx4ToFoam.C
@@ -81,7 +81,7 @@ int main(int argc, char *argv[])
         word blockName;
         label nx, ny, nz;
 
-        forAll (blocks, blockI)
+        forAll(blocks, blockI)
         {
             cfxFile >> blockName;
             cfxFile >> nx >> ny >> nz;
@@ -101,7 +101,7 @@ int main(int argc, char *argv[])
     {
         label no, blkNo, patchLabel;
 
-        forAll (cfxPatchTypes, patchI)
+        forAll(cfxPatchTypes, patchI)
         {
             // Grab patch type and name
             cfxFile >> cfxPatchTypes[patchI] >> cfxPatchNames[patchI] >> no;
@@ -110,7 +110,7 @@ int main(int argc, char *argv[])
             patchRanges[patchI].setSize(6);
             labelList& curRange = patchRanges[patchI];
 
-            forAll (curRange, rI)
+            forAll(curRange, rI)
             {
                 cfxFile >> curRange[rI];
             }
@@ -147,7 +147,7 @@ int main(int argc, char *argv[])
 
     Info<< "Reading block points" << endl;
 
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         Info<< "block " << blockI << " is a ";
         blocks[blockI].readPoints(cfxFile);
@@ -176,7 +176,7 @@ int main(int argc, char *argv[])
 
     faceListList rawPatches(npatch);
 
-    forAll (rawPatches, patchI)
+    forAll(rawPatches, patchI)
     {
         const word& patchType = cfxPatchTypes[patchI];
 
@@ -214,7 +214,7 @@ int main(int argc, char *argv[])
     // For efficiency, create merge pairs in the first pass
     labelListListList glueMergePairs(glueMasterPatches.size());
 
-    forAll (glueMasterPatches, glueI)
+    forAll(glueMasterPatches, glueI)
     {
         const label masterPatch = glueMasterPatches[glueI];
         const label slavePatch = glueSlavePatches[glueI];
@@ -247,14 +247,14 @@ int main(int argc, char *argv[])
 
         scalar sqrMergeTol = GREAT;
 
-        forAll (blockPFaces, blockPFaceLabel)
+        forAll(blockPFaces, blockPFaceLabel)
         {
             const labelList& blockPFacePoints =
                 blockPFaces[blockPFaceLabel];
 
-            forAll (blockPFacePoints, blockPFacePointI)
+            forAll(blockPFacePoints, blockPFacePointI)
             {
-                forAll (blockPFacePoints, blockPFacePointI2)
+                forAll(blockPFacePoints, blockPFacePointI2)
                 {
                     if (blockPFacePointI != blockPFacePointI2)
                     {
@@ -281,7 +281,7 @@ int main(int argc, char *argv[])
 
         // N-squared point search over all points of all faces of
         // master block over all point of all faces of slave block
-        forAll (blockPFaces, blockPFaceLabel)
+        forAll(blockPFaces, blockPFaceLabel)
         {
             const labelList& blockPFacePoints =
                 blockPFaces[blockPFaceLabel];
@@ -289,16 +289,16 @@ int main(int argc, char *argv[])
             labelList& cp = curPairs[blockPFaceLabel];
             cp.setSize(blockPFacePoints.size());
 
-        forAll (blockPFacePoints, blockPFacePointI)
+        forAll(blockPFacePoints, blockPFacePointI)
         {
             found = false;
 
-            forAll (blockNFaces, blockNFaceLabel)
+            forAll(blockNFaces, blockNFaceLabel)
             {
                 const labelList& blockNFacePoints =
                     blockNFaces[blockNFaceLabel];
 
-            forAll (blockNFacePoints, blockNFacePointI)
+            forAll(blockNFacePoints, blockNFacePointI)
             {
                 if
                 (
@@ -360,7 +360,7 @@ int main(int argc, char *argv[])
         changedPointMerge = false;
         nPasses++;
 
-        forAll (glueMasterPatches, glueI)
+        forAll(glueMasterPatches, glueI)
         {
             const label masterPatch = glueMasterPatches[glueI];
             const label slavePatch = glueSlavePatches[glueI];
@@ -372,14 +372,14 @@ int main(int argc, char *argv[])
 
             const labelListList& curPairs = glueMergePairs[glueI];
 
-            forAll (blockPFaces, blockPFaceLabel)
+            forAll(blockPFaces, blockPFaceLabel)
             {
                 const labelList& blockPFacePoints =
                     blockPFaces[blockPFaceLabel];
 
                 const labelList& cp = curPairs[blockPFaceLabel];
 
-                forAll (cp, blockPFacePointI)
+                forAll(cp, blockPFacePointI)
                 {
                     label PpointLabel =
                         blockPFacePoints[blockPFacePointI]
@@ -421,7 +421,7 @@ int main(int argc, char *argv[])
     }
 
 
-    forAll (glueMasterPatches, glueI)
+    forAll(glueMasterPatches, glueI)
     {
         const label masterPatch = glueMasterPatches[glueI];
         const label slavePatch = glueSlavePatches[glueI];
@@ -433,12 +433,12 @@ int main(int argc, char *argv[])
         const faceList& blockNFaces = rawPatches[slavePatch];
 
 
-        forAll (blockPFaces, blockPFaceLabel)
+        forAll(blockPFaces, blockPFaceLabel)
         {
             const labelList& blockPFacePoints
                 = blockPFaces[blockPFaceLabel];
 
-            forAll (blockPFacePoints, blockPFacePointI)
+            forAll(blockPFacePoints, blockPFacePointI)
             {
                 label PpointLabel =
                     blockPFacePoints[blockPFacePointI]
@@ -455,12 +455,12 @@ int main(int argc, char *argv[])
             }
         }
 
-        forAll (blockNFaces, blockNFaceLabel)
+        forAll(blockNFaces, blockNFaceLabel)
         {
             const labelList& blockNFacePoints
                 = blockNFaces[blockNFaceLabel];
 
-            forAll (blockNFacePoints, blockNFacePointI)
+            forAll(blockNFacePoints, blockNFacePointI)
             {
                 label NpointLabel =
                     blockNFacePoints[blockNFacePointI]
@@ -483,7 +483,7 @@ int main(int argc, char *argv[])
     // given old point label
     label nNewPoints = 0;
 
-    forAll (pointMergeList, pointLabel)
+    forAll(pointMergeList, pointLabel)
     {
         if (pointMergeList[pointLabel] > pointLabel)
         {
@@ -513,11 +513,11 @@ int main(int argc, char *argv[])
 
     pointField points(nMeshPoints);
 
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         const pointField& blockPoints = blocks[blockI].points();
 
-        forAll (blockPoints, blockPointLabel)
+        forAll(blockPoints, blockPointLabel)
         {
             points
             [
@@ -544,15 +544,15 @@ int main(int argc, char *argv[])
 
     label nCreatedCells = 0;
 
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         labelListList curBlockCells = blocks[blockI].blockCells();
 
-        forAll (curBlockCells, blockCellI)
+        forAll(curBlockCells, blockCellI)
         {
             labelList cellPoints(curBlockCells[blockCellI].size());
 
-            forAll (cellPoints, pointI)
+            forAll(cellPoints, pointI)
             {
                 cellPoints[pointI] =
                     pointMergeList
@@ -579,7 +579,7 @@ int main(int argc, char *argv[])
 
     label nCreatedPatches = 0;
 
-    forAll (rawPatches, patchI)
+    forAll(rawPatches, patchI)
     {
         if (rawPatches[patchI].size() && cfxPatchTypes[patchI] != "BLKBDY")
         {
@@ -610,14 +610,14 @@ int main(int argc, char *argv[])
                 label oldSize = renumberedPatch.size();
                 renumberedPatch.setSize(oldSize + curRawPatch.size());
 
-                forAll (curRawPatch, faceI)
+                forAll(curRawPatch, faceI)
                 {
                     const face& oldFace = curRawPatch[faceI];
 
                     face& newFace = renumberedPatch[oldSize + faceI];
                     newFace.setSize(oldFace.size());
 
-                    forAll (oldFace, pointI)
+                    forAll(oldFace, pointI)
                     {
                         newFace[pointI] =
                             pointMergeList
@@ -634,14 +634,14 @@ int main(int argc, char *argv[])
                 faceList& renumberedPatch = boundary[nCreatedPatches];
                 renumberedPatch.setSize(curRawPatch.size());
 
-                forAll (curRawPatch, faceI)
+                forAll(curRawPatch, faceI)
                 {
                     const face& oldFace = curRawPatch[faceI];
 
                     face& newFace = renumberedPatch[faceI];
                     newFace.setSize(oldFace.size());
 
-                    forAll (oldFace, pointI)
+                    forAll(oldFace, pointI)
                     {
                         newFace[pointI] =
                             pointMergeList
diff --git a/applications/utilities/mesh/conversion/cfx4ToFoam/hexBlock.C b/applications/utilities/mesh/conversion/cfx4ToFoam/hexBlock.C
index bf16792211..df0d7d97a4 100644
--- a/applications/utilities/mesh/conversion/cfx4ToFoam/hexBlock.C
+++ b/applications/utilities/mesh/conversion/cfx4ToFoam/hexBlock.C
@@ -58,7 +58,7 @@ hexBlock::hexBlock(const label nx, const label ny, const label nz)
 
 void hexBlock::readPoints(Istream& is)
 {
-    forAll (points_, i)
+    forAll(points_, i)
     {
         is >> points_[i].x() >> points_[i].y() >> points_[i].z();
     }
@@ -376,7 +376,7 @@ faceList hexBlock::patchFaces(const label direc, const labelList& range) const
     else if (blockHandedness_ == left)
     {
         // turn all faces inside out
-        forAll (result, faceI)
+        forAll(result, faceI)
         {
             result[faceI] = result[faceI].reverseFace();
         }
diff --git a/applications/utilities/mesh/conversion/fluent3DMeshToFoam/fluent3DMeshToFoam.L b/applications/utilities/mesh/conversion/fluent3DMeshToFoam/fluent3DMeshToFoam.L
index 4cb7406fec..333188c072 100644
--- a/applications/utilities/mesh/conversion/fluent3DMeshToFoam/fluent3DMeshToFoam.L
+++ b/applications/utilities/mesh/conversion/fluent3DMeshToFoam/fluent3DMeshToFoam.L
@@ -473,7 +473,7 @@ endOfSection               {space}")"{space}
         // set size of label list
         curFaceLabels.setSize(strtol(YYText(), &endPtr, 16));
 
-        forAll (curFaceLabels, i)
+        forAll(curFaceLabels, i)
         {
             curFaceLabels[i] = strtol(endPtr, &endPtr, 16) - 1;
         }
@@ -837,7 +837,7 @@ int main(int argc, char *argv[])
     // faces and move to neighbour
 
     boolList fm(faces.size(), false);
-    forAll (faces, facei)
+    forAll(faces, facei)
     {
         if
         (
diff --git a/applications/utilities/mesh/conversion/fluentMeshToFoam/create3DCellShape.C b/applications/utilities/mesh/conversion/fluentMeshToFoam/create3DCellShape.C
index 2645fe7509..0c23f976ac 100644
--- a/applications/utilities/mesh/conversion/fluentMeshToFoam/create3DCellShape.C
+++ b/applications/utilities/mesh/conversion/fluentMeshToFoam/create3DCellShape.C
@@ -155,7 +155,7 @@ cellShape create3DCellShape
 
     bool found = false;
 
-    forAll (localFaces, meshFaceI)
+    forAll(localFaces, meshFaceI)
     {
         if (localFaces[meshFaceI].size() == firstModelFace.size())
         {
@@ -166,7 +166,7 @@ cellShape create3DCellShape
 
             meshFaceUsed[meshFaceI] = true;
 
-            forAll (curMeshFace, pointI)
+            forAll(curMeshFace, pointI)
             {
                 pointLabels[firstModelFace[pointI]] = curMeshFace[pointI];
             }
@@ -199,7 +199,7 @@ cellShape create3DCellShape
         found = false;
 
         // Loop through mesh faces until a match is found
-        forAll (localFaces, meshFaceI)
+        forAll(localFaces, meshFaceI)
         {
             if
             (
@@ -220,7 +220,7 @@ cellShape create3DCellShape
                     // try matching the face
                     label nMatchedLabels = 0;
 
-                    forAll (meshFaceLabels, pointI)
+                    forAll(meshFaceLabels, pointI)
                     {
                         if
                         (
@@ -241,7 +241,7 @@ cellShape create3DCellShape
                     if (found)
                     {
                         // match found. Insert mesh face
-                        forAll (meshFaceLabels, pointI)
+                        forAll(meshFaceLabels, pointI)
                         {
                             pointLabels[curModelFace[pointI]] =
                                 meshFaceLabels[pointI];
diff --git a/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedQuadCellShape.C b/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedQuadCellShape.C
index 9a4ff9a9fa..2fccc72b7d 100644
--- a/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedQuadCellShape.C
+++ b/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedQuadCellShape.C
@@ -73,7 +73,7 @@ cellShape extrudedQuadCellShape
     // make a list of outward-pointing faces
     labelListList localFaces(4);
 
-    forAll (faceLabels, faceI)
+    forAll(faceLabels, faceI)
     {
         const label curFaceLabel = faceLabels[faceI];
 
@@ -100,7 +100,7 @@ cellShape extrudedQuadCellShape
         {
             // Reverse the face.  Note: it is necessary to reverse by
             // hand to preserve connectivity of a 2-D mesh.
-            // 
+            //
             localFaces[faceI].setSize(curFace.size());
 
             forAllReverse(curFace, i)
diff --git a/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedTriangleCellShape.C b/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedTriangleCellShape.C
index 401ea6225c..8b2af5980f 100644
--- a/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedTriangleCellShape.C
+++ b/applications/utilities/mesh/conversion/fluentMeshToFoam/extrudedTriangleCellShape.C
@@ -75,7 +75,7 @@ cellShape extrudedTriangleCellShape
     // make a list of outward-pointing faces
     labelListList localFaces(3);
 
-    forAll (faceLabels, faceI)
+    forAll(faceLabels, faceI)
     {
         const label curFaceLabel = faceLabels[faceI];
 
@@ -102,7 +102,7 @@ cellShape extrudedTriangleCellShape
         {
             // Reverse the face.  Note: it is necessary to reverse by
             // hand to preserve connectivity of a 2-D mesh.
-            // 
+            //
             localFaces[faceI].setSize(curFace.size());
 
             forAllReverse(curFace, i)
diff --git a/applications/utilities/mesh/conversion/fluentMeshToFoam/fluentMeshToFoam.L b/applications/utilities/mesh/conversion/fluentMeshToFoam/fluentMeshToFoam.L
index 63ac764176..79a6d75633 100644
--- a/applications/utilities/mesh/conversion/fluentMeshToFoam/fluentMeshToFoam.L
+++ b/applications/utilities/mesh/conversion/fluentMeshToFoam/fluentMeshToFoam.L
@@ -490,7 +490,7 @@ endOfSection               {space}")"{space}
         // set size of label list
         curFaceLabels.setSize(readLabel(mixedFaceStream));
 
-        forAll (curFaceLabels, i)
+        forAll(curFaceLabels, i)
         {
             curFaceLabels[i] = readHexLabel(mixedFaceStream) - 1;
         }
@@ -511,7 +511,7 @@ endOfSection               {space}")"{space}
         // for edge is 2, for triangle is 3 and for quad is 4
         curFaceLabels.setSize(faceGroupElementType);
 
-        forAll (curFaceLabels, i)
+        forAll(curFaceLabels, i)
         {
             curFaceLabels[i] = readHexLabel(mixedFaceStream) - 1;
         }
@@ -915,14 +915,14 @@ int main(int argc, char *argv[])
 
     labelList nFacesInCell(nCells, 0);
 
-    forAll (cellFaces, celli)
+    forAll(cellFaces, celli)
     {
         cellFaces[celli].setSize(fluentModelNFaces[fluentCellModelID[celli] ]);
     }
 
     // fill in owner and neighbour
 
-    forAll (owner, faceI)
+    forAll(owner, faceI)
     {
         if (owner[faceI] > -1)
         {
@@ -933,7 +933,7 @@ int main(int argc, char *argv[])
         }
     }
 
-    forAll (neighbour, faceI)
+    forAll(neighbour, faceI)
     {
         if (neighbour[faceI] > -1)
         {
@@ -976,7 +976,7 @@ int main(int argc, char *argv[])
         // points given by Fluent need to represent the FRONT plane of the
         // geometry. Therefore, the extrusion will be in -z direction
         //
-        forAll (oldPoints, pointI)
+        forAll(oldPoints, pointI)
         {
             points[nNewPoints] = oldPoints[pointI];
 
@@ -985,7 +985,7 @@ int main(int argc, char *argv[])
             nNewPoints++;
         }
 
-        forAll (oldPoints, pointI)
+        forAll(oldPoints, pointI)
         {
             points[nNewPoints] = oldPoints[pointI];
 
@@ -1000,7 +1000,7 @@ int main(int argc, char *argv[])
         // Set the number of empty faces
         frontAndBackFaces.setSize(2*nCells);
 
-        forAll (fluentCellModelID, celli)
+        forAll(fluentCellModelID, celli)
         {
             switch (fluentCellModelID[celli])
             {
@@ -1053,7 +1053,7 @@ int main(int argc, char *argv[])
         }
 
         // Create new faces
-        forAll (faces, faceI)
+        forAll(faces, faceI)
         {
 
             if (faces[faceI].size() != 2)
@@ -1078,7 +1078,7 @@ int main(int argc, char *argv[])
     {
         // 3-D shape recognition
         Info<< "Creating shapes for 3-D cells"<< endl;
-        forAll (fluentCellModelID, celli)
+        forAll(fluentCellModelID, celli)
         {
             if
             (
@@ -1116,7 +1116,7 @@ int main(int argc, char *argv[])
     // area vector points into the domain. Turn them round before making patches
     // for Foam compatibility
 
-    forAll (faces, faceI)
+    forAll(faces, faceI)
     {
         if (owner[faceI] == -1)
         {
@@ -1232,7 +1232,7 @@ int main(int argc, char *argv[])
 
         faceList patchFaces(faceGroupEndIndexIter() - faceLabel);
 
-        forAll (patchFaces, faceI)
+        forAll(patchFaces, faceI)
         {
             if
             (
diff --git a/applications/utilities/mesh/conversion/foamMeshToFluent/fluentFvMesh.C b/applications/utilities/mesh/conversion/foamMeshToFluent/fluentFvMesh.C
index 228194e4c3..00b6874b5b 100644
--- a/applications/utilities/mesh/conversion/foamMeshToFluent/fluentFvMesh.C
+++ b/applications/utilities/mesh/conversion/foamMeshToFluent/fluentFvMesh.C
@@ -113,7 +113,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
 
     const pointField& p = points();
 
-    forAll (p, pointI)
+    forAll(p, pointI)
     {
         fluentMeshFile
             << "    "
@@ -135,7 +135,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
         << "(13 (2 1 "
         << own.size() << " 2 0)" << std::endl << "(" << std::endl;
 
-    forAll (own, faceI)
+    forAll(own, faceI)
     {
         const labelList& l = fcs[faceI];
 
@@ -143,7 +143,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
 
         fluentMeshFile << l.size() << " ";
 
-        forAll (l, lI)
+        forAll(l, lI)
         {
             fluentMeshFile << l[lI] + 1 << " ";
         }
@@ -157,7 +157,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
     label nWrittenFaces = own.size();
 
     // Writing boundary faces
-    forAll (boundary(), patchI)
+    forAll(boundary(), patchI)
     {
         const unallocFaceList& patchFaces = boundaryMesh()[patchI];
 
@@ -190,7 +190,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
         fluentMeshFile
             <<" 0)" << std::endl << "(" << std::endl;
 
-        forAll (patchFaces, faceI)
+        forAll(patchFaces, faceI)
         {
             const labelList& l = patchFaces[faceI];
 
@@ -225,7 +225,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
 
     bool hasWarned = false;
 
-    forAll (cells, cellI)
+    forAll(cells, cellI)
     {
         if (cells[cellI].model() == tet)
         {
@@ -270,7 +270,7 @@ void Foam::fluentFvMesh::writeFluentMesh() const
     fluentMeshFile << "(39 (2 interior interior-1)())" << std::endl;
 
     // Writing boundary patch types
-    forAll (boundary(), patchI)
+    forAll(boundary(), patchI)
     {
         fluentMeshFile
             << "(39 (" << patchI + 10 << " ";
diff --git a/applications/utilities/mesh/conversion/gambitToFoam/gambitToFoam.L b/applications/utilities/mesh/conversion/gambitToFoam/gambitToFoam.L
index 47353f8e42..b17e1f6c99 100644
--- a/applications/utilities/mesh/conversion/gambitToFoam/gambitToFoam.L
+++ b/applications/utilities/mesh/conversion/gambitToFoam/gambitToFoam.L
@@ -583,7 +583,7 @@ mtype                 {space}"MTYPE:"{space}
         {
             scalarList patchFaceValues(nValuesForPatchFaces);
 
-            forAll (patchFaceValues, fI)
+            forAll(patchFaceValues, fI)
             {
                 patchFaceValues[fI] = readScalar(patchFacesStream);
             }
@@ -673,7 +673,7 @@ int main(int argc, char *argv[])
     // make a point mapping array
     label maxPointIndex = 0;
 
-    forAll (pointMap, pointI)
+    forAll(pointMap, pointI)
     {
         if (pointMap[pointI] > maxPointIndex)
         {
@@ -684,7 +684,7 @@ int main(int argc, char *argv[])
 
     labelList pointLookup(maxPointIndex + 1, -1);
 
-    forAll (pointMap, pointI)
+    forAll(pointMap, pointI)
     {
         pointLookup[pointMap[pointI] ] = pointI;
     }
@@ -692,7 +692,7 @@ int main(int argc, char *argv[])
     // make a cell mapping array
     label maxCellIndex = 0;
 
-    forAll (cellMap, cellI)
+    forAll(cellMap, cellI)
     {
         if (cellMap[cellI] > maxCellIndex)
         {
@@ -702,7 +702,7 @@ int main(int argc, char *argv[])
 
     labelList cellLookup(maxCellIndex + 1);
 
-    forAll (cellMap, cellI)
+    forAll(cellMap, cellI)
     {
         cellLookup[cellMap[cellI] ] = cellI;
     }
@@ -719,7 +719,7 @@ int main(int argc, char *argv[])
 
     cellShapeList cells(cellLabels.size());
 
-    forAll (cellTypes, cellI)
+    forAll(cellTypes, cellI)
     {
         const labelList& curCellLabels = cellLabels[cellI];
 
@@ -798,7 +798,7 @@ int main(int argc, char *argv[])
         faceList& patchFaces = boundary[patchI];
         patchFaces.setSize(curCells.size());
 
-        forAll (curCells, faceI)
+        forAll(curCells, faceI)
         {
             patchFaces[faceI] =
                 cells[cellLookup[curCells[faceI] ] ].faces()
diff --git a/applications/utilities/mesh/conversion/kivaToFoam/checkPatch.H b/applications/utilities/mesh/conversion/kivaToFoam/checkPatch.H
index 451c6865cc..dd4a6f2d0b 100644
--- a/applications/utilities/mesh/conversion/kivaToFoam/checkPatch.H
+++ b/applications/utilities/mesh/conversion/kivaToFoam/checkPatch.H
@@ -68,14 +68,14 @@
         quadFace[3] = pointMap[quadFace[3]];
 
         label ti = 0;
-        forAll (quadFace, i)
+        forAll(quadFace, i)
         {
             if (quadFace[i] != quadFace[(i+1)%4])
             {
                 triFace[(ti++)%3] = quadFace[i];
             }
         }
-        
+
         if (ti == 4)
         {
             pFaces[bcIndex][regionIndex].append(quadFace);
diff --git a/applications/utilities/mesh/conversion/kivaToFoam/readKivaGrid.H b/applications/utilities/mesh/conversion/kivaToFoam/readKivaGrid.H
index 92c7b7c355..71b82cfb07 100644
--- a/applications/utilities/mesh/conversion/kivaToFoam/readKivaGrid.H
+++ b/applications/utilities/mesh/conversion/kivaToFoam/readKivaGrid.H
@@ -109,7 +109,7 @@ label activeCells = 0;
 // Create and set the collocated point collapse map
 labelList pointMap(nPoints);
 
-forAll (pointMap, i)
+forAll(pointMap, i)
 {
     pointMap[i] = i;
 }
@@ -132,7 +132,7 @@ for (label i=0; i<nPoints; i++)
 
         edgeList edges = cellShapes[activeCells].edges();
 
-        forAll (edges, ei)
+        forAll(edges, ei)
         {
             if (edges[ei].mag(points) < SMALL)
             {
@@ -167,11 +167,11 @@ cellZoning.setSize(activeCells);
 
 // Map collocated points to refer to the same point and collapse cell shape
 // to the corresponding hex-degenerate.
-forAll (cellShapes, celli)
+forAll(cellShapes, celli)
 {
     cellShape& cs = cellShapes[celli];
 
-    forAll (cs, i)
+    forAll(cs, i)
     {
         cs[i] = pointMap[cs[i]];
     }
@@ -336,7 +336,7 @@ if
     {
         const face& pf = iter();
 
-        forAll (pf, pfi)
+        forAll(pf, pfi)
         {
             minz = min(minz, points[pf[pfi]].z());
         }
@@ -356,7 +356,7 @@ if
         const face& pf = iter();
 
         scalar minfz = GREAT;
-        forAll (pf, pfi)
+        forAll(pf, pfi)
         {
             minfz = min(minfz, points[pf[pfi]].z());
         }
@@ -390,7 +390,7 @@ if
         const face& pf = iter();
 
         scalar minfz = GREAT;
-        forAll (pf, pfi)
+        forAll(pf, pfi)
         {
             minfz = min(minfz, points[pf[pfi]].z());
         }
@@ -417,7 +417,7 @@ if
 label nPatches = 0;
 for (int bci=0; bci<nBCs; bci++)
 {
-    forAll (pFaces[bci], rgi)
+    forAll(pFaces[bci], rgi)
     {
         if (pFaces[bci][rgi].size())
         {
@@ -485,7 +485,7 @@ label nAddedPatches = 0;
 
 for (int bci=0; bci<nBCs; bci++)
 {
-    forAll (pFaces[bci], rgi)
+    forAll(pFaces[bci], rgi)
     {
         if (pFaces[bci][rgi].size())
         {
@@ -510,9 +510,9 @@ for (int bci=0; bci<nBCs; bci++)
 labelList pointLabels(nPoints, -1);
 
 // Scan cells for used points
-forAll (cellShapes, celli)
+forAll(cellShapes, celli)
 {
-    forAll (cellShapes[celli], i)
+    forAll(cellShapes[celli], i)
     {
         pointLabels[cellShapes[celli][i]] = 1;
     }
@@ -520,7 +520,7 @@ forAll (cellShapes, celli)
 
 // Create addressing for used points and pack points array
 label newPointi = 0;
-forAll (pointLabels, pointi)
+forAll(pointLabels, pointi)
 {
     if (pointLabels[pointi] != -1)
     {
@@ -531,24 +531,24 @@ forAll (pointLabels, pointi)
 points.setSize(newPointi);
 
 // Reset cell point labels
-forAll (cellShapes, celli)
+forAll(cellShapes, celli)
 {
     cellShape& cs = cellShapes[celli];
 
-    forAll (cs, i)
+    forAll(cs, i)
     {
         cs[i] = pointLabels[cs[i]];
     }
 }
 
 // Reset boundary-face point labels
-forAll (boundary, patchi)
+forAll(boundary, patchi)
 {
-    forAll (boundary[patchi], facei)
+    forAll(boundary[patchi], facei)
     {
         face& f = boundary[patchi][facei];
 
-        forAll (f, i)
+        forAll(f, i)
         {
             f[i] = pointLabels[f[i]];
         }
diff --git a/applications/utilities/mesh/conversion/plot3dToFoam/hexBlock.C b/applications/utilities/mesh/conversion/plot3dToFoam/hexBlock.C
index 3ed248a86e..2a36f86a22 100644
--- a/applications/utilities/mesh/conversion/plot3dToFoam/hexBlock.C
+++ b/applications/utilities/mesh/conversion/plot3dToFoam/hexBlock.C
@@ -471,7 +471,7 @@ faceList hexBlock::patchFaces(const label direc, const labelList& range) const
     else if (blockHandedness_ == left)
     {
         // turn all faces inside out
-        forAll (result, faceI)
+        forAll(result, faceI)
         {
             result[faceI] = result[faceI].reverseFace();
         }
diff --git a/applications/utilities/mesh/conversion/plot3dToFoam/plot3dToFoam.C b/applications/utilities/mesh/conversion/plot3dToFoam/plot3dToFoam.C
index 8ad25caea6..0f2c18732b 100644
--- a/applications/utilities/mesh/conversion/plot3dToFoam/plot3dToFoam.C
+++ b/applications/utilities/mesh/conversion/plot3dToFoam/plot3dToFoam.C
@@ -108,7 +108,7 @@ int main(int argc, char *argv[])
     {
         label nx, ny, nz;
 
-        forAll (blocks, blockI)
+        forAll(blocks, blockI)
         {
             if (twoDThickness > 0)
             {
@@ -132,7 +132,7 @@ int main(int argc, char *argv[])
     label sumPoints(0);
     label nMeshCells(0);
 
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         Info<< "block " << blockI << ":" << nl;
         blocks[blockI].readPoints(readBlank, twoDThickness, plot3dFile);
@@ -144,11 +144,11 @@ int main(int argc, char *argv[])
     pointField points(sumPoints);
     labelList blockOffsets(blocks.size());
     sumPoints = 0;
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         const pointField& blockPoints = blocks[blockI].points();
         blockOffsets[blockI] = sumPoints;
-        forAll (blockPoints, i)
+        forAll(blockPoints, i)
         {
             points[sumPoints++] = blockPoints[i];
         }
@@ -186,15 +186,15 @@ int main(int argc, char *argv[])
 
     label nCreatedCells = 0;
 
-    forAll (blocks, blockI)
+    forAll(blocks, blockI)
     {
         labelListList curBlockCells = blocks[blockI].blockCells();
 
-        forAll (curBlockCells, blockCellI)
+        forAll(curBlockCells, blockCellI)
         {
             labelList cellPoints(curBlockCells[blockCellI].size());
 
-            forAll (cellPoints, pointI)
+            forAll(cellPoints, pointI)
             {
                 cellPoints[pointI] =
                     oldToNew
diff --git a/applications/utilities/mesh/conversion/sammToFoam/calcPointCells.C b/applications/utilities/mesh/conversion/sammToFoam/calcPointCells.C
index c94d732234..c59fd43eda 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/calcPointCells.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/calcPointCells.C
@@ -67,14 +67,14 @@ void sammMesh::calcPointCells() const
     // do not come from the cell shape, but from the slaves of the coupled
     // match. It is also adventageous to remove the duplicates from the
     // point-cell addressing, because this removes a lot of waste later.
-    // 
+    //
 
     // For each cell
     forAll(cellShapes_, cellI)
     {
         const faceList& faces = cellFaces_[cellI];
 
-        forAll (faces, i)
+        forAll(faces, i)
         {
             // For each vertex
             const labelList& labels = faces[i];
diff --git a/applications/utilities/mesh/conversion/sammToFoam/createBoundaryFaces.C b/applications/utilities/mesh/conversion/sammToFoam/createBoundaryFaces.C
index dd4e7b944e..304d047a66 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/createBoundaryFaces.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/createBoundaryFaces.C
@@ -45,11 +45,11 @@ bool sammMesh::sammEqualFace
     // considered equal if three of the vertices are the same.
     label nEqual = 0;
 
-    forAll (cellFace, cellFaceLabelI)
+    forAll(cellFace, cellFaceLabelI)
     {
         const label curCellFaceLabel = cellFace[cellFaceLabelI];
 
-        forAll (boundaryFace, bouFaceLabelI)
+        forAll(boundaryFace, bouFaceLabelI)
         {
             if (boundaryFace[bouFaceLabelI] == curCellFaceLabel)
             {
diff --git a/applications/utilities/mesh/conversion/sammToFoam/createPolyBoundary.C b/applications/utilities/mesh/conversion/sammToFoam/createPolyBoundary.C
index f3fc9a6e1d..0309ac90b3 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/createPolyBoundary.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/createPolyBoundary.C
@@ -42,13 +42,13 @@ void sammMesh::createPolyBoundary()
 
     const labelListList& PointCells = pointCells();
 
-    forAll (boundary_, patchI)
+    forAll(boundary_, patchI)
     {
         const faceList& curShapePatch = boundary_[patchI];
 
         polyBoundaryPatchStartIndices_[patchI] = nCreatedFaces;
 
-        forAll (curShapePatch, faceI)
+        forAll(curShapePatch, faceI)
         {
             bool found = false;
 
@@ -120,7 +120,7 @@ List<polyPatch* > sammMesh::polyBoundaryPatches(const polyMesh& pMesh)
 {
     List<polyPatch* > p(boundary_.size());
 
-    forAll (boundary_, patchI)
+    forAll(boundary_, patchI)
     {
         const faceList& curShapePatch = boundary_[patchI];
 
diff --git a/applications/utilities/mesh/conversion/sammToFoam/createPolyCells.C b/applications/utilities/mesh/conversion/sammToFoam/createPolyCells.C
index 80dd2e7541..276626e099 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/createPolyCells.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/createPolyCells.C
@@ -41,13 +41,13 @@ void sammMesh::createPolyCells()
 
     label maxFaces = 0;
 
-    forAll (cellPolys_, cellI)
+    forAll(cellPolys_, cellI)
     {
         cell& curCell = cellPolys_[cellI];
 
         curCell.setSize(cellFaces_[cellI].size());
 
-        forAll (curCell, fI)
+        forAll(curCell, fI)
         {
             curCell[fI] = -1;
         }
@@ -145,7 +145,7 @@ void sammMesh::createPolyCells()
             label nextNei = -1;
             label minNei = cellPolys_.size();
 
-            forAll (neiCells, ncI)
+            forAll(neiCells, ncI)
             {
                 if (neiCells[ncI] > -1 && neiCells[ncI] < minNei)
                 {
diff --git a/applications/utilities/mesh/conversion/sammToFoam/fixCollapsedEdges.C b/applications/utilities/mesh/conversion/sammToFoam/fixCollapsedEdges.C
index a166cdefc9..e048a0a9c2 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/fixCollapsedEdges.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/fixCollapsedEdges.C
@@ -35,7 +35,7 @@ void sammMesh::fixCollapsedEdges()
 {
     cellFaces_.setSize(cellShapes_.size());
 
-    forAll (cellShapes_, cellI)
+    forAll(cellShapes_, cellI)
     {
         cellFaces_[cellI] = cellShapes_[cellI].faces();
     }
@@ -43,23 +43,23 @@ void sammMesh::fixCollapsedEdges()
     // go through the faces and find if there exist faces with duplicate
     // vertices. If so, purge the duplicates and mark the mesh as a polyMesh
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         faceList& curFaces = cellFaces_[cellI];
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             face& vertexLabels = curFaces[faceI];
 
             bool duplicatesFound = false;
 
-            forAll (vertexLabels, vI)
+            forAll(vertexLabels, vI)
             {
                 label curLabel = vertexLabels[vI];
 
                 label nFound = 0;
 
-                forAll (vertexLabels, searchI)
+                forAll(vertexLabels, searchI)
                 {
                     if (vertexLabels[searchI] == curLabel)
                     {
@@ -84,14 +84,14 @@ void sammMesh::fixCollapsedEdges()
                 // here as the shape is still needed to determine which face
                 // of the shape is used in potential couple matches.  This
                 // will be done in the end using the purgeShapes()
-                // 
+                //
 
                 // create a new face without duplicates and replace original
                 face newFace(vertexLabels.size());
 
                 label nNewVertices = 0;
 
-                forAll (vertexLabels, vI)
+                forAll(vertexLabels, vI)
                 {
                     // In order for a face to be a valid entity, duplicate
                     // vertices can only be consecutive (othervise, the
diff --git a/applications/utilities/mesh/conversion/sammToFoam/purgeCellShapes.C b/applications/utilities/mesh/conversion/sammToFoam/purgeCellShapes.C
index 3abcdd645f..0371441341 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/purgeCellShapes.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/purgeCellShapes.C
@@ -33,18 +33,18 @@ Description
 
 void sammMesh::purgeCellShapes()
 {
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         const faceList& curFaces = cellFaces_[cellI];
 
         // Get model faces
         faceList shapeFaces = cellShapes_[cellI].faces();
 
-        forAll (shapeFaces, faceI)
+        forAll(shapeFaces, faceI)
         {
             bool found = false;
 
-            forAll (curFaces, i)
+            forAll(curFaces, i)
             {
                 if (shapeFaces[faceI] == curFaces[i])
                 {
diff --git a/applications/utilities/mesh/conversion/sammToFoam/readCells.C b/applications/utilities/mesh/conversion/sammToFoam/readCells.C
index be6569c096..aa907543ed 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/readCells.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/readCells.C
@@ -163,7 +163,7 @@ void sammMesh::addSAMMcell
     // make a list of labels
     labelList sammCellLabels(curModel.nPoints(), -1);
 
-    forAll (sammCellLabels, labelI)
+    forAll(sammCellLabels, labelI)
     {
         sammCellLabels[labelI] = globalLabels[addressing[labelI]];
     }
@@ -221,7 +221,7 @@ void sammMesh::readCells()
     starCellLabelLookup_.setSize(maxLabel+1);
 
     // reset point labels to invalid value
-    forAll (starCellLabelLookup_, i)
+    forAll(starCellLabelLookup_, i)
     {
         starCellLabelLookup_[i] = -1;
     }
diff --git a/applications/utilities/mesh/conversion/sammToFoam/readCouples.C b/applications/utilities/mesh/conversion/sammToFoam/readCouples.C
index 3d9609e84d..1eb45e48d0 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/readCouples.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/readCouples.C
@@ -128,13 +128,13 @@ void sammMesh::readCouples()
         }
 
         // Once all couples are read, remove zero size faces from all cells
-        forAll (cellFaces_, cellI)
+        forAll(cellFaces_, cellI)
         {
             faceList& curFaces = cellFaces_[cellI];
 
             label zeroSizeFound = 0;
 
-            forAll (curFaces, faceI)
+            forAll(curFaces, faceI)
             {
                 if (curFaces[faceI].empty())
                 {
@@ -151,7 +151,7 @@ void sammMesh::readCouples()
 
                 label nFaces = 0;
 
-                forAll (oldFaces, faceI)
+                forAll(oldFaces, faceI)
                 {
                     if (oldFaces[faceI].size())
                     {
diff --git a/applications/utilities/mesh/conversion/sammToFoam/sammMesh.C b/applications/utilities/mesh/conversion/sammToFoam/sammMesh.C
index 8ed4ddea84..3a57153c07 100644
--- a/applications/utilities/mesh/conversion/sammToFoam/sammMesh.C
+++ b/applications/utilities/mesh/conversion/sammToFoam/sammMesh.C
@@ -116,11 +116,11 @@ void sammMesh::createPolyMeshData()
     label nProblemCells = 0;
 
     // check that there is no zeros in the cellPolys_
-    forAll (cellPolys_, cellI)
+    forAll(cellPolys_, cellI)
     {
         const labelList& curFaceLabels = cellPolys_[cellI];
 
-        forAll (curFaceLabels, faceI)
+        forAll(curFaceLabels, faceI)
         {
             if (curFaceLabels[faceI] == -1)
             {
diff --git a/applications/utilities/mesh/conversion/starToFoam/calcPointCells.C b/applications/utilities/mesh/conversion/starToFoam/calcPointCells.C
index f6e756dc15..685338a4d3 100644
--- a/applications/utilities/mesh/conversion/starToFoam/calcPointCells.C
+++ b/applications/utilities/mesh/conversion/starToFoam/calcPointCells.C
@@ -67,14 +67,14 @@ void starMesh::calcPointCells() const
     // do not come from the cell shape, but from the slaves of the coupled
     // match. It is also adventageous to remove the duplicates from the
     // point-cell addressing, because this removes a lot of waste later.
-    // 
+    //
 
     // For each cell
     forAll(cellShapes_, cellI)
     {
         const faceList& faces = cellFaces_[cellI];
 
-        forAll (faces, i)
+        forAll(faces, i)
         {
             // For each vertex
             const labelList& labels = faces[i];
diff --git a/applications/utilities/mesh/conversion/starToFoam/createBoundaryFaces.C b/applications/utilities/mesh/conversion/starToFoam/createBoundaryFaces.C
index 3875b68345..bc162011c8 100644
--- a/applications/utilities/mesh/conversion/starToFoam/createBoundaryFaces.C
+++ b/applications/utilities/mesh/conversion/starToFoam/createBoundaryFaces.C
@@ -47,11 +47,11 @@ bool starMesh::starEqualFace
 
     label nEqual = 0;
 
-    forAll (cellFace, cellFaceLabelI)
+    forAll(cellFace, cellFaceLabelI)
     {
         const label curCellFaceLabel = cellFace[cellFaceLabelI];
 
-        forAll (boundaryFace, bouFaceLabelI)
+        forAll(boundaryFace, bouFaceLabelI)
         {
             if (boundaryFace[bouFaceLabelI] == curCellFaceLabel)
             {
@@ -70,13 +70,13 @@ bool starMesh::starEqualFace
     // Boundary face is happy if all of its vertices are recognised
     bool boundaryFaceHappy = true;
 
-    forAll (boundaryFace, bouFaceLabelI)
+    forAll(boundaryFace, bouFaceLabelI)
     {
         const label curBouFaceLabel = boundaryFace[bouFaceLabelI];
 
         bool found = false;
 
-        forAll (cellFace, cellFaceLabelI)
+        forAll(cellFace, cellFaceLabelI)
         {
             if (curBouFaceLabel == cellFace[cellFaceLabelI])
             {
@@ -155,7 +155,7 @@ void starMesh::markBoundaryFaces()
                     << " Face : " << endl << curFace << endl
                     << "PROSTAR Command: vset,news,vlis";
 
-                forAll (curFace, spI)
+                forAll(curFace, spI)
                 {
                     if (curFace[spI] > -1 && curFace[spI] < starPointID_.size())
                     {
@@ -185,7 +185,7 @@ void starMesh::collectBoundaryFaces()
         const labelList& curBoundaryCellIDs = boundaryCellIDs_[patchI];
         const labelList& curBoundaryCellFaceIDs = boundaryCellFaceIDs_[patchI];
 
-        forAll (curBoundaryCellIDs, faceI)
+        forAll(curBoundaryCellIDs, faceI)
         {
             patchFaces[faceI] =
                 cellFaces_[curBoundaryCellIDs[faceI]]
diff --git a/applications/utilities/mesh/conversion/starToFoam/createCoupleMatches.C b/applications/utilities/mesh/conversion/starToFoam/createCoupleMatches.C
index d1a074c391..ba13659a18 100644
--- a/applications/utilities/mesh/conversion/starToFoam/createCoupleMatches.C
+++ b/applications/utilities/mesh/conversion/starToFoam/createCoupleMatches.C
@@ -63,7 +63,7 @@ void starMesh::createCoupleMatches()
 
     const label infoJump = max(1000, couples_.size()/20);
 
-    forAll (couples_, coupleI)
+    forAll(couples_, coupleI)
     {
         if (coupleI % infoJump == 0)
         {
@@ -168,7 +168,7 @@ void starMesh::createCoupleMatches()
 
             // Loop through all edges of the master face. For every edge,
             // intersect it with all edges of the cutting face.
-            forAll (masterEdges, masterEdgeI)
+            forAll(masterEdges, masterEdgeI)
             {
                 const edge& curMasterEdge = masterEdges[masterEdgeI];
 
@@ -186,7 +186,7 @@ void starMesh::createCoupleMatches()
                 // go through all slave edges and try to get an intersection.
                 // The point is created along the original master edge rather
                 // than its corrected direction.
-                forAll (slaveEdges, slaveEdgeI)
+                forAll(slaveEdges, slaveEdgeI)
                 {
                     const edge& curSlaveEdge = slaveEdges[slaveEdgeI];
 
@@ -464,7 +464,7 @@ void starMesh::createCoupleMatches()
 
 #           ifdef DEBUG_COUPLE_INTERSECTION
             Info<< "additional slave edge points: " << endl;
-            forAll (slaveEdgePoints, edgeI)
+            forAll(slaveEdgePoints, edgeI)
             {
                 Info<< "edge: " << edgeI << ": " << slaveEdgePoints[edgeI]
                     << endl;
@@ -498,7 +498,7 @@ void starMesh::createCoupleMatches()
             // count the number of additional points for face
             label nAdditionalMasterPoints = 0;
 
-            forAll (masterEdgePoints, edgeI)
+            forAll(masterEdgePoints, edgeI)
             {
                 nAdditionalMasterPoints += masterEdgePoints[edgeI].size();
             }
@@ -516,7 +516,7 @@ void starMesh::createCoupleMatches()
                 << endl;
 #           endif
 
-            forAll (masterEdges, masterEdgeI)
+            forAll(masterEdges, masterEdgeI)
             {
                 // Insert the starting point of the edge
                 tmpMasterFace[nTmpMasterLabels] =
@@ -637,7 +637,7 @@ void starMesh::createCoupleMatches()
             scalar masterTol =
                 cpMergePointTol_*boundBox(tmpMasterFace.points(points_)).mag();
 
-            forAll (mstEdgesToCollapse, edgeI)
+            forAll(mstEdgesToCollapse, edgeI)
             {
 #               ifdef DEBUG_FACE_ORDERING
                 Info<< "edgeI: " << edgeI << " curEdge: "
@@ -693,7 +693,7 @@ void starMesh::createCoupleMatches()
             // count the number of additional points for face
             label nAdditionalSlavePoints = 0;
 
-            forAll (slaveEdgePoints, edgeI)
+            forAll(slaveEdgePoints, edgeI)
             {
                 nAdditionalSlavePoints += slaveEdgePoints[edgeI].size();
             }
@@ -710,7 +710,7 @@ void starMesh::createCoupleMatches()
                 << "nAdditionalSlavePoints: " << nAdditionalSlavePoints << endl;
 #           endif
 
-            forAll (slaveEdges, slaveEdgeI)
+            forAll(slaveEdges, slaveEdgeI)
             {
                 // Insert the starting point of the edge
                 tmpSlaveFace[nTmpSlaveLabels] =
@@ -894,7 +894,7 @@ void starMesh::createCoupleMatches()
 
             vector masterProjDir = -newMasterFace.normal(points_);
 
-            forAll (newSlaveEdges, edgeI)
+            forAll(newSlaveEdges, edgeI)
             {
                 // Take the slave edge points and project into the master.
                 // In order to create a good intersection, move the
@@ -937,7 +937,7 @@ void starMesh::createCoupleMatches()
             {
                 vector slaveProjDir = -newSlaveFace.normal(points_);
 
-                forAll (newMasterEdges, edgeI)
+                forAll(newMasterEdges, edgeI)
                 {
                     // Take the edge master points and project into the slave.
                     // In order to create a good intersection, move the
@@ -1053,7 +1053,7 @@ void starMesh::createCoupleMatches()
                     SLList<edge> edgesToConsider;
 
                     // collect master edges
-                    forAll (newMasterEdges, edgeI)
+                    forAll(newMasterEdges, edgeI)
                     {
                         const edge& cme = newMasterEdges[edgeI];
 
@@ -1072,7 +1072,7 @@ void starMesh::createCoupleMatches()
                     }
 
                     // collect slave edges
-                    forAll (newSlaveEdges, edgeI)
+                    forAll(newSlaveEdges, edgeI)
                     {
                         const edge& cse = newSlaveEdges[edgeI];
 
@@ -1260,7 +1260,7 @@ void starMesh::createCoupleMatches()
 #               endif
 
                 // check the intersection face for duplicate points
-                forAll (intersectedFace, checkI)
+                forAll(intersectedFace, checkI)
                 {
                     for
                     (
@@ -1312,7 +1312,7 @@ void starMesh::createCoupleMatches()
             // onto the master face to ensure closedness
             vector pointProjectionNormal = -masterFace.normal(points_);
 
-            forAll (intersectedFace, intPointI)
+            forAll(intersectedFace, intPointI)
             {
 #               ifdef DEBUG_COUPLE_PROJECTION
                 Info<< "Proj: old point: "
@@ -1424,7 +1424,7 @@ void starMesh::createCoupleMatches()
         // Loop through all cells and reset faces for removal to zero size
         const labelList crfToc = cellRemovedFaces.toc();
 
-        forAll (crfToc, cellI)
+        forAll(crfToc, cellI)
         {
             const label curCell = crfToc[cellI];
 
@@ -1451,7 +1451,7 @@ void starMesh::createCoupleMatches()
         const labelList cafToc = cellAddedFaces.toc();
 
         // Insert the new faces into the list
-        forAll (cafToc, cellI)
+        forAll(cafToc, cellI)
         {
             const label curCell = cafToc[cellI];
 
@@ -1465,7 +1465,7 @@ void starMesh::createCoupleMatches()
             label nNewFaces = 0;
 
             // copy original faces that have not been removed
-            forAll (oldFaces, faceI)
+            forAll(oldFaces, faceI)
             {
                 if (oldFaces[faceI].size())
                 {
diff --git a/applications/utilities/mesh/conversion/starToFoam/createPolyBoundary.C b/applications/utilities/mesh/conversion/starToFoam/createPolyBoundary.C
index e8eb4ec479..47cb1208a8 100644
--- a/applications/utilities/mesh/conversion/starToFoam/createPolyBoundary.C
+++ b/applications/utilities/mesh/conversion/starToFoam/createPolyBoundary.C
@@ -42,13 +42,13 @@ void starMesh::createPolyBoundary()
 
     const labelListList& PointCells = pointCells();
 
-    forAll (boundary_, patchI)
+    forAll(boundary_, patchI)
     {
         const faceList& curShapePatch = boundary_[patchI];
 
         polyBoundaryPatchStartIndices_[patchI] = nCreatedFaces;
 
-        forAll (curShapePatch, faceI)
+        forAll(curShapePatch, faceI)
         {
             bool found = false;
 
@@ -100,7 +100,7 @@ void starMesh::createPolyBoundary()
                                         << endl;
 
                                     Info<< "PROSTAR Command: vset,news,vlis";
-                                    forAll (curCellFaces[cellFaceI], spI)
+                                    forAll(curCellFaces[cellFaceI], spI)
                                     {
                                         // check if the point is given by STAR
                                         // or created locally
@@ -136,7 +136,7 @@ void starMesh::createPolyBoundary()
                                         << endl;
 
                                     Info<< "PROSTAR Command: vset,news,vlis";
-                                    forAll (curCellFaces[cellFaceI], spI)
+                                    forAll(curCellFaces[cellFaceI], spI)
                                     {
                                         // check if the point is given by STAR
                                         // or created locally
@@ -181,11 +181,11 @@ void starMesh::createPolyBoundary()
     // check all cellPolys_ to see if there are any missing faces
     label nMissingFaceFound = 0;
 
-    forAll (cellPolys_, cellI)
+    forAll(cellPolys_, cellI)
     {
         const labelList& curFaces = cellPolys_[cellI];
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             if (curFaces[faceI] < 0)
             {
@@ -200,7 +200,7 @@ void starMesh::createPolyBoundary()
                 nMissingFaceFound++;
 
                 Info<< "PROSTAR Command: vset,news,vlis";
-                forAll (missingFace, spI)
+                forAll(missingFace, spI)
                 {
                     // check if the point is given by STAR or created locally
                     if
@@ -233,11 +233,11 @@ void starMesh::createPolyBoundary()
     // (faces addressed once or more than twice)
     labelList markupFaces(meshFaces_.size(), 0);
 
-    forAll (cellPolys_, cellI)
+    forAll(cellPolys_, cellI)
     {
         const labelList& curFaces = cellPolys_[cellI];
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             markupFaces[curFaces[faceI]]++;
         }
@@ -250,7 +250,7 @@ void starMesh::createPolyBoundary()
 
     label nProblemFacesFound = 0;
 
-    forAll (markupFaces, faceI)
+    forAll(markupFaces, faceI)
     {
         if (markupFaces[faceI] != 2)
         {
@@ -264,7 +264,7 @@ void starMesh::createPolyBoundary()
             nProblemFacesFound++;
 
             Info<< "PROSTAR Command: vset,news,vlis";
-            forAll (problemFace, spI)
+            forAll(problemFace, spI)
             {
                 // check if the point is given by STAR or created locally
                 if
@@ -299,7 +299,7 @@ List<polyPatch*> starMesh::polyBoundaryPatches(const polyMesh& pMesh)
 {
     List<polyPatch*> p(boundary_.size());
 
-    forAll (boundary_, patchI)
+    forAll(boundary_, patchI)
     {
         p[patchI] = polyPatch::New
         (
diff --git a/applications/utilities/mesh/conversion/starToFoam/createPolyCells.C b/applications/utilities/mesh/conversion/starToFoam/createPolyCells.C
index c4847fc390..d63a5b22d0 100644
--- a/applications/utilities/mesh/conversion/starToFoam/createPolyCells.C
+++ b/applications/utilities/mesh/conversion/starToFoam/createPolyCells.C
@@ -41,13 +41,13 @@ void starMesh::createPolyCells()
 
     label maxFaces = 0;
 
-    forAll (cellPolys_, cellI)
+    forAll(cellPolys_, cellI)
     {
         cell& curCell = cellPolys_[cellI];
 
         curCell.setSize(cellFaces_[cellI].size());
 
-        forAll (curCell, fI)
+        forAll(curCell, fI)
         {
             curCell[fI] = -1;
         }
@@ -145,7 +145,7 @@ void starMesh::createPolyCells()
             label nextNei = -1;
             label minNei = cellPolys_.size();
 
-            forAll (neiCells, ncI)
+            forAll(neiCells, ncI)
             {
                 if (neiCells[ncI] > -1 && neiCells[ncI] < minNei)
                 {
diff --git a/applications/utilities/mesh/conversion/starToFoam/fixCollapsedEdges.C b/applications/utilities/mesh/conversion/starToFoam/fixCollapsedEdges.C
index 545c98408c..d985f76495 100644
--- a/applications/utilities/mesh/conversion/starToFoam/fixCollapsedEdges.C
+++ b/applications/utilities/mesh/conversion/starToFoam/fixCollapsedEdges.C
@@ -35,7 +35,7 @@ void starMesh::fixCollapsedEdges()
 {
     cellFaces_.setSize(cellShapes_.size());
 
-    forAll (cellShapes_, cellI)
+    forAll(cellShapes_, cellI)
     {
         cellFaces_[cellI] = cellShapes_[cellI].faces();
     }
@@ -43,23 +43,23 @@ void starMesh::fixCollapsedEdges()
     // go through the faces and find if there exist faces with duplicate
     // vertices. If so, purge the duplicates and mark the mesh as a polyMesh
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         faceList& curFaces = cellFaces_[cellI];
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             face& vertexLabels = curFaces[faceI];
 
             bool duplicatesFound = false;
 
-            forAll (vertexLabels, vI)
+            forAll(vertexLabels, vI)
             {
                 label curLabel = vertexLabels[vI];
 
                 label nFound = 0;
 
-                forAll (vertexLabels, searchI)
+                forAll(vertexLabels, searchI)
                 {
                     if (vertexLabels[searchI] == curLabel)
                     {
@@ -84,13 +84,13 @@ void starMesh::fixCollapsedEdges()
                 // here as the shape is still needed to determine which face
                 // of the shape is used in potential couple matches.  This
                 // will be done in the end using the purgeShapes()
-                // 
+                //
 
                 // create a new face without duplicates and replace original
                 face newFace(vertexLabels.size());
                 label nNewVertices = 0;
 
-                forAll (vertexLabels, vI)
+                forAll(vertexLabels, vI)
                 {
                     // In order for a face to be a valid entity, duplicate
                     // vertices can only be consecutive (othervise, the
diff --git a/applications/utilities/mesh/conversion/starToFoam/mergeCoupleFacePoints.C b/applications/utilities/mesh/conversion/starToFoam/mergeCoupleFacePoints.C
index 73ca9628d8..78c350ab59 100644
--- a/applications/utilities/mesh/conversion/starToFoam/mergeCoupleFacePoints.C
+++ b/applications/utilities/mesh/conversion/starToFoam/mergeCoupleFacePoints.C
@@ -63,7 +63,7 @@ void starMesh::mergeCoupleFacePoints()
 
     label nMergeSets = 0;
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         const faceList& curFaces = cellFaces_[cellI];
 
@@ -74,13 +74,13 @@ void starMesh::mergeCoupleFacePoints()
 
         scalar pointMergeTol = GREAT;
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             nPointsInCell += curFaces[faceI].size();
 
             edgeList curEdges = curFaces[faceI].edges();
 
-            forAll (curEdges, edgeI)
+            forAll(curEdges, edgeI)
             {
                 scalar length = curEdges[edgeI].mag(points_);
 
@@ -98,11 +98,11 @@ void starMesh::mergeCoupleFacePoints()
         labelList cellPoints(nPointsInCell);
         label nAddedPoints = 0;
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             const face& f = curFaces[faceI];
 
-            forAll (f, fI)
+            forAll(f, fI)
             {
                 cellPoints[nAddedPoints] = f[fI];
                 nAddedPoints++;
@@ -210,7 +210,7 @@ void starMesh::mergeCoupleFacePoints()
                                 << minMerge << " and " << maxMerge << endl;
 #                           endif
 
-                            forAll (renumberPoints, elimI)
+                            forAll(renumberPoints, elimI)
                             {
                                 if (renumberPoints[elimI] == maxMerge)
                                 {
@@ -235,7 +235,7 @@ void starMesh::mergeCoupleFacePoints()
 
     // Insert the primary point renumbering into the list
     // Take care of possibly unused points in the list
-    forAll (renumberPoints, pointI)
+    forAll(renumberPoints, pointI)
     {
         if (renumberPoints[pointI] < 0)
         {
@@ -253,17 +253,17 @@ void starMesh::mergeCoupleFacePoints()
     // renumbering of all faces.  This will only be used to see which
     // points are still used!
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         faceList& prelimFaces = cellFaces_[cellI];
 
-        forAll (prelimFaces, faceI)
+        forAll(prelimFaces, faceI)
         {
             face oldFacePoints = prelimFaces[faceI];
 
             face& prelimFacePoints = prelimFaces[faceI];
 
-            forAll (prelimFacePoints, pointI)
+            forAll(prelimFacePoints, pointI)
             {
                 if (renumberPoints[oldFacePoints[pointI]] < 0)
                 {
@@ -284,26 +284,26 @@ void starMesh::mergeCoupleFacePoints()
     // re-create the point list and renumber the whole lot
     renumberPoints = 0;
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         const faceList& curFaces = cellFaces_[cellI];
 
-        forAll (curFaces, faceI)
+        forAll(curFaces, faceI)
         {
             const face& curFacePoints = curFaces[faceI];
 
-            forAll (curFacePoints, pointI)
+            forAll(curFacePoints, pointI)
             {
                 renumberPoints[curFacePoints[pointI]]++;
             }
         }
     }
 
-    forAll (cellShapes_, cellI)
+    forAll(cellShapes_, cellI)
     {
         const labelList& curLabels = cellShapes_[cellI];
 
-        forAll (curLabels, pointI)
+        forAll(curLabels, pointI)
         {
             if (renumberPoints[curLabels[pointI]] == 0)
             {
@@ -320,7 +320,7 @@ void starMesh::mergeCoupleFacePoints()
 
     label nUsedPoints = 0;
 
-    forAll (renumberPoints, pointI)
+    forAll(renumberPoints, pointI)
     {
         if (renumberPoints[pointI] > 0)
         {
@@ -347,17 +347,17 @@ void starMesh::mergeCoupleFacePoints()
 
     Info<< "Renumbering all faces" << endl;
 
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         faceList& newFaces = cellFaces_[cellI];
 
-        forAll (newFaces, faceI)
+        forAll(newFaces, faceI)
         {
             face oldFacePoints = newFaces[faceI];
 
             face& newFacePoints = newFaces[faceI];
 
-            forAll (newFacePoints, pointI)
+            forAll(newFacePoints, pointI)
             {
                 if (renumberPoints[oldFacePoints[pointI]] < 0)
                 {
@@ -376,13 +376,13 @@ void starMesh::mergeCoupleFacePoints()
 
     Info<< "Renumbering all cell shapes" << endl;
 
-    forAll (cellShapes_, cellI)
+    forAll(cellShapes_, cellI)
     {
         labelList oldLabels = cellShapes_[cellI];
 
         labelList& curLabels = cellShapes_[cellI];
 
-        forAll (curLabels, pointI)
+        forAll(curLabels, pointI)
         {
             if (renumberPoints[curLabels[pointI]] < 0)
             {
@@ -408,7 +408,7 @@ void starMesh::mergeCoupleFacePoints()
 
     starPointID_ = -1;
 
-    forAll (starPointID_, pointI)
+    forAll(starPointID_, pointI)
     {
         if (renumberPoints[pointI] > -1)
         {
diff --git a/applications/utilities/mesh/conversion/starToFoam/purgeCellShapes.C b/applications/utilities/mesh/conversion/starToFoam/purgeCellShapes.C
index 92dedc00e6..ba208afb0e 100644
--- a/applications/utilities/mesh/conversion/starToFoam/purgeCellShapes.C
+++ b/applications/utilities/mesh/conversion/starToFoam/purgeCellShapes.C
@@ -33,18 +33,18 @@ Description
 
 void starMesh::purgeCellShapes()
 {
-    forAll (cellFaces_, cellI)
+    forAll(cellFaces_, cellI)
     {
         const faceList& curFaces = cellFaces_[cellI];
 
         // Get model faces
         faceList shapeFaces = cellShapes_[cellI].faces();
 
-        forAll (shapeFaces, faceI)
+        forAll(shapeFaces, faceI)
         {
             bool found = false;
 
-            forAll (curFaces, i)
+            forAll(curFaces, i)
             {
                 if (shapeFaces[faceI] == curFaces[i])
                 {
diff --git a/applications/utilities/mesh/conversion/starToFoam/readCells.C b/applications/utilities/mesh/conversion/starToFoam/readCells.C
index 9e2d5efd6a..821145f06f 100644
--- a/applications/utilities/mesh/conversion/starToFoam/readCells.C
+++ b/applications/utilities/mesh/conversion/starToFoam/readCells.C
@@ -108,7 +108,7 @@ void starMesh::addRegularCell
     // get reference to the addressing list
     const label* addressing = regularAddressingTable[regularTypeFlag];
 
-    forAll (regularCellLabels, labelI)
+    forAll(regularCellLabels, labelI)
     {
         regularCellLabels[labelI] = labels[addressing[labelI]];
     }
@@ -125,7 +125,7 @@ void starMesh::addSAMMcell
 {
     // get type, reg and permutation flag
     label typeFlag = labels[21];
-//     label regularityFlag = labels[22];  // Not used. 
+//     label regularityFlag = labels[22];  // Not used.
     label permutationFlag = labels[23];
 
     // grab the shape from the table
@@ -194,7 +194,7 @@ void starMesh::addSAMMcell
     // get reference to the addressing list
     const label* addressing = sammAddressingTable[sammTypeFlag];
 
-    forAll (sammCellLabels, labelI)
+    forAll(sammCellLabels, labelI)
     {
         sammCellLabels[labelI] = labels[addressing[labelI]];
     }
@@ -239,7 +239,7 @@ void starMesh::readCells()
                 }
 
                 // backward compatibility: number of trailing rubbish in
-                // STAR is unknown. 
+                // STAR is unknown.
                 // Fixed to cope with missing \n on last line.
                 readToNl(cellsFile);
             }
@@ -259,15 +259,15 @@ void starMesh::readCells()
     starCellPermutation_.setSize(nCells);
 
     // reset permutation to invalid value
-    forAll (starCellPermutation_, i)
+    forAll(starCellPermutation_, i)
     {
         starCellPermutation_[i] = -1;
     }
 
     starCellLabelLookup_.setSize(maxLabel+1);
-    
+
     // reset point labels to invalid value
-    forAll (starCellLabelLookup_, i)
+    forAll(starCellLabelLookup_, i)
     {
         starCellLabelLookup_[i] = -1;
     }
@@ -285,8 +285,8 @@ void starMesh::readCells()
 
             label addOnToCell = 0;
 
-            // reset the labels to -1. Debugging. 
-            forAll (labels, i)
+            // reset the labels to -1. Debugging.
+            forAll(labels, i)
             {
                 labels[i] = -1;
             }
@@ -342,7 +342,7 @@ void starMesh::readCells()
                 }
 
                 // backward compatibility: number of trailing rubbish in
-                // STAR is unknown. 
+                // STAR is unknown.
                 readToNl(cellsFile);
 
                 addOnToCell--;
@@ -367,13 +367,13 @@ void starMesh::readCells()
             // check cell labels
             const labelList& curShapeLabels = cellShapes_[cellI];
 
-            forAll (curShapeLabels, i)
+            forAll(curShapeLabels, i)
             {
                 if (curShapeLabels[i] < 0)
                 {
                     FatalErrorIn("starMesh::readCells()")
                         << "Invalid vertex found in cell " << cellI
-                        << ". STAR cell no: " << lineLabel 
+                        << ". STAR cell no: " << lineLabel
                         << " labels: " << curShapeLabels
                         << abort(FatalError);
                 }
diff --git a/applications/utilities/mesh/generation/blockMesh/mergePatchPairs.H b/applications/utilities/mesh/generation/blockMesh/mergePatchPairs.H
index 2fd2fe2dc1..c8a25abeec 100644
--- a/applications/utilities/mesh/generation/blockMesh/mergePatchPairs.H
+++ b/applications/utilities/mesh/generation/blockMesh/mergePatchPairs.H
@@ -7,7 +7,7 @@
             List<faceZone*> fz(3*mergePatchPairs.size());
             List<cellZone*> cz(0);
 
-            forAll (mergePatchPairs, pairI)
+            forAll(mergePatchPairs, pairI)
             {
                 const word mergeName
                 (
@@ -34,7 +34,7 @@
 
                 labelList isf(masterPatch.size());
 
-                forAll (isf, i)
+                forAll(isf, i)
                 {
                     isf[i] = masterPatch.start() + i;
                 }
@@ -58,7 +58,7 @@
 
                 labelList osf(slavePatch.size());
 
-                forAll (osf, i)
+                forAll(osf, i)
                 {
                     osf[i] = slavePatch.start() + i;
                 }
@@ -91,7 +91,7 @@
             attachPolyTopoChanger polyMeshAttacher(mesh);
             polyMeshAttacher.setSize(mergePatchPairs.size());
 
-            forAll (mergePatchPairs, pairI)
+            forAll(mergePatchPairs, pairI)
             {
                 const word mergeName
                 (
diff --git a/applications/utilities/mesh/manipulation/mergeMeshes/mergePolyMesh.C b/applications/utilities/mesh/manipulation/mergeMeshes/mergePolyMesh.C
index d1b633a4b9..08df4619fa 100644
--- a/applications/utilities/mesh/manipulation/mergeMeshes/mergePolyMesh.C
+++ b/applications/utilities/mesh/manipulation/mergeMeshes/mergePolyMesh.C
@@ -48,7 +48,7 @@ Foam::label Foam::mergePolyMesh::patchIndex(const polyPatch& p)
 
     bool nameFound = false;
 
-    forAll (patchNames_, patchI)
+    forAll(patchNames_, patchI)
     {
         if (patchNames_[patchI] == pName)
         {
@@ -98,7 +98,7 @@ Foam::label Foam::mergePolyMesh::zoneIndex
     const word& curName
 )
 {
-    forAll (names, zoneI)
+    forAll(names, zoneI)
     {
         if (names[zoneI] == curName)
         {
@@ -129,7 +129,7 @@ Foam::mergePolyMesh::mergePolyMesh(const IOobject& io)
     wordList curPatchTypes = boundaryMesh().types();
     wordList curPatchNames = boundaryMesh().names();
 
-    forAll (curPatchTypes, patchI)
+    forAll(curPatchTypes, patchI)
     {
         patchTypes_.append(curPatchTypes[patchI]);
         patchNames_.append(curPatchNames[patchI]);
@@ -144,7 +144,7 @@ Foam::mergePolyMesh::mergePolyMesh(const IOobject& io)
         pointZoneNames_.setCapacity(2*curPointZoneNames.size());
     }
 
-    forAll (curPointZoneNames, zoneI)
+    forAll(curPointZoneNames, zoneI)
     {
         pointZoneNames_.append(curPointZoneNames[zoneI]);
     }
@@ -156,7 +156,7 @@ Foam::mergePolyMesh::mergePolyMesh(const IOobject& io)
     {
         faceZoneNames_.setCapacity(2*curFaceZoneNames.size());
     }
-    forAll (curFaceZoneNames, zoneI)
+    forAll(curFaceZoneNames, zoneI)
     {
         faceZoneNames_.append(curFaceZoneNames[zoneI]);
     }
@@ -168,7 +168,7 @@ Foam::mergePolyMesh::mergePolyMesh(const IOobject& io)
     {
         cellZoneNames_.setCapacity(2*curCellZoneNames.size());
     }
-    forAll (curCellZoneNames, zoneI)
+    forAll(curCellZoneNames, zoneI)
     {
         cellZoneNames_.append(curCellZoneNames[zoneI]);
     }
@@ -194,12 +194,12 @@ void Foam::mergePolyMesh::addMesh(const polyMesh& m)
     const pointZoneMesh& pz = m.pointZones();
     labelList pointZoneIndices(pz.size());
 
-    forAll (pz, zoneI)
+    forAll(pz, zoneI)
     {
         pointZoneIndices[zoneI] = zoneIndex(pointZoneNames_, pz[zoneI].name());
     }
 
-    forAll (p, pointI)
+    forAll(p, pointI)
     {
         // Grab zone ID.  If a point is not in a zone, it will return -1
         zoneID = pz.whichZone(pointI);
@@ -231,12 +231,12 @@ void Foam::mergePolyMesh::addMesh(const polyMesh& m)
     const cellZoneMesh& cz = m.cellZones();
     labelList cellZoneIndices(cz.size());
 
-    forAll (cz, zoneI)
+    forAll(cz, zoneI)
     {
         cellZoneIndices[zoneI] = zoneIndex(cellZoneNames_, cz[zoneI].name());
     }
 
-    forAll (c, cellI)
+    forAll(c, cellI)
     {
         // Grab zone ID.  If a cell is not in a zone, it will return -1
         zoneID = cz.whichZone(cellI);
@@ -267,7 +267,7 @@ void Foam::mergePolyMesh::addMesh(const polyMesh& m)
     // Gather the patch indices
     labelList patchIndices(bm.size());
 
-    forAll (patchIndices, patchI)
+    forAll(patchIndices, patchI)
     {
         patchIndices[patchI] = patchIndex(bm[patchI]);
     }
@@ -281,7 +281,7 @@ void Foam::mergePolyMesh::addMesh(const polyMesh& m)
     const faceZoneMesh& fz = m.faceZones();
     labelList faceZoneIndices(fz.size());
 
-    forAll (fz, zoneI)
+    forAll(fz, zoneI)
     {
         faceZoneIndices[zoneI] = zoneIndex(faceZoneNames_, fz[zoneI].name());
     }
@@ -295,13 +295,13 @@ void Foam::mergePolyMesh::addMesh(const polyMesh& m)
     label newOwn, newNei, newPatch, newZone;
     bool newZoneFlip;
 
-    forAll (f, faceI)
+    forAll(f, faceI)
     {
         const face& curFace = f[faceI];
 
         face newFace(curFace.size());
 
-        forAll (curFace, pointI)
+        forAll(curFace, pointI)
         {
             newFace[pointI] = renumberPoints[curFace[pointI]];
         }
diff --git a/applications/utilities/mesh/manipulation/mirrorMesh/mirrorFvMesh.C b/applications/utilities/mesh/manipulation/mirrorMesh/mirrorFvMesh.C
index 0c97671590..52735eb41a 100644
--- a/applications/utilities/mesh/manipulation/mirrorMesh/mirrorFvMesh.C
+++ b/applications/utilities/mesh/manipulation/mirrorMesh/mirrorFvMesh.C
@@ -82,13 +82,13 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     labelList mirrorPointLookup(oldPoints.size(), -1);
 
     // Grab the old points
-    forAll (oldPoints, pointI)
+    forAll(oldPoints, pointI)
     {
         newPoints[nNewPoints] = oldPoints[pointI];
         nNewPoints++;
     }
 
-    forAll (oldPoints, pointI)
+    forAll(oldPoints, pointI)
     {
         scalar alpha =
             mirrorPlane.normalIntersect
@@ -145,7 +145,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
 
     const unallocLabelList& oldOwnerStart = lduAddr().ownerStartAddr();
 
-    forAll (newCellFaces, cellI)
+    forAll(newCellFaces, cellI)
     {
         labelList& curFaces = newCellFaces[cellI];
 
@@ -154,7 +154,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
 
         curFaces.setSize(e - s);
 
-        forAll (curFaces, i)
+        forAll(curFaces, i)
         {
             curFaces[i] = s + i;
         }
@@ -164,7 +164,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     // as internal
     boolListList insertedBouFace(oldPatches.size());
 
-    forAll (oldPatches, patchI)
+    forAll(oldPatches, patchI)
     {
         const polyPatch& curPatch = oldPatches[patchI];
 
@@ -186,7 +186,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
         const unallocLabelList& curFaceCells = curPatch.faceCells();
         const label curStart = curPatch.start();
 
-        forAll (curPatch, faceI)
+        forAll(curPatch, faceI)
         {
             // Find out if the mirrored face is identical to the
             // original.  If so, the face needs to become internal and
@@ -194,7 +194,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
             const face& origFace = curPatch[faceI];
 
             face mirrorFace(origFace.size());
-            forAll (mirrorFace, pointI)
+            forAll(mirrorFace, pointI)
             {
                 mirrorFace[pointI] = mirrorPointLookup[origFace[pointI]];
             }
@@ -225,11 +225,11 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     label nNewFaces = 0;
 
     // Insert original (internal) faces
-    forAll (newCellFaces, cellI)
+    forAll(newCellFaces, cellI)
     {
         const labelList& curCellFaces = newCellFaces[cellI];
 
-        forAll (curCellFaces, cfI)
+        forAll(curCellFaces, cfI)
         {
             newFaces[nNewFaces] = oldFaces[curCellFaces[cfI]];
             masterFaceLookup[curCellFaces[cfI]] = nNewFaces;
@@ -264,7 +264,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     labelList newPatchStarts(boundary().size(), -1);
     label nNewPatches = 0;
 
-    forAll (boundaryMesh(), patchI)
+    forAll(boundaryMesh(), patchI)
     {
         const label curPatchSize = boundaryMesh()[patchI].size();
         const label curPatchStart = boundaryMesh()[patchI].start();
@@ -346,14 +346,14 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     label nNewCells = 0;
 
     // Grab the original cells.  Take care of face renumbering.
-    forAll (oldCells, cellI)
+    forAll(oldCells, cellI)
     {
         const cell& oc = oldCells[cellI];
 
         cell& nc = newCells[nNewCells];
         nc.setSize(oc.size());
 
-        forAll (oc, i)
+        forAll(oc, i)
         {
             nc[i] = masterFaceLookup[oc[i]];
         }
@@ -362,14 +362,14 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     }
 
     // Mirror the cells
-    forAll (oldCells, cellI)
+    forAll(oldCells, cellI)
     {
         const cell& oc = oldCells[cellI];
 
         cell& nc = newCells[nNewCells];
         nc.setSize(oc.size());
 
-        forAll (oc, i)
+        forAll(oc, i)
         {
             nc[i] = mirrorFaceLookup[oc[i]];
         }
@@ -394,7 +394,7 @@ Foam::mirrorFvMesh::mirrorFvMesh(const IOobject& io)
     // Add the boundary patches
     List<polyPatch*> p(newPatchTypes.size());
 
-    forAll (p, patchI)
+    forAll(p, patchI)
     {
         p[patchI] = polyPatch::New
         (
diff --git a/applications/utilities/mesh/manipulation/renumberMesh/renumberMesh.C b/applications/utilities/mesh/manipulation/renumberMesh/renumberMesh.C
index 56d63f74c9..099cdcdead 100644
--- a/applications/utilities/mesh/manipulation/renumberMesh/renumberMesh.C
+++ b/applications/utilities/mesh/manipulation/renumberMesh/renumberMesh.C
@@ -125,7 +125,7 @@ labelList regionFaceOrder
 
     label prevRegion = -1;
 
-    forAll (cellOrder, newCellI)
+    forAll(cellOrder, newCellI)
     {
         label oldCellI = cellOrder[newCellI];
 
diff --git a/applications/utilities/mesh/manipulation/stitchMesh/stitchMesh.C b/applications/utilities/mesh/manipulation/stitchMesh/stitchMesh.C
index a6b9e2cb3e..70927ca8f3 100644
--- a/applications/utilities/mesh/manipulation/stitchMesh/stitchMesh.C
+++ b/applications/utilities/mesh/manipulation/stitchMesh/stitchMesh.C
@@ -284,7 +284,7 @@ int main(int argc, char *argv[])
 
         labelList osf(slavePatch.size());
 
-        forAll (osf, i)
+        forAll(osf, i)
         {
             osf[i] = slavePatch.start() + i;
         }
diff --git a/applications/utilities/miscellaneous/foamFormatConvert/foamFormatConvert.C b/applications/utilities/miscellaneous/foamFormatConvert/foamFormatConvert.C
index 0140ba4e23..7e4363ffb0 100644
--- a/applications/utilities/miscellaneous/foamFormatConvert/foamFormatConvert.C
+++ b/applications/utilities/miscellaneous/foamFormatConvert/foamFormatConvert.C
@@ -65,7 +65,7 @@ int main(int argc, char *argv[])
 #   include "createTime.H"
     Foam::instantList timeDirs = Foam::timeSelector::select0(runTime, args);
 
-    forAll (timeDirs, timeI)
+    forAll(timeDirs, timeI)
     {
         runTime.setTime(timeDirs[timeI], timeI);
         Info<< "Time = " << runTime.timeName() << endl;
diff --git a/applications/utilities/miscellaneous/foamInfoExec/foamInfoExec.C b/applications/utilities/miscellaneous/foamInfoExec/foamInfoExec.C
index 6f78a733b8..41abe8dcdb 100644
--- a/applications/utilities/miscellaneous/foamInfoExec/foamInfoExec.C
+++ b/applications/utilities/miscellaneous/foamInfoExec/foamInfoExec.C
@@ -56,7 +56,7 @@ int main(int argc, char *argv[])
             Foam::Time::findTimes(args.rootPath()/args.caseName())
         );
 
-        forAll (times, i)
+        forAll(times, i)
         {
             Info<< times[i].name() << endl;
         }
diff --git a/applications/utilities/parallelProcessing/decomposePar/domainDecomposition.C b/applications/utilities/parallelProcessing/decomposePar/domainDecomposition.C
index 392594cc62..a286a2facf 100644
--- a/applications/utilities/parallelProcessing/decomposePar/domainDecomposition.C
+++ b/applications/utilities/parallelProcessing/decomposePar/domainDecomposition.C
@@ -118,7 +118,7 @@ bool Foam::domainDecomposition::writeDecomposition()
     // Make a lookup map for globally shared points
     Map<label> sharedPointLookup(2*globallySharedPoints_.size());
 
-    forAll (globallySharedPoints_, pointi)
+    forAll(globallySharedPoints_, pointi)
     {
         sharedPointLookup.insert(globallySharedPoints_[pointi], pointi);
     }
@@ -174,7 +174,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
         pointField procPoints(curPointLabels.size());
 
-        forAll (curPointLabels, pointi)
+        forAll(curPointLabels, pointi)
         {
             procPoints[pointi] = meshPoints[curPointLabels[pointi]];
 
@@ -190,7 +190,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
         faceList procFaces(curFaceLabels.size());
 
-        forAll (curFaceLabels, facei)
+        forAll(curFaceLabels, facei)
         {
             // Mark the original face as used
             // Remember to decrement the index by one (turning index)
@@ -217,7 +217,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
             procFaceLabels.setSize(origFaceLabels.size());
 
-            forAll (origFaceLabels, pointi)
+            forAll(origFaceLabels, pointi)
             {
                 procFaceLabels[pointi] = pointLookup[origFaceLabels[pointi]];
             }
@@ -230,7 +230,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
         cellList procCells(curCellLabels.size());
 
-        forAll (curCellLabels, celli)
+        forAll(curCellLabels, celli)
         {
             const labelList& origCellLabels = meshCells[curCellLabels[celli]];
 
@@ -238,7 +238,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
             curCell.setSize(origCellLabels.size());
 
-            forAll (origCellLabels, cellFaceI)
+            forAll(origCellLabels, cellFaceI)
             {
                 curCell[cellFaceI] = faceLookup[origCellLabels[cellFaceI]];
             }
@@ -305,7 +305,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
         label nPatches = 0;
 
-        forAll (curPatchSizes, patchi)
+        forAll(curPatchSizes, patchi)
         {
             // Get the face labels consistent with the field mapping
             // (reuse the patch field mappers)
@@ -336,7 +336,7 @@ bool Foam::domainDecomposition::writeDecomposition()
             nPatches++;
         }
 
-        forAll (curProcessorPatchSizes, procPatchI)
+        forAll(curProcessorPatchSizes, procPatchI)
         {
             procPatches[nPatches] =
                 new processorPolyPatch
@@ -377,7 +377,7 @@ bool Foam::domainDecomposition::writeDecomposition()
 
             // Use the pointToZone map to find out the single zone (if any),
             // use slow search only for shared points.
-            forAll (curPointLabels, pointi)
+            forAll(curPointLabels, pointi)
             {
                 label curPoint = curPointLabels[pointi];
 
@@ -448,7 +448,7 @@ bool Foam::domainDecomposition::writeDecomposition()
             // Go through all the zoned faces and find out if they
             // belong to a zone.  If so, add it to the zone as
             // necessary
-            forAll (curFaceLabels, facei)
+            forAll(curFaceLabels, facei)
             {
                 // Remember to decrement the index by one (turning index)
                 //
@@ -535,7 +535,7 @@ bool Foam::domainDecomposition::writeDecomposition()
                 zoneCells[zoneI].setCapacity(cz[zoneI].size() / nProcs_);
             }
 
-            forAll (curCellLabels, celli)
+            forAll(curCellLabels, celli)
             {
                 label curCellI = curCellLabels[celli];
 
@@ -598,7 +598,7 @@ bool Foam::domainDecomposition::writeDecomposition()
         label nProcPatches = 0;
         label nProcFaces = 0;
 
-        forAll (procMesh.boundaryMesh(), patchi)
+        forAll(procMesh.boundaryMesh(), patchi)
         {
             if
             (
diff --git a/applications/utilities/parallelProcessing/decomposePar/domainDecompositionMesh.C b/applications/utilities/parallelProcessing/decomposePar/domainDecompositionMesh.C
index 1b31d9edd5..644655f0a7 100644
--- a/applications/utilities/parallelProcessing/decomposePar/domainDecompositionMesh.C
+++ b/applications/utilities/parallelProcessing/decomposePar/domainDecompositionMesh.C
@@ -65,7 +65,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
     {
         List<SLList<label> > procCellList(nProcs_);
 
-        forAll (cellToProc_, celli)
+        forAll(cellToProc_, celli)
         {
             if (cellToProc_[celli] >= nProcs_)
             {
@@ -81,7 +81,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
         }
 
         // Convert linked lists into normal lists
-        forAll (procCellList, procI)
+        forAll(procCellList, procI)
         {
             procCellAddressing_[procI] = procCellList[procI];
         }
@@ -98,7 +98,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
     {
         List<SLList<label> > procFaceList(nProcs_);
 
-        forAll (neighbour, facei)
+        forAll(neighbour, facei)
         {
             if (cellToProc_[owner[facei]] == cellToProc_[neighbour[facei]])
             {
@@ -117,7 +117,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
         List<SLList<label> > procPatchIndex(nProcs_);
 
-        forAll (neighbour, facei)
+        forAll(neighbour, facei)
         {
             if (cellToProc_[owner[facei]] != cellToProc_[neighbour[facei]])
             {
@@ -225,16 +225,16 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
         // for all processors, set the size of start index and patch size
         // lists to the number of patches in the mesh
-        forAll (procPatchSize_, procI)
+        forAll(procPatchSize_, procI)
         {
             procPatchSize_[procI].setSize(patches.size());
             procPatchStartIndex_[procI].setSize(patches.size());
         }
 
-        forAll (patches, patchi)
+        forAll(patches, patchi)
         {
             // Reset size and start index for all processors
-            forAll (procPatchSize_, procI)
+            forAll(procPatchSize_, procI)
             {
                 procPatchSize_[procI][patchi] = 0;
                 procPatchStartIndex_[procI][patchi] =
@@ -251,7 +251,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
                 const unallocLabelList& patchFaceCells =
                     patches[patchi].faceCells();
 
-                forAll (patchFaceCells, facei)
+                forAll(patchFaceCells, facei)
                 {
                     const label curProc = cellToProc_[patchFaceCells[facei]];
 
@@ -284,7 +284,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
                     cycOffset
                 );
 
-                forAll (firstFaceCells, facei)
+                forAll(firstFaceCells, facei)
                 {
                     if
                     (
@@ -439,7 +439,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
                 // Ordering in cyclic boundaries is important.
                 // Add the other half of cyclic faces for cyclic boundaries
                 // that remain on the processor
-                forAll (secondFaceCells, facei)
+                forAll(secondFaceCells, facei)
                 {
                     if
                     (
@@ -463,7 +463,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
         // Convert linked lists into normal lists
         // Add inter-processor boundaries and remember start indices
-        forAll (procFaceList, procI)
+        forAll(procFaceList, procI)
         {
             // Get internal and regular boundary processor faces
             SLList<label>& curProcFaces = procFaceList[procI];
@@ -604,7 +604,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
     // are omitted. For processor patches, set index to -1.
     // At the same time, filter the procPatchSize_ and procPatchStartIndex_
     // lists to exclude zero-size patches
-    forAll (procPatchSize_, procI)
+    forAll(procPatchSize_, procI)
     {
         // Make a local copy of old lists
         const labelList oldPatchSizes = procPatchSize_[procI];
@@ -628,7 +628,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
         label nPatches = 0;
 
-        forAll (oldPatchSizes, patchi)
+        forAll(oldPatchSizes, patchi)
         {
             if (!filterEmptyPatches || oldPatchSizes[patchi] > 0)
             {
@@ -646,7 +646,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
         curPatchSizes.setSize(nPatches);
         curPatchStarts.setSize(nPatches);
 
-        forAll (curProcessorPatchSizes, procPatchI)
+        forAll(curProcessorPatchSizes, procPatchI)
         {
             curBoundaryAddressing[nPatches] = -1;
 
@@ -662,19 +662,19 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
     // used for the processor. Collect the list of used points for the
     // processor.
 
-    forAll (procPointAddressing_, procI)
+    forAll(procPointAddressing_, procI)
     {
         boolList pointLabels(nPoints(), false);
 
         // Get reference to list of used faces
         const labelList& procFaceLabels = procFaceAddressing_[procI];
 
-        forAll (procFaceLabels, facei)
+        forAll(procFaceLabels, facei)
         {
             // Because of the turning index, some labels may be negative
             const labelList& facePoints = fcs[mag(procFaceLabels[facei]) - 1];
 
-            forAll (facePoints, pointi)
+            forAll(facePoints, pointi)
             {
                 // Mark the point as used
                 pointLabels[facePoints[pointi]] = true;
@@ -688,7 +688,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
         label nUsedPoints = 0;
 
-        forAll (pointLabels, pointi)
+        forAll(pointLabels, pointi)
         {
             if (pointLabels[pointi])
             {
@@ -734,7 +734,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
             // Reset the lookup list
             pointsUsage = 0;
 
-            forAll (curProcessorPatchStarts, patchi)
+            forAll(curProcessorPatchStarts, patchi)
             {
                 const label curStart = curProcessorPatchStarts[patchi];
                 const label curEnd = curStart + curProcessorPatchSizes[patchi];
@@ -753,7 +753,7 @@ void Foam::domainDecomposition::decomposeMesh(const bool filterEmptyPatches)
 
                     const face& f = fcs[curF];
 
-                    forAll (f, pointi)
+                    forAll(f, pointi)
                     {
                         if (pointsUsage[f[pointi]] == 0)
                         {
diff --git a/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposerDecomposeFields.C b/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposerDecomposeFields.C
index afea303fbc..eb2b34577a 100644
--- a/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposerDecomposeFields.C
+++ b/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposerDecomposeFields.C
@@ -43,7 +43,7 @@ Foam::fvFieldDecomposer::decomposeField
     // Create and map the patch field values
     PtrList<fvPatchField<Type> > patchFields(boundaryAddressing_.size());
 
-    forAll (boundaryAddressing_, patchi)
+    forAll(boundaryAddressing_, patchi)
     {
         if (boundaryAddressing_[patchi] >= 0)
         {
@@ -115,7 +115,7 @@ Foam::fvFieldDecomposer::decomposeField
             procMesh_.nInternalFaces()
         )
     );
-    forAll (mapAddr, i)
+    forAll(mapAddr, i)
     {
         mapAddr[i] -= 1;
     }
@@ -134,18 +134,18 @@ Foam::fvFieldDecomposer::decomposeField
     // (i.e. using slices)
     Field<Type> allFaceField(field.mesh().nFaces());
 
-    forAll (field.internalField(), i)
+    forAll(field.internalField(), i)
     {
         allFaceField[i] = field.internalField()[i];
     }
 
-    forAll (field.boundaryField(), patchi)
+    forAll(field.boundaryField(), patchi)
     {
         const Field<Type> & p = field.boundaryField()[patchi];
 
         const label patchStart = field.mesh().boundaryMesh()[patchi].start();
 
-        forAll (p, i)
+        forAll(p, i)
         {
             allFaceField[patchStart + i] = p[i];
         }
@@ -154,7 +154,7 @@ Foam::fvFieldDecomposer::decomposeField
     // Create and map the patch field values
     PtrList<fvsPatchField<Type> > patchFields(boundaryAddressing_.size());
 
-    forAll (boundaryAddressing_, patchi)
+    forAll(boundaryAddressing_, patchi)
     {
         if (boundaryAddressing_[patchi] >= 0)
         {
@@ -217,7 +217,7 @@ void Foam::fvFieldDecomposer::decomposeFields
     const PtrList<GeoField>& fields
 ) const
 {
-    forAll (fields, fieldI)
+    forAll(fields, fieldI)
     {
         decomposeField(fields[fieldI])().write();
     }
diff --git a/applications/utilities/parallelProcessing/decomposePar/lagrangianFieldDecomposerDecomposeFields.C b/applications/utilities/parallelProcessing/decomposePar/lagrangianFieldDecomposerDecomposeFields.C
index 36a5888d45..14ec7ff396 100644
--- a/applications/utilities/parallelProcessing/decomposePar/lagrangianFieldDecomposerDecomposeFields.C
+++ b/applications/utilities/parallelProcessing/decomposePar/lagrangianFieldDecomposerDecomposeFields.C
@@ -104,7 +104,7 @@ void Foam::lagrangianFieldDecomposer::decomposeFields
 {
     if (particleIndices_.size())
     {
-        forAll (fields, fieldI)
+        forAll(fields, fieldI)
         {
             decomposeField(cloudName, fields[fieldI])().write();
         }
diff --git a/applications/utilities/parallelProcessing/decomposePar/pointFieldDecomposerDecomposeFields.C b/applications/utilities/parallelProcessing/decomposePar/pointFieldDecomposerDecomposeFields.C
index 4f22046b39..b1c78e2486 100644
--- a/applications/utilities/parallelProcessing/decomposePar/pointFieldDecomposerDecomposeFields.C
+++ b/applications/utilities/parallelProcessing/decomposePar/pointFieldDecomposerDecomposeFields.C
@@ -48,7 +48,7 @@ Foam::pointFieldDecomposer::decomposeField
     );
 
     // Create and map the patch field values
-    forAll (boundaryAddressing_, patchi)
+    forAll(boundaryAddressing_, patchi)
     {
         if (patchFieldDecomposerPtrs_[patchi])
         {
@@ -117,7 +117,7 @@ void Foam::pointFieldDecomposer::decomposeFields
     const PtrList<GeoField>& fields
 ) const
 {
-    forAll (fields, fieldI)
+    forAll(fields, fieldI)
     {
         decomposeField(fields[fieldI])().write();
     }
diff --git a/applications/utilities/parallelProcessing/reconstructPar/checkFaceAddressingComp.H b/applications/utilities/parallelProcessing/reconstructPar/checkFaceAddressingComp.H
index 02e513a2ee..733d64f5ff 100644
--- a/applications/utilities/parallelProcessing/reconstructPar/checkFaceAddressingComp.H
+++ b/applications/utilities/parallelProcessing/reconstructPar/checkFaceAddressingComp.H
@@ -8,11 +8,11 @@
 
     PtrList<labelIOList>& faceProcAddressing = procMeshes.faceProcAddressing();
 
-    forAll (faceProcAddressing, procI)
+    forAll(faceProcAddressing, procI)
     {
         const labelList& curFaceAddr = faceProcAddressing[procI];
 
-        forAll (curFaceAddr, faceI)
+        forAll(curFaceAddr, faceI)
         {
             if (mag(curFaceAddr[faceI]) < minFaceIndex)
             {
@@ -33,11 +33,11 @@
             << "the current version fo decomposePar"
             << endl;
 
-        forAll (faceProcAddressing, procI)
+        forAll(faceProcAddressing, procI)
         {
             labelList& curFaceAddr = faceProcAddressing[procI];
 
-            forAll (curFaceAddr, faceI)
+            forAll(curFaceAddr, faceI)
             {
                 curFaceAddr[faceI] += sign(curFaceAddr[faceI]);
             }
diff --git a/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C b/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C
index 8def808e45..c412010cbd 100644
--- a/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C
+++ b/applications/utilities/parallelProcessing/reconstructPar/reconstructPar.C
@@ -81,7 +81,7 @@ int main(int argc, char *argv[])
     // Create the processor databases
     PtrList<Time> databases(nProcs);
 
-    forAll (databases, procI)
+    forAll(databases, procI)
     {
         databases.set
         (
@@ -118,7 +118,7 @@ int main(int argc, char *argv[])
     }
 
     // Set all times on processor meshes equal to reconstructed mesh
-    forAll (databases, procI)
+    forAll(databases, procI)
     {
         databases[procI].setTime(runTime.timeName(), runTime.timeIndex());
     }
@@ -132,7 +132,7 @@ int main(int argc, char *argv[])
 #   include "checkFaceAddressingComp.H"
 
     // Loop over all times
-    forAll (timeDirs, timeI)
+    forAll(timeDirs, timeI)
     {
         // Set time for global database
         runTime.setTime(timeDirs[timeI], timeI);
@@ -140,7 +140,7 @@ int main(int argc, char *argv[])
         Info<< "Time = " << runTime.timeName() << endl << endl;
 
         // Set time for all databases
-        forAll (databases, procI)
+        forAll(databases, procI)
         {
             databases[procI].setTime(timeDirs[timeI], timeI);
         }
@@ -247,7 +247,7 @@ int main(int argc, char *argv[])
             pointMesh pMesh(mesh);
             PtrList<pointMesh> pMeshes(procMeshes.meshes().size());
 
-            forAll (pMeshes, procI)
+            forAll(pMeshes, procI)
             {
                 pMeshes.set(procI, new pointMesh(procMeshes.meshes()[procI]));
             }
@@ -283,7 +283,7 @@ int main(int argc, char *argv[])
         {
             HashTable<IOobjectList> cloudObjects;
 
-            forAll (databases, procI)
+            forAll(databases, procI)
             {
                 fileNameList cloudDirs
                 (
@@ -294,7 +294,7 @@ int main(int argc, char *argv[])
                     )
                 );
 
-                forAll (cloudDirs, i)
+                forAll(cloudDirs, i)
                 {
                     // Check if we already have cloud objects for this cloudname
                     HashTable<IOobjectList>::const_iterator iter =
diff --git a/applications/utilities/parallelProcessing/reconstructParMesh/reconstructParMesh.C b/applications/utilities/parallelProcessing/reconstructParMesh/reconstructParMesh.C
index 04810a4e0c..a143356450 100644
--- a/applications/utilities/parallelProcessing/reconstructParMesh/reconstructParMesh.C
+++ b/applications/utilities/parallelProcessing/reconstructParMesh/reconstructParMesh.C
@@ -371,7 +371,7 @@ int main(int argc, char *argv[])
     // Read all databases.
     PtrList<Time> databases(nProcs);
 
-    forAll (databases, procI)
+    forAll(databases, procI)
     {
         Info<< "Reading database "
             << args.caseName()/fileName(word("processor") + name(procI))
diff --git a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/foamDataToFluent.C b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/foamDataToFluent.C
index 41df27a4c8..77781311f3 100644
--- a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/foamDataToFluent.C
+++ b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/foamDataToFluent.C
@@ -79,7 +79,7 @@ int main(int argc, char *argv[])
         // Writing number of faces
         label nFaces = mesh.nFaces();
 
-        forAll (mesh.boundary(), patchI)
+        forAll(mesh.boundary(), patchI)
         {
             nFaces += mesh.boundary()[patchI].size();
         }
diff --git a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentScalarField.C b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentScalarField.C
index 162915c329..cced2b0e32 100644
--- a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentScalarField.C
+++ b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentScalarField.C
@@ -59,7 +59,7 @@ void writeFluentField
         << "1 " << phiInternal.size() // Start and end of list
         << ")(" << endl;
 
-    forAll (phiInternal, cellI)
+    forAll(phiInternal, cellI)
     {
         stream << phiInternal[cellI] << endl;
     }
@@ -70,7 +70,7 @@ void writeFluentField
     label nWrittenFaces = phiInternal.size();
 
     // Writing boundary faces
-    forAll (phi.boundaryField(), patchI)
+    forAll(phi.boundaryField(), patchI)
     {
         if (isType<emptyFvPatchScalarField>(phi.boundaryField()[patchI]))
         {
@@ -95,7 +95,7 @@ void writeFluentField
 
             nWrittenFaces += emptyFaceCells.size();
 
-            forAll (emptyFaceCells, faceI)
+            forAll(emptyFaceCells, faceI)
             {
                 stream << phiInternal[emptyFaceCells[faceI]] << endl;
             }
@@ -123,7 +123,7 @@ void writeFluentField
 
             nWrittenFaces += patchPhi.size();
 
-            forAll (patchPhi, faceI)
+            forAll(patchPhi, faceI)
             {
                 stream << patchPhi[faceI] << endl;
             }
diff --git a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentVectorField.C b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentVectorField.C
index 188b58451c..b3285326da 100644
--- a/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentVectorField.C
+++ b/applications/utilities/postProcessing/dataConversion/foamDataToFluent/writeFluentVectorField.C
@@ -58,7 +58,7 @@ void writeFluentField
         << "1 " << phiInternal.size() // Start and end of list
         << ")(" << endl;
 
-    forAll (phiInternal, cellI)
+    forAll(phiInternal, cellI)
     {
         stream
             << phiInternal[cellI].x() << " "
@@ -73,7 +73,7 @@ void writeFluentField
     label nWrittenFaces = phiInternal.size();
 
     // Writing boundary faces
-    forAll (phi.boundaryField(), patchI)
+    forAll(phi.boundaryField(), patchI)
     {
         const vectorField& patchPhi = phi.boundaryField()[patchI];
 
@@ -90,7 +90,7 @@ void writeFluentField
 
         nWrittenFaces += patchPhi.size();
 
-        forAll (patchPhi, faceI)
+        forAll(patchPhi, faceI)
         {
             stream
                 << patchPhi[faceI].x() << " "
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightField.C b/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightField.C
index 90a597e578..a7e1e4dbd7 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightField.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightField.C
@@ -199,7 +199,7 @@ void writeAllFaceData
             {
                 writeData(map(pf, prims, cmpt), ensightFile);
 
-                forAll (patchProcessors, i)
+                forAll(patchProcessors, i)
                 {
                     if (patchProcessors[i] != 0)
                     {
@@ -245,7 +245,7 @@ void writeAllFaceDataBinary
             {
                 writeEnsDataBinary(map(pf, prims, cmpt), ensightFile);
 
-                forAll (patchProcessors, i)
+                forAll(patchProcessors, i)
                 {
                     if (patchProcessors[i] != 0)
                     {
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightMesh.C b/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightMesh.C
index 4796f4b4f1..c58cd63e85 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightMesh.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsight/ensightMesh.C
@@ -112,7 +112,7 @@ Foam::ensightMesh::ensightMesh
 
         allPatchProcs_.setSize(allPatchNames_.size());
 
-        forAll (allPatchProcs_, patchi)
+        forAll(allPatchProcs_, patchi)
         {
             if (mesh_.boundary()[patchi].size())
             {
@@ -132,7 +132,7 @@ Foam::ensightMesh::ensightMesh
                 patchNameList = allPatchNames_;
             }
 
-            forAll (patchNameList, i)
+            forAll(patchNameList, i)
             {
                 patchNames_.insert(patchNameList[i]);
             }
@@ -217,7 +217,7 @@ Foam::ensightMesh::ensightMesh
 
     if (!args.optionFound("noPatches"))
     {
-        forAll (mesh.boundary(), patchi)
+        forAll(mesh.boundary(), patchi)
         {
             if (mesh.boundary()[patchi].size())
             {
@@ -945,7 +945,7 @@ Foam::faceList Foam::ensightMesh::map
 {
     faceList ppf(prims.size());
 
-    forAll (prims, i)
+    forAll(prims, i)
     {
         ppf[i] = patchFaces[prims[i]];
     }
@@ -981,7 +981,7 @@ void Foam::ensightMesh::writeAllFacePrims
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1067,7 +1067,7 @@ void Foam::ensightMesh::writeAllNSided
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1102,7 +1102,7 @@ void Foam::ensightMesh::writeAllNSided
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1191,7 +1191,7 @@ void Foam::ensightMesh::writeAllNSidedBinary
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1226,7 +1226,7 @@ void Foam::ensightMesh::writeAllNSidedBinary
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1280,7 +1280,7 @@ void Foam::ensightMesh::writeAllFacePrimsBinary
                 );
             }
 
-            forAll (patchProcessors, i)
+            forAll(patchProcessors, i)
             {
                 if (patchProcessors[i] != 0)
                 {
@@ -1531,7 +1531,7 @@ void Foam::ensightMesh::writeAscii
 
                         patchPointOffsets = 0;
 
-                        forAll (patchProcessors, i)
+                        forAll(patchProcessors, i)
                         {
                             if (patchProcessors[i] != 0)
                             {
@@ -1820,7 +1820,7 @@ void Foam::ensightMesh::writeBinary
                         patchPointOffsets = 0;
 
 
-                        forAll (patchProcessors, i)
+                        forAll(patchProcessors, i)
                         {
                             if (patchProcessors[i] != 0)
                             {
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputCase.H b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputCase.H
index 0e3530b9c0..3ccdbbe91b 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputCase.H
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputCase.H
@@ -144,7 +144,7 @@ if (fieldTimesUsed.size())
         << "filename numbers:" << nl;
 
     label count = 0;
-    forAll (fieldTimesUsed, i)
+    forAll(fieldTimesUsed, i)
     {
         caseFile
             << " " << setw(12) << fieldTimesUsed[i];
@@ -159,7 +159,7 @@ if (fieldTimesUsed.size())
         << nl << "time values:" << nl;
 
     count = 0;
-    forAll (fieldTimesUsed, i)
+    forAll(fieldTimesUsed, i)
     {
         caseFile
             << " " << setw(12)
@@ -188,7 +188,7 @@ forAllConstIter(HashTable<DynamicList<label> >, cloudTimesUsed, cloudIter)
             << "filename numbers:" << nl;
 
         label count = 0;
-        forAll (timesUsed, i)
+        forAll(timesUsed, i)
         {
             caseFile
                 << " " << setw(12) << timesUsed[i];
@@ -203,7 +203,7 @@ forAllConstIter(HashTable<DynamicList<label> >, cloudTimesUsed, cloudIter)
             << nl << "time values:" << nl;
 
         count = 0;
-        forAll (timesUsed, i)
+        forAll(timesUsed, i)
         {
             caseFile
                 << " " << setw(12)
diff --git a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputFunctions.C b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputFunctions.C
index 4f7f8fdeda..c6008b0e12 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputFunctions.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToEnsightParts/ensightOutputFunctions.C
@@ -117,12 +117,12 @@ void Foam::ensightParticlePositions
     // binary write is Ensight6 - first ids, then positions
     if (format == IOstream::BINARY)
     {
-        forAll (parcels, i)
+        forAll(parcels, i)
         {
             os.write(i+1);
         }
 
-        forAllIter(Cloud<passiveParticle>, parcels, elmnt)
+        forAllConstIter(Cloud<passiveParticle>, parcels, elmnt)
         {
             const vector& p = elmnt().position();
 
@@ -135,7 +135,7 @@ void Foam::ensightParticlePositions
     {
         label nParcels = 0;
 
-        forAllIter(Cloud<passiveParticle>, parcels, elmnt)
+        forAllConstIter(Cloud<passiveParticle>, parcels, elmnt)
         {
             const vector& p = elmnt().position();
 
diff --git a/applications/utilities/postProcessing/dataConversion/foamToFieldview9/fieldviewTopology.C b/applications/utilities/postProcessing/dataConversion/foamToFieldview9/fieldviewTopology.C
index 473ab6f7af..cb929f9749 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToFieldview9/fieldviewTopology.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToFieldview9/fieldviewTopology.C
@@ -52,7 +52,7 @@ Foam::labelList Foam::fieldviewTopology::calcFaceAddressing
     const label cellI
 )
 {
-    // return value. 
+    // return value.
     labelList shapeToMesh(shape.nFaces(), -1);
 
     const faceList modelFaces(shape.faces());
@@ -122,10 +122,10 @@ Foam::fieldviewTopology::fieldviewTopology
 
     if (setWallInfo)
     {
-        forAll (mesh.boundaryMesh(), patchI)
+        forAll(mesh.boundaryMesh(), patchI)
         {
             const polyPatch& currPatch = mesh.boundaryMesh()[patchI];
-            if 
+            if
             (
                 isA<wallPolyPatch>(currPatch)
              || isA<symmetryPolyPatch>(currPatch)
@@ -424,7 +424,7 @@ Foam::fieldviewTopology::fieldviewTopology
                 faceLabels[labelI++] = patchFace[1] + 1;
                 faceLabels[labelI++] = patchFace[2] + 1;
                 faceLabels[labelI++] = patchFace[3] + 1;
-            } 
+            }
         }
 
         faceLabels.setSize(labelI);
diff --git a/applications/utilities/postProcessing/dataConversion/foamToFieldview9/foamToFieldview9.C b/applications/utilities/postProcessing/dataConversion/foamToFieldview9/foamToFieldview9.C
index 1e3b781a7e..9f06a7fb6c 100644
--- a/applications/utilities/postProcessing/dataConversion/foamToFieldview9/foamToFieldview9.C
+++ b/applications/utilities/postProcessing/dataConversion/foamToFieldview9/foamToFieldview9.C
@@ -379,7 +379,7 @@ int main(int argc, char *argv[])
         // num patches
         writeInt(fvFile, mesh.boundary().size());
 
-        forAll (mesh.boundary(), patchI)
+        forAll(mesh.boundary(), patchI)
         {
             const fvPatch& currPatch = mesh.boundary()[patchI];
 
@@ -694,7 +694,7 @@ int main(int argc, char *argv[])
             {
                 const volScalarField& vsf = *volFieldPtrs[fieldI];
 
-                forAll (mesh.boundary(), patchI)
+                forAll(mesh.boundary(), patchI)
                 {
                     writeFaceData
                     (
@@ -709,7 +709,7 @@ int main(int argc, char *argv[])
             }
             else
             {
-                forAll (mesh.boundaryMesh(), patchI)
+                forAll(mesh.boundaryMesh(), patchI)
                 {
                     // Dummy value.
                     floatField fField
@@ -735,7 +735,7 @@ int main(int argc, char *argv[])
             {
                 const surfaceScalarField& ssf = *surfFieldPtrs[fieldI];
 
-                forAll (mesh.boundary(), patchI)
+                forAll(mesh.boundary(), patchI)
                 {
                     writeFaceData
                     (
@@ -750,7 +750,7 @@ int main(int argc, char *argv[])
             }
             else
             {
-                forAll (mesh.boundaryMesh(), patchI)
+                forAll(mesh.boundaryMesh(), patchI)
                 {
                     // Dummy value.
                     floatField fField
@@ -784,7 +784,7 @@ int main(int argc, char *argv[])
                 const volScalarField& vsf = *volFieldPtrs[fieldI];
 
                 // All non-empty patches
-                forAll (mesh.boundary(), patchI)
+                forAll(mesh.boundary(), patchI)
                 {
                     writeFaceData
                     (
@@ -799,7 +799,7 @@ int main(int argc, char *argv[])
             }
             else
             {
-                forAll (mesh.boundary(), patchI)
+                forAll(mesh.boundary(), patchI)
                 {
                     // Dummy value.
                     floatField fField(topo.nPolyFaces()[patchI], 0.0);
@@ -836,7 +836,7 @@ int main(int argc, char *argv[])
             }
             else
             {
-                forAll (mesh.boundaryMesh(), patchI)
+                forAll(mesh.boundaryMesh(), patchI)
                 {
                     // Dummy value.
                     floatField fField
diff --git a/applications/utilities/postProcessing/dataConversion/smapToFoam/smapToFoam.C b/applications/utilities/postProcessing/dataConversion/smapToFoam/smapToFoam.C
index 1a29e87652..ee450c1200 100644
--- a/applications/utilities/postProcessing/dataConversion/smapToFoam/smapToFoam.C
+++ b/applications/utilities/postProcessing/dataConversion/smapToFoam/smapToFoam.C
@@ -49,7 +49,7 @@ int main(int argc, char *argv[])
 
     fileNameList fieldNames = readDir(runTime.timePath(), fileName::FILE);
     dictionary fieldNameDict;
-    forAll (fieldNames, i)
+    forAll(fieldNames, i)
     {
         fieldNameDict.add(fieldNames[i], word(fieldNames[i]));
     }
@@ -168,7 +168,7 @@ int main(int argc, char *argv[])
 
         label cell;
         scalar value;
-        forAll (mesh.cells(), celli)
+        forAll(mesh.cells(), celli)
         {
             if (celli > 0)
             {
diff --git a/applications/utilities/postProcessing/velocityField/streamFunction/streamFunction.C b/applications/utilities/postProcessing/velocityField/streamFunction/streamFunction.C
index 2b4e3ce58d..113b1ff612 100644
--- a/applications/utilities/postProcessing/velocityField/streamFunction/streamFunction.C
+++ b/applications/utilities/postProcessing/velocityField/streamFunction/streamFunction.C
@@ -101,7 +101,7 @@ int main(int argc, char *argv[])
             );
 
             labelList visitedPoint(mesh.nPoints());
-            forAll (visitedPoint, pointI)
+            forAll(visitedPoint, pointI)
             {
                 visitedPoint[pointI] = 0;
             }
@@ -128,13 +128,13 @@ int main(int argc, char *argv[])
             {
                 found = false;
 
-                forAll (patches, patchI)
+                forAll(patches, patchI)
                 {
                     const primitivePatch& bouFaces = patches[patchI];
 
                     if (!isType<emptyPolyPatch>(patches[patchI]))
                     {
-                        forAll (bouFaces, faceI)
+                        forAll(bouFaces, faceI)
                         {
                             if
                             (
@@ -147,7 +147,7 @@ int main(int argc, char *argv[])
                                 // Zero flux face found
                                 found = true;
 
-                                forAll (zeroPoints, pointI)
+                                forAll(zeroPoints, pointI)
                                 {
                                     if (visitedPoint[zeroPoints[pointI]] == 1)
                                     {
@@ -161,7 +161,7 @@ int main(int argc, char *argv[])
                                     Info<< "Zero face: patch: " << patchI
                                         << "    face: " << faceI << endl;
 
-                                    forAll (zeroPoints, pointI)
+                                    forAll(zeroPoints, pointI)
                                     {
                                         streamFunction[zeroPoints[pointI]] = 0;
                                         visitedPoint[zeroPoints[pointI]] = 1;
@@ -185,13 +185,13 @@ int main(int argc, char *argv[])
 
                     const cellList& c = mesh.cells();
 
-                    forAll (c, cI)
+                    forAll(c, cI)
                     {
                         labelList zeroPoints = c[cI].labels(mesh.faces());
 
                         bool found = true;
 
-                        forAll (zeroPoints, pointI)
+                        forAll(zeroPoints, pointI)
                         {
                             if (visitedPoint[zeroPoints[pointI]] == 1)
                             {
@@ -202,7 +202,7 @@ int main(int argc, char *argv[])
 
                         if (found)
                         {
-                            forAll (zeroPoints, pointI)
+                            forAll(zeroPoints, pointI)
                             {
                                 streamFunction[zeroPoints[pointI]] = 0.0;
                                 visitedPoint[zeroPoints[pointI]] = 1;
@@ -242,7 +242,7 @@ int main(int argc, char *argv[])
                          scalar currentBStream = 0.0;
                          vector currentBStreamPoint(0, 0, 0);
 
-                         forAll (curBPoints, pointI)
+                         forAll(curBPoints, pointI)
                          {
                              // Check if the point has been visited
                              if (visitedPoint[curBPoints[pointI]] == 1)
@@ -262,7 +262,7 @@ int main(int argc, char *argv[])
                          if (bPointFound)
                          {
                              // Sort out other points on the face
-                             forAll (curBPoints, pointI)
+                             forAll(curBPoints, pointI)
                              {
                                  // Check if the point has been visited
                                  if (visitedPoint[curBPoints[pointI]] == 0)
@@ -367,7 +367,7 @@ int main(int argc, char *argv[])
                          scalar currentStream = 0.0;
                          point currentStreamPoint(0, 0, 0);
 
-                         forAll (curPoints, pointI)
+                         forAll(curPoints, pointI)
                          {
                              // Check if the point has been visited
                              if (visitedPoint[curPoints[pointI]] == 1)
@@ -386,7 +386,7 @@ int main(int argc, char *argv[])
                          if (pointFound)
                          {
                              // Sort out other points on the face
-                             forAll (curPoints, pointI)
+                             forAll(curPoints, pointI)
                              {
                                  // Check if the point has been visited
                                  if (visitedPoint[curPoints[pointI]] == 0)
diff --git a/applications/utilities/preProcessing/mapFields/mapFields.C b/applications/utilities/preProcessing/mapFields/mapFields.C
index ccaf04d4a5..c93f563e85 100644
--- a/applications/utilities/preProcessing/mapFields/mapFields.C
+++ b/applications/utilities/preProcessing/mapFields/mapFields.C
@@ -194,12 +194,12 @@ wordList addProcessorPatches
 {
     // Add the processor patches to the cutting list
     HashTable<label> cuttingPatchTable;
-    forAll (cuttingPatches, i)
+    forAll(cuttingPatches, i)
     {
         cuttingPatchTable.insert(cuttingPatches[i], i);
     }
 
-    forAll (meshTarget.boundary(), patchi)
+    forAll(meshTarget.boundary(), patchi)
     {
         if (isA<processorFvPatch>(meshTarget.boundary()[patchi]))
         {
diff --git a/applications/utilities/preProcessing/mapFields/mapLagrangian.C b/applications/utilities/preProcessing/mapFields/mapLagrangian.C
index 575c6fa396..591fb8ff24 100644
--- a/applications/utilities/preProcessing/mapFields/mapLagrangian.C
+++ b/applications/utilities/preProcessing/mapFields/mapLagrangian.C
@@ -223,7 +223,7 @@ void mapLagrangian(const meshToMesh& meshToMeshInterp)
                         {
                             unmappedSource.erase(sourceParticleI);
                             addParticles.append(sourceParticleI);
-                            iter().cell()=targetCell;
+                            iter().cell() = targetCell;
                             targetParcels.addParticle
                             (
                                 sourceParcels.remove(&iter())
diff --git a/applications/utilities/surface/surfaceSubset/surfaceSubset.C b/applications/utilities/surface/surfaceSubset/surfaceSubset.C
index e32677fe77..dd8a6ede00 100644
--- a/applications/utilities/surface/surfaceSubset/surfaceSubset.C
+++ b/applications/utilities/surface/surfaceSubset/surfaceSubset.C
@@ -121,7 +121,7 @@ int main(int argc, char *argv[])
 
         // pick up cells sharing the point
 
-        forAll (markedPoints, pointI)
+        forAll(markedPoints, pointI)
         {
             if
             (
@@ -140,7 +140,7 @@ int main(int argc, char *argv[])
             const labelList& curFaces =
                 surf1.pointFaces()[markedPoints[pointI]];
 
-            forAll (curFaces, i)
+            forAll(curFaces, i)
             {
                 facesToSubset[curFaces[i]] =  true;
             }
@@ -159,7 +159,7 @@ int main(int argc, char *argv[])
 
         // pick up cells sharing the edge
 
-        forAll (markedEdges, edgeI)
+        forAll(markedEdges, edgeI)
         {
             if
             (
@@ -177,7 +177,7 @@ int main(int argc, char *argv[])
 
             const labelList& curFaces = surf1.edgeFaces()[markedEdges[edgeI]];
 
-            forAll (curFaces, i)
+            forAll(curFaces, i)
             {
                 facesToSubset[curFaces[i]] =  true;
             }
@@ -298,7 +298,7 @@ int main(int argc, char *argv[])
         Info<< "Found " << markedFaces.size() << " marked face(s)." << endl;
 
         // Check and mark faces to pick up
-        forAll (markedFaces, faceI)
+        forAll(markedFaces, faceI)
         {
             if
             (
@@ -322,7 +322,7 @@ int main(int argc, char *argv[])
                 const labelList& curFaces =
                     surf1.faceFaces()[markedFaces[faceI]];
 
-                forAll (curFaces, i)
+                forAll(curFaces, i)
                 {
                     label faceI = curFaces[i];