remove trailing space from some files

applications/solvers/combustion/PDRFoam/XiEqns

@@ -22,7 +22,7 @@
         volScalarField GEta = GEtaCoef/tauEta;
         volScalarField XiEqEta = 1.0 + XiCoef*sqrt(up/(Su + SuMin))*Reta;
 
-        volScalarField R = 
+        volScalarField R =
             GEta*XiEqEta/(XiEqEta - 0.999) + GIn*XiIn/(XiIn - 0.999);
 
         volScalarField XiEqStar = R/(R - GEta - GIn);
@@ -42,7 +42,7 @@
         volScalarField GEta = GEtaCoef/tauEta;
         volScalarField XiEqEta = 1.0 + XiCoef*sqrt(up/(Su + SuMin))*Reta;
 
-        volScalarField R = 
+        volScalarField R =
             GEta*XiEqEta/(XiEqEta - 0.999) + GIn*XiIn/(XiIn - 0.999);
 
         volScalarField XiEqStar = R/(R - GEta - GIn);

applications/solvers/combustion/PDRFoam/XiModels/algebraic/algebraic.C

@@ -91,7 +91,7 @@ bool Foam::XiModels::algebraic::read(const dictionary& XiProperties)
     XiModel::read(XiProperties);
 
     XiModelCoeffs_.lookup("XiShapeCoef") >> XiShapeCoef;
-    
+
     return true;
 }
 

applications/solvers/combustion/PDRFoam/laminarFlameSpeed/SCOPE/SCOPELaminarFlameSpeed.C

@@ -117,7 +117,7 @@ inline Foam::scalar Foam::laminarFlameSpeedModels::SCOPE::polyPhi
 {
     scalar x = phi - 1.0;
 
-    return 
+    return
         a[0]
        *(
            scalar(1)

applications/solvers/combustion/PDRFoam/setDeltaT.H

@@ -45,7 +45,7 @@ if (adjustTimeStep)
             maxDeltaT
         )
     );
-    
+
     Info<< "deltaT = " <<  runTime.deltaTValue() << endl;
 }
 

applications/solvers/compressible/rhoCentralFoam/Allwmake

