& selectSigmaCoord,
INTEGER selectSigmaCoord
ELSEIF ( selectSigmaCoord.NE.0 ) THEN
IF ( selectSigmaCoord.EQ.0 ) THEN
useFVgradPhi = selectSigmaCoord.NE.0
IF ( select_rStar.GE.1 .OR. selectSigmaCoord.GE.1 ) THEN
IF ( selectSigmaCoord.NE.0 ) THEN
& ( usingPCoords .OR. deepAtmosphere .OR. selectSigmaCoord.NE.0
CALL WRITE_0D_I( selectSigmaCoord, INDEX_NONE,
ELSEIF ( nonlinFreeSurf.GT.0 .AND. selectSigmaCoord.EQ.0 ) THEN
& selectSigmaCoord, rSigmaBnd, hybSigmFile,
IF ( selectSigmaCoord .EQ. 0 ) THEN
IF ( nonlinFreeSurf.GT.0 .AND. selectSigmaCoord.NE.0 ) THEN
selectSigmaCoord = 0