10 |
SUBROUTINE MOMENTUM_CORRECTION_STEP(myTime, myIter, myThid) |
SUBROUTINE MOMENTUM_CORRECTION_STEP(myTime, myIter, myThid) |
11 |
C !DESCRIPTION: \bv |
C !DESCRIPTION: \bv |
12 |
C *==========================================================* |
C *==========================================================* |
13 |
C | SUBROUTINE MOMENTUM_CORRECTION_STEP |
C | SUBROUTINE MOMENTUM_CORRECTION_STEP |
14 |
C *==========================================================* |
C *==========================================================* |
15 |
C |1rst Part : Update U,V. |
C |1rst Part : Update U,V. |
16 |
C | |
C | |
41 |
#ifdef ALLOW_SHAP_FILT |
#ifdef ALLOW_SHAP_FILT |
42 |
#include "SHAP_FILT.h" |
#include "SHAP_FILT.h" |
43 |
#endif |
#endif |
44 |
#ifdef ALLOW_ZONAL_FILT |
#ifdef ALLOW_ZONAL_FILT |
45 |
#include "ZONAL_FILT.h" |
#include "ZONAL_FILT.h" |
46 |
#endif |
#endif |
47 |
|
|
117 |
C-- Filter (and exchange) |
C-- Filter (and exchange) |
118 |
#ifdef ALLOW_SHAP_FILT |
#ifdef ALLOW_SHAP_FILT |
119 |
IF (useSHAP_FILT) THEN |
IF (useSHAP_FILT) THEN |
120 |
IF ( .NOT.shap_filt_uvStar ) |
IF ( .NOT.shap_filt_uvStar ) THEN |
121 |
& CALL SHAP_FILT_APPLY_UV( uVel, vVel, myTime, myIter, myThid ) |
CALL TIMER_START('SHAP_FILT_UV [MOM_CORR_STEP]',myThid) |
122 |
|
CALL SHAP_FILT_APPLY_UV( uVel, vVel, myTime, myIter, myThid ) |
123 |
|
CALL TIMER_STOP ('SHAP_FILT_UV [MOM_CORR_STEP]',myThid) |
124 |
|
ENDIF |
125 |
ENDIF |
ENDIF |
126 |
#endif |
#endif |
127 |
#ifdef ALLOW_ZONAL_FILT |
#ifdef ALLOW_ZONAL_FILT |
128 |
IF (useZONAL_FILT) THEN |
IF (useZONAL_FILT) THEN |
129 |
IF ( .NOT.zonal_filt_uvStar ) |
IF ( .NOT.zonal_filt_uvStar ) THEN |
130 |
& CALL ZONAL_FILT_APPLY_UV( uVel, vVel, myThid ) |
CALL TIMER_START('ZONAL_FILT_UV [MOM_CORR_STEP]',myThid) |
131 |
|
CALL ZONAL_FILT_APPLY_UV( uVel, vVel, myThid ) |
132 |
|
CALL TIMER_STOP ('ZONAL_FILT_UV [MOM_CORR_STEP]',myThid) |
133 |
|
ENDIF |
134 |
ENDIF |
ENDIF |
135 |
#endif |
#endif |
136 |
|
|
137 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
138 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |