65 |
C phiHyd - Hydrostatic part of the potential phi. |
C phiHyd - Hydrostatic part of the potential phi. |
66 |
C In z coords phiHyd is the hydrostatic pressure anomaly |
C In z coords phiHyd is the hydrostatic pressure anomaly |
67 |
C In p coords phiHyd is the geopotential surface height anomaly. |
C In p coords phiHyd is the geopotential surface height anomaly. |
68 |
|
C etaSurfX, etaSurfY - Holds surface elevation gradient in X and Y. |
69 |
C iMin, iMax - Ranges and sub-block indices on which calculations |
C iMin, iMax - Ranges and sub-block indices on which calculations |
70 |
C jMin, jMax are applied. |
C jMin, jMax are applied. |
71 |
C bi, bj |
C bi, bj |
98 |
_RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL buoyKM1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
99 |
_RL buoyK (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL buoyK (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
100 |
_RL rhotmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rhotmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
101 |
_RL pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL etaSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
102 |
_RL pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL etaSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
103 |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K13 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
104 |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K23 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
105 |
_RL K33 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL K33 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
230 |
C-- Calculate gradient of surface pressure |
C-- Calculate gradient of surface pressure |
231 |
CALL CALC_GRAD_ETA_SURF( |
CALL CALC_GRAD_ETA_SURF( |
232 |
I bi,bj,iMin,iMax,jMin,jMax, |
I bi,bj,iMin,iMax,jMin,jMax, |
233 |
O pSurfX,pSurfY, |
O etaSurfX,etaSurfY, |
234 |
I myThid) |
I myThid) |
235 |
|
|
236 |
C-- Update fields in top level according to tendency terms |
C-- Update fields in top level according to tendency terms |
237 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
238 |
I bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myTime,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,K,etaSurfX,etaSurfY,myTime,myThid) |
239 |
|
|
240 |
IF ( .NOT. BOTTOM_LAYER ) THEN |
IF ( .NOT. BOTTOM_LAYER ) THEN |
241 |
C-- Update fields in layer below according to tendency terms |
C-- Update fields in layer below according to tendency terms |
242 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
243 |
I bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,K+1,etaSurfX,etaSurfY,myTime,myThid) |
244 |
ENDIF |
ENDIF |
245 |
|
|
246 |
C-- Density of 1st level (below W(1)) reference to level 1 |
C-- Density of 1st level (below W(1)) reference to level 1 |
286 |
IF ( .NOT. BOTTOM_LAYER ) THEN |
IF ( .NOT. BOTTOM_LAYER ) THEN |
287 |
C-- Update fields in layer below according to tendency terms |
C-- Update fields in layer below according to tendency terms |
288 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
289 |
I bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,K+1,etaSurfX,etaSurfY,myTime,myThid) |
290 |
ENDIF |
ENDIF |
291 |
|
|
292 |
C-- Density of K level (below W(K)) reference to K level |
C-- Density of K level (below W(K)) reference to K level |