243 |
|
|
244 |
C Eddy component of vertical flux (interior component only) -> vrF |
C Eddy component of vertical flux (interior component only) -> vrF |
245 |
IF (momViscosity.AND..NOT.implicitViscosity) THEN |
IF (momViscosity.AND..NOT.implicitViscosity) THEN |
246 |
IF ( k.LT.Nr ) THEN |
CALL MOM_U_RVISCFLUX(bi,bj,k+1,uVel,KappaRU,vrF,myThid) |
|
CALL MOM_U_RVISCFLUX(bi,bj,k+1,uVel,KappaRU,vrF,myThid) |
|
|
ENDIF |
|
247 |
|
|
248 |
C Combine fluxes |
C Combine fluxes |
249 |
DO j=jMin,jMax |
DO j=jMin,jMax |
303 |
|
|
304 |
C Eddy component of vertical flux (interior component only) -> vrF |
C Eddy component of vertical flux (interior component only) -> vrF |
305 |
IF (momViscosity.AND..NOT.implicitViscosity) THEN |
IF (momViscosity.AND..NOT.implicitViscosity) THEN |
306 |
IF ( k.LT.Nr ) THEN |
CALL MOM_V_RVISCFLUX(bi,bj,k+1,vVel,KappaRV,vrF,myThid) |
|
CALL MOM_V_RVISCFLUX(bi,bj,k+1,vVel,KappaRV,vrF,myThid) |
|
|
ENDIF |
|
307 |
|
|
308 |
C Combine fluxes -> fVerV |
C Combine fluxes -> fVerV |
309 |
DO j=jMin,jMax |
DO j=jMin,jMax |
456 |
ENDIF |
ENDIF |
457 |
|
|
458 |
#ifdef ALLOW_TIMEAVE |
#ifdef ALLOW_TIMEAVE |
|
#ifndef MINIMAL_TAVE_OUTPUT |
|
459 |
IF (taveFreq.GT.0.) THEN |
IF (taveFreq.GT.0.) THEN |
460 |
CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock, |
CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock, |
461 |
& Nr, k, bi, bj, myThid) |
& Nr, k, bi, bj, myThid) |
462 |
CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock, |
CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock, |
463 |
& Nr, k, bi, bj, myThid) |
& Nr, k, bi, bj, myThid) |
464 |
ENDIF |
ENDIF |
|
#endif /* ndef MINIMAL_TAVE_OUTPUT */ |
|
465 |
#endif /* ALLOW_TIMEAVE */ |
#endif /* ALLOW_TIMEAVE */ |
466 |
|
|
467 |
C-- Vertical shear terms (-w*du/dr & -w*dv/dr) |
C-- Vertical shear terms (-w*du/dr & -w*dv/dr) |