Info<< "Reading field U\n" << endl;
volVectorField U
(
    IOobject
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    ),
    mesh
);

#include <finiteVolume/createPhi.H>

singlePhaseTransportModel laminarTransport(U, phi);

autoPtr<incompressible::LESModel> sgsModel
(
    incompressible::LESModel::New(U, phi, laminarTransport)
);

volScalarField::GeometricBoundaryField d = nearWallDist(mesh).y();

// ************************ vim: set sw=4 sts=4 et: ************************ //
