IF ( useOBCS .AND. nonlinFreeSurf.GT.0 )
IF (nonlinFreeSurf.GT.0) THEN
addSurfPhiAnom = select_rStar.EQ.0 .AND. nonlinFreeSurf.GT.3
IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
IF ( nonlinFreeSurf.GT.2) THEN
ELSEIF ( nonlinFreeSurf.GT.0) THEN
IF ( nonlinFreeSurf.GT.0) THEN
IF ( momStepping .AND. nonlinFreeSurf.GT.2 ) THEN
IF ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
CALL WRITE_0D_I( nonlinFreeSurf, INDEX_NONE,
IF (useRealFreshWaterFlux .AND. nonlinFreeSurf.GT.0) THEN
IF (nonlinFreeSurf.NE.0) THEN
IF (nonlinFreeSurf.NE.0 .AND. .NOT.exactConserv) THEN
IF (nonlinFreeSurf.NE.0 .AND. nonHydrostatic) THEN
IF ( nonlinFreeSurf.NE.0 .AND.
& .AND. nonlinFreeSurf.GT.1) THEN
IF ( nonlinFreeSurf.GE.0)
IF (staggerTimeStep.AND.nonlinFreeSurf.GT.0) THEN
& exactConserv,uniformLin_PhiSurf,nonlinFreeSurf,hFacInf,hFacSup,
& nonlinFreeSurf, select_rStar,
INTEGER nonlinFreeSurf
IF ( nonlinFreeSurf.GT.0 .AND.
IF (nonlinFreeSurf .GE. 0) THEN
IF (nonlinFreeSurf .GT. 0) THEN
IF (select_rStar.GE.2 .AND. nonlinFreeSurf.GE.1 ) THEN
IF (select_rStar.GE.2 .AND. nonlinFreeSurf.GE.4 ) THEN
ELSEIF (select_rStar.GE.1 .AND. nonlinFreeSurf.GE.4 ) THEN
nonlinFreeSurf = 0