| 252 |
C-- endif buoyancyRelation |
C-- endif buoyancyRelation |
| 253 |
ENDIF |
ENDIF |
| 254 |
|
|
| 255 |
|
C-- fill-in phiRef array (presently not used) |
| 256 |
|
IF ( buoyancyRelation.EQ.'OCEANIC' ) THEN |
| 257 |
|
c phiRef(1) = gravitySign*gravity*(rF(1)-Ro_SeaLevel) |
| 258 |
|
phiRef(1) = 0. _d 0 |
| 259 |
|
DO k=1,Nr |
| 260 |
|
phiRef(2*k) = gravitySign*gravity*(rC(k) - Ro_SeaLevel) |
| 261 |
|
phiRef(2*k+1) = gravitySign*gravity*(rF(k+1)-Ro_SeaLevel) |
| 262 |
|
ENDDO |
| 263 |
|
ELSEIF ( buoyancyRelation.EQ.'OCEANICP' ) THEN |
| 264 |
|
C should be : phiRef = phi_Origin - (rC - Ro_SeaLevel)/rhoConst |
| 265 |
|
C- but since the reference geopotential "phi_Origin" @ p = rF(k=1) |
| 266 |
|
C is not currently stored, we only get a relative geopotential: |
| 267 |
|
phiRef(1) = -recip_rhoConst*rF(1) |
| 268 |
|
DO k=1,Nr |
| 269 |
|
phiRef(2*k) = -recip_rhoConst*rC(k) |
| 270 |
|
phiRef(2*k+1) = -recip_rhoConst*rF(k+1) |
| 271 |
|
ENDDO |
| 272 |
|
ENDIF |
| 273 |
|
|
| 274 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
| 275 |
|
|
| 276 |
IF ( usingZCoords .AND. rhoRefFile .NE. ' ' ) THEN |
IF ( usingZCoords .AND. rhoRefFile .NE. ' ' ) THEN |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch