& usingPCoords, usingZCoords, useDynP_inEos_Zc,
LOGICAL usingPCoords
IF ( usingPCoords ) THEN
IF ( usingPCoords .AND. k.NE.1 .AND.
IF ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
ELSEIF ( usingPCoords ) THEN
& useRealFreshWaterFlux .AND. usingPCoords ) THEN
IF ( usingPCoords ) k = Nr
& (usingPCoords .AND. K .EQ. Nr) ) THEN
IF ( usingPCoords ) rUnits = 'Pa'
CALL WRITE_0D_L( usingPCoords, INDEX_NONE,
IF ( ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
IF ( usingPCoords .AND. bathyFile.NE.' ' ) THEN
IF (usingPCoords) THEN
IF (usingPCoords .AND. bathyFile .NE. ' '
IF (usingPCoords .OR. bathyFile .EQ. ' ') THEN
IF ( equationOfState .EQ. 'JMD95Z' .AND. usingPCoords ) THEN
IF ( pLoadFile .NE. ' ' .AND. usingPCoords ) THEN
usingPCoords = .FALSE.
usingPCoords = .TRUE.
IF ( usingPCoords ) rEmpty = rF(Nr+1)
IF ( usingPCoords ) rSigBndRS = rF(Nr+1)
IF ( usingPCoords .AND. setInterFDr ) THEN
IF ( usingPCoords .AND. setCenterDr ) THEN
IF ( k.EQ.Nr .AND. myTime.NE.baseTime .AND. usingPCoords
IF ( usingPCoords ) ks = Nr
if ( (kSize.EQ.1).AND.(usingPCoords) ) then
& ( .NOT.nonHydrostatic .OR. usingPCoords
IF ( usingPCoords ) diagUnits = 's^2/m^2 '
IF ( usingPCoords ) downward = -1