258 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
259 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
260 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
261 |
|
#ifdef ALLOW_OBCS |
262 |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K, myThid ) |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K, myThid ) |
263 |
|
#endif |
264 |
IF ( .NOT. BOTTOM_LAYER ) THEN |
IF ( .NOT. BOTTOM_LAYER ) THEN |
265 |
C-- Update fields in layer below according to tendency terms |
C-- Update fields in layer below according to tendency terms |
266 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
267 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
268 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
269 |
|
#ifdef ALLOW_OBCS |
270 |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
271 |
|
#endif |
272 |
ENDIF |
ENDIF |
273 |
#endif |
#endif |
274 |
C-- Density of 1st level (below W(1)) reference to level 1 |
C-- Density of 1st level (below W(1)) reference to level 1 |
325 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
326 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
327 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
328 |
|
#ifdef ALLOW_OBCS |
329 |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
330 |
|
#endif |
331 |
ENDIF |
ENDIF |
332 |
#endif |
#endif |
333 |
C-- Density of K level (below W(K)) reference to K level |
C-- Density of K level (below W(K)) reference to K level |
464 |
CALL TIMESTEP( |
CALL TIMESTEP( |
465 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
466 |
I myThid) |
I myThid) |
467 |
|
#ifdef ALLOW_OBCS |
468 |
C-- Apply open boundary conditions |
C-- Apply open boundary conditions |
469 |
IF (openBoundaries) CALL APPLY_OBCS2( bi, bj, K, myThid ) |
IF (openBoundaries) CALL APPLY_OBCS2( bi, bj, K, myThid ) |
470 |
|
#endif |
471 |
C-- Freeze water |
C-- Freeze water |
472 |
IF (allowFreezing) |
IF (allowFreezing) |
473 |
& CALL FREEZE( bi, bj, iMin, iMax, jMin, jMax, K, myThid ) |
& CALL FREEZE( bi, bj, iMin, iMax, jMin, jMax, K, myThid ) |