pkg/my82/my82_calc_visc.F

08be60903a Mart*0001 #include "MY82_OPTIONS.h"
                0002
fb62f539dc Jean*0003       SUBROUTINE MY82_CALC_VISC(
08be60903a Mart*0004      I        bi,bj,iMin,iMax,jMin,jMax,K,
                0005      U        KappaRU,KappaRV,
                0006      I        myThid)
                0007
                0008 CBOP
5e48dccc42 Jean*0009 C     *==========================================================*
                0010 C     | SUBROUTINE MY82_CALC_VISC
                0011 C     | o Add contrubution to net viscosity from MY82 mixing
                0012 C     *==========================================================*
08be60903a Mart*0013       IMPLICIT NONE
                0014
                0015 C     == GLobal variables ==
                0016 #include "SIZE.h"
                0017 #include "EEPARAMS.h"
                0018 #include "PARAMS.h"
                0019 #include "GRID.h"
fb62f539dc Jean*0020 #include "MY82.h"
08be60903a Mart*0021
                0022 C     == Routine arguments ==
                0023 C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
                0024 C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
                0025 C
                0026       INTEGER bi,bj,iMin,iMax,jMin,jMax,K
                0027       _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
                0028       _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
                0029       INTEGER myThid
                0030
                0031 #ifdef ALLOW_MY82
                0032
                0033 C     == Local variables ==
                0034 C     I, J, K - Loop counters
                0035       INTEGER i,j
                0036 CEOP
                0037       DO j=jMin,jMax
                0038        DO i=iMin,iMax
                0039         KappaRU(i,j,k) = max ( KappaRU(i,j,k),
5e48dccc42 Jean*0040      &         KappaRU(i,j,k) - viscArNr(k) + _maskW(i,j,k,bi,bj) *
fb62f539dc Jean*0041      &    0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) )
08be60903a Mart*0042        ENDDO
                0043       ENDDO
                0044
                0045       DO j=jMin,jMax
                0046        DO i=iMin,iMax
                0047         KappaRV(i,j,k) = max ( KappaRV(i,j,k),
5e48dccc42 Jean*0048      &         KappaRV(i,j,k) - viscArNr(k) + _maskS(i,j,k,bi,bj) *
fb62f539dc Jean*0049      &    0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) )
08be60903a Mart*0050        ENDDO
                0051       ENDDO
                0052
                0053 #endif /* ALLOW_MY82 */
                0054
                0055       RETURN
                0056       END