From b1c3e76352e1a0c65e2fbb82b0d22491cf3ca57d Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Thu, 10 Oct 2013 10:50:28 +0100
Subject: [PATCH] Formatting and stability changes

---
 applications/solvers/combustion/XiFoam/XiFoam.C       | 8 ++++----
 applications/solvers/combustion/XiFoam/bEqn.H         | 9 +++++----
 applications/solvers/combustion/XiFoam/createFields.H | 1 -
 3 files changed, 9 insertions(+), 9 deletions(-)

diff --git a/applications/solvers/combustion/XiFoam/XiFoam.C b/applications/solvers/combustion/XiFoam/XiFoam.C
index b5c2b2e61..25b74f745 100644
--- a/applications/solvers/combustion/XiFoam/XiFoam.C
+++ b/applications/solvers/combustion/XiFoam/XiFoam.C
@@ -90,8 +90,8 @@ int main(int argc, char *argv[])
         // --- Pressure-velocity PIMPLE corrector loop
         for (int oCorr=0; oCorr<nOuterCorr; oCorr++)
         {
-#       include "rhoEqn.H"
-#       include "UEqn.H"
+#           include "rhoEqn.H"
+#           include "UEqn.H"
 
 #           include "ftEqn.H"
 #           include "bEqn.H"
@@ -106,8 +106,8 @@ int main(int argc, char *argv[])
             // --- PISO loop
             for (int corr=1; corr<=nCorr; corr++)
             {
-#           include "pEqn.H"
-        }
+#              include "pEqn.H"
+            }
 
             turbulence->correct();
         }
diff --git a/applications/solvers/combustion/XiFoam/bEqn.H b/applications/solvers/combustion/XiFoam/bEqn.H
index 74ba01501..715b87310 100644
--- a/applications/solvers/combustion/XiFoam/bEqn.H
+++ b/applications/solvers/combustion/XiFoam/bEqn.H
@@ -51,7 +51,7 @@ if (ign.ignited())
         fvm::ddt(rho, b)
       + mvConvection->fvmDiv(phi, b)
       + fvm::div(phiSt, b, "div(phiSt,b)")
-      - fvm::Sp(fvc::div(phiSt), b)
+      + fvm::SuSp(-fvc::div(phiSt), b)
       - fvm::laplacian(turbulence->alphaEff(), b)
     );
 
@@ -65,7 +65,8 @@ if (ign.ignited())
     // ~~~~~~~~~~~
     bEqn.solve();
 
-    Info<< "min(b) = " << min(b).value() << endl;
+    Info<< "min(b) = " << min(b).value()
+        << " max(b) = " << max(b).value() << endl;
 
 
     // Calculate coefficients for Gulder's flame speed correlation
@@ -138,7 +139,7 @@ if (ign.ignited())
         (
             fvm::ddt(rho, Su)
           + fvm::div(phi + phiXi, Su, "div(phiXi,Su)")
-          - fvm::Sp(fvc::div(phiXi), Su)
+          + fvm::SuSp(-fvc::div(phiXi), Su)
           ==
           - fvm::SuSp(-rho*Rc*Su0/Su, Su)
           - fvm::SuSp(rho*(sigmas + Rc), Su)
@@ -203,7 +204,7 @@ if (ign.ignited())
         (
             fvm::ddt(rho, Xi)
           + fvm::div(phi + phiXi, Xi, "div(phiXi,Xi)")
-          - fvm::Sp(fvc::div(phiXi), Xi)
+          + fvm::SuSp(-fvc::div(phiXi), Xi)
          ==
             rho*R
           - fvm::Sp(rho*(R - G), Xi)
diff --git a/applications/solvers/combustion/XiFoam/createFields.H b/applications/solvers/combustion/XiFoam/createFields.H
index ef16bd615..93aa61ea0 100644
--- a/applications/solvers/combustion/XiFoam/createFields.H
+++ b/applications/solvers/combustion/XiFoam/createFields.H
@@ -30,7 +30,6 @@
 
     const volScalarField& T = thermo.T();
 
-
     Info<< "\nReading field U\n" << endl;
     volVectorField U
     (