79 |
_RL pH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
_RL pH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nz) |
80 |
_RL rhokm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rhokm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
81 |
_RL rhokp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rhokp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
82 |
|
_RL pSurfX(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
83 |
|
_RL pSurfY(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
84 |
INTEGER iMin, iMax |
INTEGER iMin, iMax |
85 |
INTEGER jMin, jMax |
INTEGER jMin, jMax |
86 |
INTEGER bi, bj |
INTEGER bi, bj |
142 |
jMin = 1-OLy+1 |
jMin = 1-OLy+1 |
143 |
jMax = sNy+OLy |
jMax = sNy+OLy |
144 |
|
|
145 |
C-- Update fields according to tendency terms |
C-- Calculate gradient of surface pressure |
146 |
|
CALL GRAD_PSURF( |
147 |
|
I bi,bj,iMin,iMax,jMin,jMax, |
148 |
|
O pSurfX,pSurfY, |
149 |
|
I myThid) |
150 |
|
|
151 |
|
C-- Update fields in top level according to tendency terms |
152 |
CALL TIMESTEP( |
CALL TIMESTEP( |
153 |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,1,pSurfX,pSurfY,myThid) |
154 |
|
|
155 |
DO K=2,Nz |
DO K=2,Nz |
156 |
|
C-- Update fields in Kth level according to tendency terms |
157 |
|
CALL TIMESTEP( |
158 |
|
I bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myThid) |
159 |
C Density of K-1 level (above W(K)) reference to K level |
C Density of K-1 level (above W(K)) reference to K level |
160 |
CALL FIND_RHO( |
CALL FIND_RHO( |
161 |
I bi, bj, iMin, iMax, jMin, jMax, K-1, K, 'LINEAR', |
I bi, bj, iMin, iMax, jMin, jMax, K-1, K, 'LINEAR', |