@@ -2,7 +2,7 @@
 cd ${0%/*} || exit 1    # run from this directory
 set -x
 
-wmake libso BCs 
+wmake libso BCs
 wmake
 
 # ----------------------------------------------------------------- end-of-file

applications/solvers/compressible/rhoSonicFoam/rhoSonicFoam.C

@@ -84,7 +84,7 @@ int main(int argc, char *argv[])
         (
             fvm::ddt(rhoU)
           + fvm::div(phiv, rhoU)
-         == 
+         ==
           - fvc::grad(p)
         );
 

applications/solvers/compressible/rhopSonicFoam/Allwmake

@@ -2,7 +2,7 @@
 cd ${0%/*} || exit 1    # run from this directory
 set -x
 
-wmake libso BCs 
+wmake libso BCs
 wmake
 
 # ----------------------------------------------------------------- end-of-file

applications/solvers/compressible/rhopSonicFoam/BCs/p/inviscidWallPFvPatchScalarField.C

@@ -131,7 +131,7 @@ void inviscidWallPFvPatchScalarField::updateCoeffs()
     const fvPatchField<vector>& rhoUp =
         patch().lookupPatchField<volVectorField, vector>("rhoU");
 
-    const fvsPatchField<scalar>& phip = 
+    const fvsPatchField<scalar>& phip =
         patch().lookupPatchField<surfaceScalarField, scalar>("phi");
 
     const fvsPatchField<scalar>& rAp =
@@ -147,7 +147,7 @@ void inviscidWallPFvPatchScalarField::updateCoeffs()
 void inviscidWallPFvPatchScalarField::write(Ostream& os) const
 {
     fixedGradientFvPatchScalarField::write(os);
-    os.writeKeyword("fluxFraction") 
+    os.writeKeyword("fluxFraction")
         << fluxFraction_ << token::END_STATEMENT << nl;
     writeEntry("value", os);
 }

applications/solvers/compressible/rhopSonicFoam/BCs/rhoE/mixedRhoEFvPatchScalarField.C

@@ -145,7 +145,7 @@ void mixedRhoEFvPatchScalarField::updateCoeffs()
 
     const volScalarField& T = db().lookupObject<volScalarField>("T");
     const label patchi = patch().index();
-    fvPatchScalarField& Tp = 
+    fvPatchScalarField& Tp =
         const_cast<fvPatchScalarField&>(T.boundaryField()[patchi]);
 
     Tp.evaluate();
@@ -164,7 +164,7 @@ void mixedRhoEFvPatchScalarField::updateCoeffs()
     refGrad() =
         rhop*Cv.value()*Tp.snGrad()
       + (
-            refValue() 
+            refValue()
           - (0.5*rhop.patchInternalField()*
                 magSqr(rhoUp.patchInternalField()/rhop.patchInternalField()))
         )*patch().deltaCoeffs();

applications/solvers/compressible/rhopSonicFoam/BCs/rhoU/fixedRhoUFvPatchVectorField.H

@@ -66,7 +66,7 @@ public:
             const dictionary&
         );
 
-        //- Construct by mapping given fixedRhoUFvPatchVectorField 
+        //- Construct by mapping given fixedRhoUFvPatchVectorField
         //  onto a new patch
         fixedRhoUFvPatchVectorField
         (

applications/solvers/electromagnetics/electrostaticFoam/Make/files

-electrostaticFoam.C  
+electrostaticFoam.C
 
 EXE = $(FOAM_APPBIN)/electrostaticFoam

applications/solvers/incompressible/icoFoam/Make/files

-icoFoam.C  
+icoFoam.C
 
 EXE = $(FOAM_APPBIN)/icoFoam

applications/solvers/incompressible/nonNewtonianIcoFoam/Make/files

-nonNewtonianIcoFoam.C  
+nonNewtonianIcoFoam.C
 
 EXE = $(FOAM_APPBIN)/nonNewtonianIcoFoam

applications/solvers/incompressible/simpleFoam/Make/files

-simpleFoam.C 
+simpleFoam.C
 
 EXE = $(FOAM_APPBIN)/simpleFoam

applications/solvers/multiphase/bubbleFoam/UEqns.H

@@ -41,7 +41,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
        - fvc::interpolate(nuEffb)
         *mesh.magSf()*fvc::snGrad(beta)/fvc::interpolate(beta + scalar(0.001));
 
-    UbEqn = 
+    UbEqn =
     (
         (scalar(1) + Cvm*rhob*alpha/rhob)*
         (

applications/solvers/multiphase/cavitatingFoam/CourantNo.H

@@ -36,7 +36,7 @@ scalar acousticCoNum = 0.0;
 
 if (mesh.nInternalFaces())
 {
-    surfaceScalarField SfUfbyDelta = 
+    surfaceScalarField SfUfbyDelta =
         mesh.surfaceInterpolation::deltaCoeffs()*mag(phiv);
 
     CoNum = max(SfUfbyDelta/mesh.magSf())

applications/solvers/multiphase/cavitatingFoam/continuityErrs.H

@@ -3,12 +3,12 @@
 
     dimensionedScalar totalMass = fvc::domainIntegrate(rho);
 
-    scalar sumLocalContErr = 
+    scalar sumLocalContErr =
     (
         fvc::domainIntegrate(mag(rho - thermoRho))/totalMass
     ).value();
 
-    scalar globalContErr = 
+    scalar globalContErr =
     (
         fvc::domainIntegrate(rho - thermoRho)/totalMass
     ).value();

applications/solvers/multiphase/cavitatingFoam/setDeltaT.H

@@ -47,7 +47,7 @@ if (adjustTimeStep)
             maxDeltaT
         )
     );
-    
+
     Info<< "deltaT = " <<  runTime.deltaTValue() << endl;
 }
 

applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/alphaContactAngle/alphaContactAngleFvPatchScalarField.C

@@ -63,7 +63,7 @@ Ostream& operator<<
 )
 {
     os  << tp.theta0_ << token::SPACE
-        << tp.uTheta_ << token::SPACE 
+        << tp.uTheta_ << token::SPACE
         << tp.thetaA_ << token::SPACE
         << tp.thetaR_;
 

applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/phase/phase.C

@@ -28,7 +28,7 @@ License
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-Foam::phase::phase 
+Foam::phase::phase
 (
     const word& name,
     const dictionary& phaseDict,

applications/solvers/multiphase/settlingFoam/correctViscosity.H

 {
-    mul = muc + 
+    mul = muc +
         plasticViscosity
         (
             plasticViscosityCoeff,

applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H

@@ -56,7 +56,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
             -fvc::interpolate(nuEffb)*mesh.magSf()*fvc::snGrad(beta)
             /fvc::interpolate(beta + scalar(0.001));
 
-        UbEqn = 
+        UbEqn =
         (
             (scalar(1) + Cvm*rhob*alpha/rhob)*
             (

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/conductivityModel/conductivityModel/newConductivityModel.C

@@ -35,9 +35,9 @@ Foam::autoPtr<Foam::conductivityModel> Foam::conductivityModel::New
 {
     word conductivityModelType(dict.lookup("conductivityModel"));
 
-    Info<< "Selecting conductivityModel " 
+    Info<< "Selecting conductivityModel "
         << conductivityModelType << endl;
-        
+
     dictionaryConstructorTable::iterator cstrIter =
         dictionaryConstructorTablePtr_->find(conductivityModelType);
 

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C

@@ -73,7 +73,7 @@ frictionalPressure
 ) const
 {
 
-    return 
+    return
         Fr*pow(max(alpha - alphaMinFriction, scalar(0)), eta)
        /pow(max(alphaMax - alpha, scalar(5.0e-2)), p);
 }
@@ -104,9 +104,9 @@ Foam::tmp<Foam::volScalarField> Foam::JohnsonJacksonFrictionalStress::muf
     const dimensionedScalar& alphaMax,
     const volScalarField& pf,
     const volTensorField& D,
-    const dimensionedScalar& phi            
+    const dimensionedScalar& phi
 ) const
-{    
+{
     return dimensionedScalar("0.5", dimTime, 0.5)*pf*sin(phi);
 }
 

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/frictionalStressModel/frictionalStressModel/newFrictionalStressModel.C

@@ -35,9 +35,9 @@ Foam::autoPtr<Foam::frictionalStressModel> Foam::frictionalStressModel::New
 {
     word frictionalStressModelType(dict.lookup("frictionalStressModel"));
 
-    Info<< "Selecting frictionalStressModel " 
+    Info<< "Selecting frictionalStressModel "
         << frictionalStressModelType << endl;
-        
+
     dictionaryConstructorTable::iterator cstrIter =
         dictionaryConstructorTablePtr_->find(frictionalStressModelType);
 

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/granularPressureModel/granularPressureModel/newGranularPressureModel.C

@@ -35,9 +35,9 @@ Foam::autoPtr<Foam::granularPressureModel> Foam::granularPressureModel::New
 {
     word granularPressureModelType(dict.lookup("granularPressureModel"));
 
-    Info<< "Selecting granularPressureModel " 
+    Info<< "Selecting granularPressureModel "
         << granularPressureModelType << endl;
-        
+
     dictionaryConstructorTable::iterator cstrIter =
         dictionaryConstructorTablePtr_->find(granularPressureModelType);
 

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/radialModel/CarnahanStarling/CarnahanStarlingRadial.C

@@ -65,7 +65,7 @@ Foam::tmp<Foam::volScalarField> Foam::CarnahanStarlingRadial::g0
 ) const
 {
 
-    return 
+    return
         1.0/(1.0 - alpha)
       + 3.0*alpha/(2.0*sqr(1.0 - alpha))
       + sqr(alpha)/(2.0*pow(1.0 - alpha, 3));
@@ -78,8 +78,8 @@ Foam::tmp<Foam::volScalarField> Foam::CarnahanStarlingRadial::g0prime
     const dimensionedScalar& alphaMax
 ) const
 {
-    return 
-        - alpha/sqr(1.0 - alpha) 
+    return
+        - alpha/sqr(1.0 - alpha)
         + (3.0*(1.0 - alpha) + 6.0*sqr(alpha))/(2.0*(1.0 - alpha))
         + (2.0*alpha*(1.0 - alpha) + 3.0*pow(alpha, 3))
          /(2.0*pow(1.0 - alpha, 4));

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/radialModel/Gidaspow/GidaspowRadial.C

@@ -74,7 +74,7 @@ Foam::tmp<Foam::volScalarField> Foam::GidaspowRadial::g0prime
     const dimensionedScalar& alphaMax
 ) const
 {
-    return 
+    return
         (-1.0/5.0)*pow(alpha/alphaMax, -2.0/3.0)
        /(alphaMax*sqr(1.0 - pow(alpha/alphaMax, 1.0/3.0)));
 }

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/radialModel/SinclairJackson/SinclairJacksonRadial.C

@@ -74,7 +74,7 @@ Foam::tmp<Foam::volScalarField> Foam::SinclairJacksonRadial::g0prime
     const dimensionedScalar& alphaMax
 ) const
 {
-    return 
+    return
         (-1.0/3.0)*pow(alpha/alphaMax, -2.0/3.0)
        /(alphaMax*sqr(1.0 - pow(alpha/alphaMax, 1.0/3.0)));
 }

applications/solvers/multiphase/twoPhaseEulerFoam/kineticTheoryModels/radialModel/radialModel/newRadialModel.C

@@ -35,9 +35,9 @@ Foam::autoPtr<Foam::radialModel> Foam::radialModel::New
 {
     word radialModelType(dict.lookup("radialModel"));
 
-    Info<< "Selecting radialModel " 
+    Info<< "Selecting radialModel "
         << radialModelType << endl;
-        
+
     dictionaryConstructorTable::iterator cstrIter =
         dictionaryConstructorTablePtr_->find(radialModelType);
 

applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H

@@ -42,7 +42,7 @@
 
     phib = (fvc::interpolate(Ub) & mesh.Sf()) + fvc::ddtPhiCorr(rUbA, Ub, phib)
          + phiDragb;
-    
+
     phi = alphaf*phia + betaf*phib;
 
     surfaceScalarField Dp("(rho*(1|A(U)))", alphaf*rUaAf/rhoa + betaf*rUbAf/rhob);

applications/solvers/stressAnalysis/solidDisplacementFoam/tractionDisplacement/tractionDisplacementFvPatchVectorField.C

@@ -177,7 +177,7 @@ void tractionDisplacementFvPatchVectorField::updateCoeffs()
     gradient() =
     (
         (traction_ + pressure_*n)/rho.value()
-      + twoMuLambda*fvPatchField<vector>::snGrad() - (n & sigmaD) 
+      + twoMuLambda*fvPatchField<vector>::snGrad() - (n & sigmaD)
     )/twoMuLambda;
 
     Switch thermalStress(thermalProperties.lookup("thermalStress"));
@@ -187,7 +187,7 @@ void tractionDisplacementFvPatchVectorField::updateCoeffs()
         dimensionedScalar alpha(thermalProperties.lookup("alpha"));
         dimensionedScalar threeKalpha = threeK*alpha;
 
-        const fvPatchField<scalar>& T = 
+        const fvPatchField<scalar>& T =
             patch().lookupPatchField<volScalarField, scalar>("T");
 
         gradient() += n*threeKalpha.value()*T/twoMuLambda;

applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/solidEquilibriumDisplacementFoam.C

@@ -72,7 +72,7 @@ int main(int argc, char *argv[])
         {
             volTensorField gradDcorr = fvc::grad(Dcorr);
 
-            sigmaExp = 
+            sigmaExp =
                 (lambda - mu)*gradDcorr + mu*gradDcorr.T()
               + (lambda*I)*tr(gradDcorr);
 

applications/solvers/stressAnalysis/solidEquilibriumDisplacementFoam/tractionDisplacementCorrection/tractionDisplacementCorrectionFvPatchVectorField.C

@@ -171,7 +171,7 @@ void tractionDisplacementCorrectionFvPatchVectorField::updateCoeffs()
     const fvPatchField<symmTensor>& sigmaD =
         patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD");
 
-    const fvPatchField<tensor>& sigmaExp = 
+    const fvPatchField<tensor>& sigmaExp =
         patch().lookupPatchField<volTensorField, tensor>("sigmaExp");
 
     gradient() =

applications/utilities/mesh/advanced/selectCells/edgeStats.H

@@ -67,7 +67,7 @@ class edgeStats
 
         //- If 2d get component of normal dir.
         direction getNormalDir(const twoDPointCorrector*) const;
-        
+
 
         //- Disallow default bitwise copy construct
         edgeStats(const edgeStats&);

applications/utilities/mesh/conversion/fluentMeshToFoam/README.txt

@@ -2,35 +2,35 @@ Notes for fluentMeshToFoam with zone preservation
 #################################################
 
 1. New option added:
-    - writeSets: 
+    - writeSets:
     Writes all Fluent boundaries faceSets preserving Fluent names
     Writes all Fluent regions to cellSets preserving Fluent names
     lines: 1375 - 1393 & 1673 - 1741
     sets are useful for post-processing using foamToVTK with the "-faceSet
     <name>" and "-cellSet <name>" options.
-    
-    - writeZones: 
+
+    - writeZones:
     Writes all regions to cellZones preserving Fluent names
-    Writes all region internal face to faceZones preserving Fluent names 
+    Writes all region internal face to faceZones preserving Fluent names
     lines: 1545 - 1667
     Zones are usefull for porous media and MRF calculations
 
 2. Zone Access
-    - Zones are simple lists of label lists that can be accessed from polyMesh 
+    - Zones are simple lists of label lists that can be accessed from polyMesh
     with the cellZones(), faceZones() and pointZones() member functions
 
     - Example (Members from polyMesh.H and ZoneMesh.H):
     label thisCellZoneID = mesh.cellZones().findZoneID("thisZoneName");
     const labelList& thisCellZone = mesh.cellZones()[thisCellZoneID];
-    
+
     - Zone integrity is preserved during mesh modification and decompomposition.
-    
-    - Once created via addZones, zones allow modification through non-const 
+
+    - Once created via addZones, zones allow modification through non-const
     access
 
 3. Fluent boundary types.
     - All internal and baffle elements are ignored during conversion
-    
+
     - Boundary faces labelled as internal (i.e. interior, interface, internal,
     solid, fan, radiator, porous-jump) but that are in fact external boundaries
     will be added to a default wall boundary.

applications/utilities/mesh/conversion/gmshToFoam/Make/options

@@ -4,4 +4,4 @@ EXE_INC = \
 
 EXE_LIBS = \
     -ldynamicMesh \
-    -lmeshTools 
+    -lmeshTools

applications/utilities/mesh/conversion/sammToFoam/readBoundary.C

@@ -188,7 +188,7 @@ void sammMesh::readBoundary()
                 boundary_[patchLabel][nPatchFaces[patchLabel]]
                     = face(pointLabels);
             }
-            
+
             // increment counter of faces in current patch
             nPatchFaces[patchLabel]++;
         }

applications/utilities/mesh/conversion/starToFoam/readBoundary.C

@@ -186,7 +186,7 @@ void starMesh::readBoundary()
                 boundary_[patchLabel][nPatchFaces[patchLabel]]
                     = face(pointLabels);
             }
-            
+
             // increment counter of faces in current patch
             nPatchFaces[patchLabel]++;
         }

applications/utilities/mesh/conversion/starToFoam/readPoints.C

@@ -95,7 +95,7 @@ void starMesh::readPoints(const scalar scaleFactor)
                 pointLabel = readVtxLabel(pointsFile);
 
                 if (!pointsFile) break;
-                
+
                 maxLabel = max(maxLabel, pointLabel);
 
                 x = readVtxCmpt(pointsFile);

applications/utilities/mesh/conversion/writeMeshObj/Make/options

@@ -2,4 +2,4 @@ EXE_INC = \
     -I$(LIB_SRC)/meshTools/lnInclude
 
 EXE_LIBS = \
-    -lmeshTools 
+    -lmeshTools

applications/utilities/mesh/generation/extrudeMesh/extrudedMesh/extrudedMesh.C

@@ -95,7 +95,7 @@ Foam::Xfer<Foam::faceList> Foam::extrudedMesh::extrudedFaces
 
     const label nLayers = model.nLayers();
 
-    label nFaces = 
+    label nFaces =
         (nLayers + 1)*surfaceFaces.size() + nLayers*surfaceEdges.size();
 
     faceList eFaces(nFaces);

applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C

@@ -192,7 +192,7 @@ int main(int argc, char *argv[])
         autoRefineDriver::debug = debug;
         autoSnapDriver::debug = debug;
         autoLayerDriver::debug = debug;
-    }        
+    }
 
 
     // Read geometry

applications/utilities/mesh/manipulation/checkMesh/printMeshStats.C

@@ -129,7 +129,7 @@ void Foam::printMeshStats(const polyMesh& mesh, const bool allTopology)
     }
 
     reduce(nHex,sumOp<label>());
-    reduce(nPrism,sumOp<label>()); 
+    reduce(nPrism,sumOp<label>());
     reduce(nWedge,sumOp<label>());
     reduce(nPyr,sumOp<label>());
     reduce(nTetWedge,sumOp<label>());

applications/utilities/mesh/manipulation/moveEngineMesh/moveEngineMesh.C

@@ -26,7 +26,7 @@ Application
     moveEngineMesh
 
 Description
-    Solver for moving meshes for engine calculations. 
+    Solver for moving meshes for engine calculations.
 
 \*---------------------------------------------------------------------------*/
 
@@ -54,7 +54,7 @@ int main(int argc, char *argv[])
         mesh.move();
 
         runTime.write();
- 
+
         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
             << nl << endl;

applications/utilities/mesh/manipulation/objToVTK/Make/options

-/* 
+/*
 EXE_INC = \
     -I$(LIB_SRC)/triSurface/lnInclude
 

applications/utilities/mesh/manipulation/setSet/writePatch.C

@@ -67,7 +67,7 @@ void writePatch
     //------------------------------------------------------------------
     //