[bye bye my blue]

in yPlus.C,

forAll(patches, patchi)
{
const fvPatch& patch = patches[patchi];

if (isA<nutWallFunctionFvPatchScalarField>(nutBf[patchi]))
{
const nutWallFunctionFvPatchScalarField& nutPf =
dynamic_cast<const nutWallFunctionFvPatchScalarField&>
(
nutBf[patchi]
);

yPlusBf[patchi] = nutPf.yPlus();
}
else if (isA<wallFvPatch>(patch))
{
yPlusBf[patchi] =
d[patchi]
*sqrt
(
nuEffBf[patchi]
*mag(turbModel.U().boundaryField()[patchi].snGrad())
)/nuBf[patchi];
}
}
nuEffBf[patchi] means eddy viscosity?
and mag(turbModel.U().boundaryField()[patchi].snGrad()) = abs(dU/dy). is it right??

I saw that openfoam yPlus is actually yStar. But I think openfoam yPlus is real yPlus in openfoam4.0