{
    p.boundaryField().updateCoeffs();

    volScalarField A
    (
  	IOobject
  	(
    	 "A",
    	 runTime.timeName(),
    	 mesh,
    	 IOobject::NO_READ,
    	 IOobject::NO_WRITE
  	),
	UrelEqn().A()
    );
    A.write();

    volVectorField H
    (
  	IOobject
  	(
    	 "H",
    	 runTime.timeName(),
    	 mesh,
    	 IOobject::NO_READ,
    	 IOobject::NO_WRITE
  	),
	UrelEqn().H()
    );
    H.write();

    volScalarField H1
    (
  	IOobject
  	(
    	 "H1",
    	 runTime.timeName(),
    	 mesh,
    	 IOobject::NO_READ,
    	 IOobject::NO_WRITE
  	),
	UrelEqn().H1()
    );
    H1.write();

    volScalarField rAUrel(1.0/UrelEqn().A());
    Urel = rAUrel*UrelEqn().H();
    UrelEqn.clear();

    phi = fvc::interpolate(Urel, "interpolate(HbyA)") & mesh.Sf();
    adjustPhi(phi, Urel, p);

    // Non-orthogonal pressure corrector loop
    while (simple.correctNonOrthogonal())
    {
        fvScalarMatrix pEqn
        (
            fvm::laplacian(rAUrel, p) == fvc::div(phi)
        );

        pEqn.setReference(pRefCell, pRefValue);

        pEqn.solve();

        if (simple.finalNonOrthogonalIter())
        {
            phi -= pEqn.flux();
        }
    }

    #include "continuityErrs.H"

    // Explicitly relax pressure for momentum corrector
    p.relax();

    // Momentum corrector
    Urel -= rAUrel*fvc::grad(p);
    Urel.correctBoundaryConditions();
    sources.correct(Urel);
}

    // Relative momentum predictor

    volVectorField divPhiUrel
    (
  	IOobject
  	(
    	 "divPhiUrel",
    	 runTime.timeName(),
    	 mesh,
    	 IOobject::NO_READ,
    	 IOobject::NO_WRITE
  	),
	fvc::div(phi, Urel)
    );
    divPhiUrel.write();

    volVectorField gradP
    (
  	IOobject
  	(
    	 "gradP",
    	 runTime.timeName(),
    	 mesh,
    	 IOobject::NO_READ,
    	 IOobject::NO_WRITE
  	),
	fvc::grad(p)
    );
    gradP.write();

    tmp<fvVectorMatrix> UrelEqn
    (
        fvm::div(phi, Urel)
      + turbulence->divDevReff(Urel)
      + SRF->Su()
     ==
        sources(Urel)
    );

    UrelEqn().relax();

    sources.constrain(UrelEqn());

    solve(UrelEqn() == -fvc::grad(p));

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-92812fe8fd9e
Exec   : SRFSimpleFoamTest -parallel
Date   : Feb 13 2013
Time   : 14:49:31
Host   : "x470172"
PID    : 7162
Case   : /c/home/dh/parallel_test/10_simple
nProcs : 4
Slaves : 
3
(
"x470172.7163"
"x470172.7164"
"x470172.7165"
)

Pstream initialized with:
    floatTransfer     : 0
    nProcsSimpleSum   : 0
    commsType         : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 6000

// using new solver syntax:
p
{
    solver          GAMG;
    tolerance       1e-16;
    relTol          0.01;
    smoother        GaussSeidel;
    cacheAgglomeration true;
    nCellsInCoarsestLevel 10;
    agglomerator    faceAreaPair;
    mergeLevels     1;
    minIter         0;
    maxIter         1000;
    nPreSweeps      0;
    nPostSweeps     2;
    nFinestSweeps   2;
    scaleCorrection true;
    directSolveCoarsest false;
}

// using new solver syntax:
Urel
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
k
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
epsilon
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
omega
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

Reading field p

Reading field Urel

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
Creating SRF model

Selecting SRFModel rpm
No field sources present


SIMPLE: no convergence criteria found. Calculations will run for 1e+15 steps.


Starting time loop

