251 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
252 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
253 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
254 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K, myThid ) |
255 |
IF ( .NOT. BOTTOM_LAYER ) THEN |
IF ( .NOT. BOTTOM_LAYER ) THEN |
256 |
C-- Update fields in layer below according to tendency terms |
C-- Update fields in layer below according to tendency terms |
257 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
258 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
259 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
260 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
261 |
ENDIF |
ENDIF |
262 |
#endif |
#endif |
263 |
C-- Density of 1st level (below W(1)) reference to level 1 |
C-- Density of 1st level (below W(1)) reference to level 1 |
310 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
311 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
312 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
313 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
314 |
ENDIF |
ENDIF |
315 |
#endif |
#endif |
316 |
C-- Density of K level (below W(K)) reference to K level |
C-- Density of K level (below W(K)) reference to K level |
431 |
CALL TIMESTEP( |
CALL TIMESTEP( |
432 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
433 |
I myThid) |
I myThid) |
434 |
|
IF (openBoundaries) CALL APPLY_OBCS2( bi, bj, K, myThid ) |
435 |
C-- Diagnose barotropic divergence of predicted fields |
C-- Diagnose barotropic divergence of predicted fields |
436 |
CALL CALC_DIV_GHAT( |
CALL CALC_DIV_GHAT( |
437 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |