addSurfPhiAnom = select_rStar.EQ.0 .AND. nonlinFreeSurf.GT.3
if(usefizhi .and. (select_rStar.GT.0) )then
IF ( select_rStar.GT.0 ) THEN
IF ( select_rStar.NE.0 )
IF ( select_rStar.NE.0 ) THEN
CALL WRITE_0D_I( select_rStar, INDEX_NONE,
IF (select_rStar .NE. 0) THEN
IF (select_rStar.NE.0 .AND. .NOT.exactConserv) THEN
ELSEIF (select_rStar .NE. 0) THEN
IF (select_rStar.GT.0) THEN
& select_rStar,
& nonlinFreeSurf, select_rStar,
INTEGER select_rStar
if(select_rStar.GT.0)then
IF ( momAdvection .AND. select_rStar.GT.0 ) THEN
IF ( momAdvection .AND. select_rStar.LT.0 ) THEN
IF ( select_rStar.EQ.0 .AND. .NOT.rigidLid ) 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
select_rStar = 0
IF (select_rStar.NE.0) THEN