phi: phi
Compressible: 0
Turbulent: 0
LES: 0
--> FOAM Warning : 
    From function probes::read()
    in file patch/patchFieldFunctionObject/patchFieldFunctionObject.C at line 91
    Unknown field ptot when reading dictionary "::PsiTotal"
    Can only probe registered volScalar, volVector, volSphericalTensor, volSymmTensor and volTensor fields

Time = 6001

smoothSolver:  Solving for Urelx, Initial residual = 0.00191058, Final residual = 1.30178e-05, No Iterations 5
smoothSolver:  Solving for Urely, Initial residual = 0.000450372, Final residual = 3.12678e-06, No Iterations 5
smoothSolver:  Solving for Urelz, Initial residual = 0.00112577, Final residual = 6.27861e-06, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0844304, Final residual = 0.000679339, No Iterations 3
time step continuity errors : sum local = 0.543478, global = -0.0155972, cumulative = -0.0155972
smoothSolver:  Solving for omega, Initial residual = 0.000577573, Final residual = 3.99261e-06, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 0.000372194, Final residual = 1.60936e-06, No Iterations 5
ExecutionTime = 5.64 s  ClockTime = 6 s

forces output:
    forces(pressure, viscous)((-76.5138 -26.4733 213.565) (0.338827 0.0744777 0.108567))
    moment(pressure, viscous)((12.245 -18.9689 -1.86084) (0.00468325 0.013836 -0.0316924))

fieldMinMax output:
    min(p) = 957.104
    max(p) = 2348.45

fieldMinMax output:
    min(mag(Urel)) = 0
    max(mag(Urel)) = 19.4621

End

Finalising parallel run

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-92812fe8fd9e
Exec   : SRFSimpleFoamTest
Date   : Feb 13 2013
Time   : 14:50:33
Host   : "x470172"
PID    : 7189
Case   : /c/home/dh/parallel_test/10_simple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 6000

// using new solver syntax:
p
{
    solver          GAMG;
    tolerance       1e-16;
    relTol          0.01;
    smoother        GaussSeidel;
    cacheAgglomeration true;
    nCellsInCoarsestLevel 10;
    agglomerator    faceAreaPair;
    mergeLevels     1;
    minIter         0;
    maxIter         1000;
    nPreSweeps      0;
    nPostSweeps     2;
    nFinestSweeps   2;
    scaleCorrection true;
    directSolveCoarsest false;
}

// using new solver syntax:
Urel
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
k
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
epsilon
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

// using new solver syntax:
omega
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    nSweeps         1;
    tolerance       1e-16;
    relTol          0.01;
    maxIter         500;
}

Reading field p

Reading field Urel

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
Creating SRF model

Selecting SRFModel rpm
No field sources present


SIMPLE: no convergence criteria found. Calculations will run for 1e+15 steps.


Starting time loop

phi: phi
Compressible: 0
Turbulent: 0
LES: 0
--> FOAM Warning : 
    From function probes::read()
    in file patch/patchFieldFunctionObject/patchFieldFunctionObject.C at line 91
    Unknown field ptot when reading dictionary "::PsiTotal"
    Can only probe registered volScalar, volVector, volSphericalTensor, volSymmTensor and volTensor fields

Time = 6001

smoothSolver:  Solving for Urelx, Initial residual = 0.00191839, Final residual = 1.27509e-05, No Iterations 5
smoothSolver:  Solving for Urely, Initial residual = 0.000452133, Final residual = 3.0692e-06, No Iterations 5
smoothSolver:  Solving for Urelz, Initial residual = 0.00114632, Final residual = 6.2392e-06, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0912255, Final residual = 0.000644178, No Iterations 4
time step continuity errors : sum local = 0.518156, global = 0.032982, cumulative = 0.032982
smoothSolver:  Solving for omega, Initial residual = 0.000587208, Final residual = 4.03108e-06, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 0.000380995, Final residual = 1.63151e-06, No Iterations 5
ExecutionTime = 18.25 s  ClockTime = 20 s

forces output:
    forces(pressure, viscous)((-76.7439 -26.8078 213.407) (0.338762 0.0744753 0.108581))
    moment(pressure, viscous)((12.2519 -18.9781 -1.87774) (0.00468457 0.0138324 -0.0316869))

fieldMinMax output:
    min(p) = 956.365
    max(p) = 2348.45

fieldMinMax output:
    min(mag(Urel)) = 0
    max(mag(Urel)) = 19.4621

End