& freeSurfFac, implicSurfPress, implicDiv2Dflow, implicitNHPress,
_RL implicitNHPress
& .AND. implicitNHPress.EQ.1. _d 0 ) THEN
IF ( implicitNHPress*implicSurfPress*implicDiv2Dflow.NE.1.
& implicitNHPress.NE.implicSurfPress ) THEN
CALL WRITE_0D_RL( implicitNHPress, INDEX_NONE,
nhFac = pfFacMom*implicitNHPress
& *implicitNHPress*implicDiv2DFlow
& implicSurfPress, implicDiv2Dflow, implicitNHPress,
implicitNHPress = UNSET_RL
IF ( implicitNHPress.EQ.UNSET_RL )
& implicitNHPress = implicSurfPress
& implicitNHPress.LT.1. _d 0 ) 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)