if ( equationOfState .eq. 'JMD95Z' .and. usingPCoords ) then
IF (usingPCoords) K = Nr
IF (usingPCoords) THEN
IF ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
IF ( usingPCoords ) THEN
CALL WRITE_0D_L( usingPCoords, INDEX_NONE,
IF ( usingPCoords .AND. bathyFile.NE.' ' ) THEN
IF (usingPCoords .AND. bathyFile .NE. ' '
IF (usingPCoords .OR. bathyFile .EQ. ' ') THEN
ELSEIF ( usingPCoords ) THEN
usingPCoords = .FALSE.
usingPCoords = .TRUE.
& usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
LOGICAL usingPCoords
IF ( usingPCoords ) km = Nr
IF ( usingPCoords ) diagUnits = 's^2/m^2 '
& useRealFreshWaterFlux .AND. usingPCoords ) THEN