& + implicSurfPress*implicDiv2DFlow
& *implicitNHPress*implicDiv2DFlow
& +implicDiv2DFlow*deltaTfreesurf
& * implicSurfPress*implicDiv2DFlow
wFacKp = implicDiv2DFlow*deepFac2F(kp1)*rhoFacF(kp1)
wFacKm = implicDiv2DFlow*deepFac2F( k )*rhoFacF( k )
& tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow
& *implicSurfPress*implicDiv2DFlow