& freeSurfFac, implicSurfPress, implicDiv2Dflow, implicitNHPress,
_RL implicSurfPress
IF ( implicitNHPress*implicSurfPress*implicDiv2Dflow.NE.1.
& implicitNHPress.NE.implicSurfPress ) THEN
CALL WRITE_0D_RL( implicSurfPress, INDEX_NONE,
psFac = pfFacMom*implicSurfPress
IF (implicSurfPress.NE.1.) THEN
& + implicSurfPress*implicDiv2DFlow
& implicSurfPress, implicDiv2Dflow, implicitNHPress,
& implicitNHPress = implicSurfPress
& * implicSurfPress*implicDiv2DFlow
implicSurfPress = 1. _d 0
psFac = pfFacMom*(1. _d 0 - implicSurfPress)
& *implicSurfPress*implicDiv2DFlow