| 350 |
C--- endif momAdvection & k=1 |
C--- endif momAdvection & k=1 |
| 351 |
ENDIF |
ENDIF |
| 352 |
|
|
|
|
|
| 353 |
C--- Calculate vertical transports (at k+1) below U & V points : |
C--- Calculate vertical transports (at k+1) below U & V points : |
| 354 |
IF (momAdvection) THEN |
IF (momAdvection) THEN |
| 355 |
CALL MOM_CALC_RTRANS( k+1, bi, bj, |
CALL MOM_CALC_RTRANS( k+1, bi, bj, |
| 368 |
#endif /* MOM_BOUNDARY_CONSERVE */ |
#endif /* MOM_BOUNDARY_CONSERVE */ |
| 369 |
|
|
| 370 |
IF (momViscosity) THEN |
IF (momViscosity) THEN |
| 371 |
CALL MOM_CALC_VISC( bi, bj, k, |
DO j=1-OLy,sNy+OLy |
| 372 |
O viscAh_Z, viscAh_D, viscA4_Z, viscA4_D, |
DO i=1-OLx,sNx+OLx |
| 373 |
I hDiv, vort3, tension, strain, KE, hFacZ, |
viscAh_D(i,j) = viscAhD |
| 374 |
I myThid ) |
viscAh_Z(i,j) = viscAhZ |
| 375 |
|
viscA4_D(i,j) = viscA4D |
| 376 |
|
viscA4_Z(i,j) = viscA4Z |
| 377 |
|
ENDDO |
| 378 |
|
ENDDO |
| 379 |
|
IF ( useVariableVisc ) THEN |
| 380 |
|
CALL MOM_CALC_VISC( bi, bj, k, |
| 381 |
|
O viscAh_Z, viscAh_D, viscA4_Z, viscA4_D, |
| 382 |
|
I hDiv, vort3, tension, strain, KE, hFacZ, |
| 383 |
|
I myThid ) |
| 384 |
|
ENDIF |
| 385 |
ENDIF |
ENDIF |
| 386 |
|
|
| 387 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |