IOobject alphaPhi10Header(IOobject::groupName("alphaPhi0", alpha1.group()), runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE)
surfaceScalarField alphaPhiUn(IOobject("alphaPhiUn", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedScalar(phi.dimensions(), Zero))