dimensionedScalar nu
(
"nu",
dimViscosity,
transportProperties
);
wordList patchTypes(mesh.boundaryMesh().size(),"fixedValue");
volScalarField nu0
(
IOobject
(
"nu0",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
nu,
patchTypes
);
const scalar nuWanted(1000.0);
labelUList neiCells(nu0.boundaryFieldRef()[0].patch().faceCells());
forAll(neiCells, cellI)
{
nu0.primitiveFieldRef()[neiCells[cellI]]=nuWanted;
}