& use3Dsolver, implicitIntGravWave, staggerTimeStep,
LOGICAL staggerTimeStep
IF ( .NOT.(staggerTimeStep.OR.implicitIntGravWave) ) THEN
IF ( staggerTimeStep .AND. LS_whenToSample.LT.2 ) THEN
IF ( selectAddFluid.GE.1 .AND. .NOT.staggerTimeStep ) THEN
IF ( implicitIntGravWave .AND. staggerTimeStep ) THEN
CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,
IF ( nonlinFreeSurf.GT.0 .AND. staggerTimeStep ) THEN
IF ( .NOT.staggerTimeStep ) THEN
IF ( staggerTimeStep ) THEN
IF (staggerTimeStep) iterNb = myIter - 1
IF (staggerTimeStep.AND.nonlinFreeSurf.GT.0) THEN
IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
& implicitIntGravWave, staggerTimeStep, doResetHFactors,
IF ( addFrictionHeating .AND. staggerTimeStep ) THEN
IF ( addFrictionHeating .AND. .NOT.staggerTimeStep ) THEN
staggerTimeStep = .FALSE.
IF ( .NOT.staggerTimeStep ) wrTime = myTime + deltaTClock
IF ( .NOT.staggerTimeStep .AND. .NOT. implicitIntGravWave ) THEN
IF ( staggerTimeStep.AND.select_rStar.GT.0 ) THEN
IF (staggerTimeStep) THEN