CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
& buoyancyRelation, eosType,
IF ( buoyancyRelation .NE. 'OCEANIC' ) THEN
IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN
ELSEIF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
ELSEIF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN
& .AND. buoyancyRelation.EQ.'OCEANICP' ) THEN
CALL WRITE_0D_C( buoyancyRelation, -1, INDEX_NONE,
IF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN
ELSEIF ( buoyancyRelation.EQ.'OCEANICP' ) THEN
& rhoConst, thetaConst, rhoConstFresh, buoyancyRelation,
IF ( buoyancyRelation.EQ.'ATMOSPHERIC' ) THEN
ELSEIF ( buoyancyRelation.EQ.'OCEANICP') THEN
ELSEIF ( buoyancyRelation.EQ.'OCEANIC' ) THEN
& buoyancyRelation(1:11)
buoyancyRelation = 'OCEANIC'
& buoyancyRelation.EQ.'OCEANICP' ) THEN
ELSEIF ( buoyancyRelation.EQ.'ATMOSPHERIC' ) THEN
IF ( buoyancyRelation.EQ.'OCEANIC' ) THEN