20 |
C | C*P* comments indicating place holders for which code is | |
C | C*P* comments indicating place holders for which code is | |
21 |
C | presently being developed. | |
C | presently being developed. | |
22 |
C \==========================================================/ |
C \==========================================================/ |
23 |
|
IMPLICIT NONE |
24 |
|
|
25 |
C == Global variables === |
C == Global variables === |
26 |
#include "SIZE.h" |
#include "SIZE.h" |
252 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
253 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
254 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
255 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K, myThid ) |
256 |
IF ( .NOT. BOTTOM_LAYER ) THEN |
IF ( .NOT. BOTTOM_LAYER ) THEN |
257 |
C-- Update fields in layer below according to tendency terms |
C-- Update fields in layer below according to tendency terms |
258 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
259 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
260 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
261 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
262 |
ENDIF |
ENDIF |
263 |
#endif |
#endif |
264 |
C-- Density of 1st level (below W(1)) reference to level 1 |
C-- Density of 1st level (below W(1)) reference to level 1 |
311 |
CALL CORRECTION_STEP( |
CALL CORRECTION_STEP( |
312 |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
I bi,bj,iMin,iMax,jMin,jMax,K+1, |
313 |
I etaSurfX,etaSurfY,myTime,myThid) |
I etaSurfX,etaSurfY,myTime,myThid) |
314 |
|
IF (openBoundaries) CALL APPLY_OBCS1( bi, bj, K+1, myThid ) |
315 |
ENDIF |
ENDIF |
316 |
#endif |
#endif |
317 |
C-- Density of K level (below W(K)) reference to K level |
C-- Density of K level (below W(K)) reference to K level |
432 |
CALL TIMESTEP( |
CALL TIMESTEP( |
433 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
434 |
I myThid) |
I myThid) |
435 |
|
IF (openBoundaries) CALL APPLY_OBCS2( bi, bj, K, myThid ) |
436 |
C-- Diagnose barotropic divergence of predicted fields |
C-- Diagnose barotropic divergence of predicted fields |
437 |
CALL CALC_DIV_GHAT( |
CALL CALC_DIV_GHAT( |
438 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |