IF ( usingPCoords ) ks = Nr
& usingPCoords, usingZCoords,
LOGICAL usingPCoords
IF ( usingPCoords ) THEN
ELSEIF ( usingPCoords ) THEN
IF ( (nonlinFreeSurf.GT.0 .OR. usingPCoords)
& useRealFreshWaterFlux .AND. usingPCoords ) THEN
IF ( usingPCoords .AND. gravityFile .NE. ' ' ) THEN
& ( usingPCoords .OR. deepAtmosphere .OR. selectSigmaCoord.NE.0
IF ( usingPCoords ) k = Nr
CALL WRITE_0D_L( usingPCoords, INDEX_NONE,
IF ( usingPCoords ) rUnits = 'Pa'
& (usingPCoords .AND. K .EQ. Nr) ) THEN
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
IF ( usingPCoords .AND. .NOT.uniformLin_PhiSurf ) THEN
IF ( usingPCoords .AND. top_Pres.NE.UNSET_RL ) 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. myIter.NE.0 .AND. usingPCoords
IF ( oldFreeSurfTerm .OR. usingPCoords ) THEN
if ( (kSize.EQ.1).AND.(usingPCoords) ) then
& ( .NOT.nonHydrostatic .OR. usingPCoords
IF ( usingPCoords ) diagUnits = 's^2/m^2 '
IF ( usingPCoords ) downward = -1