& freeSurfFac, implicSurfPress, implicDiv2DFlow, implicitNHPress,
_RL implicDiv2DFlow
& .AND. implicDiv2DFlow.EQ.0. _d 0
& .AND. implicDiv2DFlow.NE.1. ) THEN
IF ( implicitNHPress*implicSurfPress*implicDiv2DFlow.NE.1.
IF ( implicitNHPress*implicSurfPress*implicDiv2DFlow.EQ.0.
CALL WRITE_0D_RL( implicDiv2DFlow, INDEX_NONE,
IF (implicDiv2DFlow.LT.1.) THEN
& + implicSurfPress*implicDiv2DFlow
& *implicitNHPress*implicDiv2DFlow
& *implicitNHPress*implicDiv2DFlow )
& implicSurfPress, implicDiv2DFlow, implicitNHPress,
& *implicDiv2DFlow*deltaTFreeSurf
tmpFac = freeSurfFac*implicDiv2DFlow/deltaTMom
& * implicSurfPress*implicDiv2DFlow
wFacKp = implicDiv2DFlow*deepFac2F(kp1)*rhoFacF(kp1)
wFacKm = implicDiv2DFlow*deepFac2F( k )*rhoFacF( k )
IF ( oldFreeSurfTerm .OR. implicDiv2DFlow.EQ.zeroRL ) THEN
& / ( implicDiv2DFlow*deltaTFreeSurf )
tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow/deltaTMom
implicDiv2DFlow = 1. _d 0
tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow
& *implicSurfPress*implicDiv2DFlow