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