& freeSurfFac, implicSurfPress, implicDiv2DFlow, implicitNHPress,
_RL implicitNHPress
& .AND. implicitNHPress.EQ.1. _d 0 ) THEN
IF ( implicitNHPress*implicSurfPress*implicDiv2DFlow.NE.1.
IF ( implicitNHPress*implicSurfPress*implicDiv2DFlow.EQ.0.
& implicitNHPress.NE.implicSurfPress ) THEN
CALL WRITE_0D_RL( implicitNHPress, INDEX_NONE,
nhFac = pfFacMom*implicitNHPress
& *implicitNHPress*implicDiv2DFlow
& *implicitNHPress*implicDiv2DFlow )
& implicSurfPress, implicDiv2DFlow, implicitNHPress,
implicitNHPress = UNSET_RL
IF ( implicitNHPress.EQ.UNSET_RL )
& implicitNHPress = implicSurfPress
& ( implicitNHPress.LT.oneRL .OR. selectP_inEOS_Zc.EQ.3 ) ) THEN
IF ( use3Dsolver .AND. implicitNHPress.NE.1. _d 0 ) THEN
nhFac = pfFacMom*(1. _d 0 - implicitNHPress)
IF ( implicitNHPress.NE.1. _d 0 ) THEN
tmpFac = pfFacMom*(1. _d 0 - implicitNHPress)