Clean-up and formatting

applications/solvers/coupled/conjugateHeatFoam/solveEnergy.H

@@ -20,9 +20,6 @@
           + 3.0*radiation->Rp()*pow4(T)
         );
 
-        //Done automatically in 1.6.1-ext
-        //TFluidEqn->boundaryManipulate(T.boundaryField());
-
         fvScalarMatrix* TSolidEqn = new fvScalarMatrix
         (
             fvm::ddt(rhoCpsolid, Tsolid)
@@ -30,10 +27,6 @@
           + fvm::SuSp(-solidThermo.S()/Tsolid, Tsolid)
         );
 
-        //Done automatically in 1.6.1-ext
-        //TSolidEqn->boundaryManipulate(Tsolid.boundaryField());
-
-
         // Add fluid equation
         TEqns.set(0, TFluidEqn);
 

applications/solvers/coupled/pUCoupledFoam/createFields.H

@@ -29,6 +29,7 @@ volVectorField U
 #   include "createPhi.H"
 
     singlePhaseTransportModel laminarTransport(U, phi);
+
     autoPtr<incompressible::RASModel> turbulence
     (
         incompressible::RASModel::New(U, phi, laminarTransport)
@@ -67,3 +68,4 @@ volScalarField rAU
     );
 
     mesh.schemesDict().setFluxRequired(p.name());
+    mesh.schemesDict().setFluxRequired("pcorr");

applications/solvers/incompressible/simpleFoam/pEqn.H

@@ -3,7 +3,7 @@
 
     // Prepare clean 1/Ap without contribution from under-relaxation
     // HJ, 26/Oct/2015
-    volScalarField rUA
+    volScalarField rAU
     (
         "(1|A(U))",
         1/HUEqn().A()
@@ -13,7 +13,7 @@
     // the flux precursor
     U.storePrevIter();
 
-    U = rUA*HUEqn().H();
+    U = rAU*HUEqn().H();
     HUEqn.clear();
 
     phi = fvc::interpolate(U) & mesh.Sf();