8 |
I myThid) |
I myThid) |
9 |
|
|
10 |
CBOP |
CBOP |
11 |
C /==========================================================\ |
C *==========================================================* |
12 |
C | SUBROUTINE MY82_CALC_VISC | |
C | SUBROUTINE MY82_CALC_VISC |
13 |
C | o Add contrubution to net viscosity from MY82 mixing | |
C | o Add contrubution to net viscosity from MY82 mixing |
14 |
C \==========================================================/ |
C *==========================================================* |
15 |
IMPLICIT NONE |
IMPLICIT NONE |
16 |
|
|
17 |
C == GLobal variables == |
C == GLobal variables == |
39 |
DO j=jMin,jMax |
DO j=jMin,jMax |
40 |
DO i=iMin,iMax |
DO i=iMin,iMax |
41 |
KappaRU(i,j,k) = max ( KappaRU(i,j,k), |
KappaRU(i,j,k) = max ( KappaRU(i,j,k), |
42 |
& KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) * |
& KappaRU(i,j,k) - viscArNr(k) + _maskW(i,j,k,bi,bj) * |
43 |
& 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) ) |
& 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) ) |
44 |
ENDDO |
ENDDO |
45 |
ENDDO |
ENDDO |
47 |
DO j=jMin,jMax |
DO j=jMin,jMax |
48 |
DO i=iMin,iMax |
DO i=iMin,iMax |
49 |
KappaRV(i,j,k) = max ( KappaRV(i,j,k), |
KappaRV(i,j,k) = max ( KappaRV(i,j,k), |
50 |
& KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) * |
& KappaRV(i,j,k) - viscArNr(k) + _maskS(i,j,k,bi,bj) * |
51 |
& 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) ) |
& 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) ) |
52 |
ENDDO |
ENDDO |
53 |
ENDDO |
ENDDO |