247 |
CEOP |
CEOP |
248 |
|
|
249 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
250 |
IF ( debugLevel .GE. debLevB ) |
IF (debugMode) CALL DEBUG_ENTER( 'DYNAMICS', myThid ) |
|
& CALL DEBUG_ENTER( 'DYNAMICS', myThid ) |
|
251 |
#endif |
#endif |
252 |
|
|
253 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
558 |
|
|
559 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
560 |
C-- Apply open boundary conditions |
C-- Apply open boundary conditions |
561 |
IF (useOBCS) THEN |
c IF (useOBCS) THEN |
562 |
CALL OBCS_APPLY_UV( bi, bj, k, gU, gV, myThid ) |
c CALL OBCS_APPLY_UV( bi, bj, k, gU, gV, myThid ) |
563 |
ENDIF |
c ENDIF |
564 |
#endif /* ALLOW_OBCS */ |
#endif /* ALLOW_OBCS */ |
565 |
|
|
566 |
ENDIF |
ENDIF |
567 |
|
|
|
|
|
568 |
C-- end of dynamics k loop (1:Nr) |
C-- end of dynamics k loop (1:Nr) |
569 |
ENDDO |
ENDDO |
570 |
|
|
602 |
|
|
603 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
604 |
C-- Apply open boundary conditions |
C-- Apply open boundary conditions |
605 |
IF ( useOBCS .AND.(implicitViscosity.OR.momImplVertAdv) ) THEN |
c IF ( useOBCS .AND.(implicitViscosity.OR.momImplVertAdv) ) THEN |
606 |
DO K=1,Nr |
c DO K=1,Nr |
607 |
CALL OBCS_APPLY_UV( bi, bj, k, gU, gV, myThid ) |
c CALL OBCS_APPLY_UV( bi, bj, k, gU, gV, myThid ) |
608 |
ENDDO |
c ENDDO |
609 |
|
IF ( useOBCS ) THEN |
610 |
|
CALL OBCS_APPLY_UV( bi, bj, 0, gU, gV, myThid ) |
611 |
ENDIF |
ENDIF |
612 |
#endif /* ALLOW_OBCS */ |
#endif /* ALLOW_OBCS */ |
613 |
|
|
639 |
C-- Step forward W field in N-H algorithm |
C-- Step forward W field in N-H algorithm |
640 |
IF ( nonHydrostatic ) THEN |
IF ( nonHydrostatic ) THEN |
641 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
642 |
IF ( debugLevel .GE. debLevB ) |
IF (debugMode) CALL DEBUG_CALL('CALC_GW', myThid ) |
|
& CALL DEBUG_CALL('CALC_GW', myThid ) |
|
643 |
#endif |
#endif |
644 |
CALL TIMER_START('CALC_GW [DYNAMICS]',myThid) |
CALL TIMER_START('CALC_GW [DYNAMICS]',myThid) |
645 |
CALL CALC_GW( |
CALL CALC_GW( |
689 |
#endif /* ALLOW_DIAGNOSTICS */ |
#endif /* ALLOW_DIAGNOSTICS */ |
690 |
|
|
691 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
692 |
If ( debugLevel .GE. debLevB ) THEN |
IF ( debugLevel .GE. debLevB ) THEN |
693 |
CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid) |
CALL DEBUG_STATS_RL(1,EtaN,'EtaN (DYNAMICS)',myThid) |
694 |
CALL DEBUG_STATS_RL(Nr,uVel,'Uvel (DYNAMICS)',myThid) |
CALL DEBUG_STATS_RL(Nr,uVel,'Uvel (DYNAMICS)',myThid) |
695 |
CALL DEBUG_STATS_RL(Nr,vVel,'Vvel (DYNAMICS)',myThid) |
CALL DEBUG_STATS_RL(Nr,vVel,'Vvel (DYNAMICS)',myThid) |
719 |
#endif |
#endif |
720 |
|
|
721 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
722 |
IF ( debugLevel .GE. debLevB ) |
IF (debugMode) CALL DEBUG_LEAVE( 'DYNAMICS', myThid ) |
|
& CALL DEBUG_LEAVE( 'DYNAMICS', myThid ) |
|
723 |
#endif |
#endif |
724 |
|
|
725 |
RETURN |
RETURN |