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