pkg/kl10/kl10_calc_visc.F

d8d1486ca1 Jean*0001 #include "KL10_OPTIONS.h"
                0002
                0003       SUBROUTINE KL10_CALC_VISC(
                0004      I        bi,bj,iMin,iMax,jMin,jMax,k,
                0005      U        KappaRU,KappaRV,
                0006      I        myThid)
                0007
                0008 CBOP
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE KL10_CALC_VISC                                |
                0011 C     | o Add contrubution to net viscosity from KL10 mixing     |
                0012 C     *==========================================================*
                0013       IMPLICIT NONE
                0014
                0015 C     == GLobal variables ==
                0016 #include "SIZE.h"
                0017 #include "EEPARAMS.h"
                0018 #include "PARAMS.h"
                0019 #include "DYNVARS.h"
                0020 #include "GRID.h"
                0021 #include "KL10.h"
                0022
                0023 C     == Routine arguments ==
                0024 C     bi, bj :: tile numbers
                0025 C     iMin, iMax, jMin, jMax :: Range of points for which calculation
                0026 C     maskUp :: land/water mask for Wvel points (above tracer level)
                0027 C     myThid :: Instance number for this innvocation of KL10_CALC_VISC
                0028 C
                0029       INTEGER bi,bj,iMin,iMax,jMin,jMax,k
                0030       _RL KappaRU(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
                0031       _RL KappaRV(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
                0032       INTEGER myThid
                0033
                0034 #ifdef ALLOW_KL10
                0035
                0036 C     == Local variables ==
                0037 C     i, j    :: Loop counters
                0038       INTEGER i,j
                0039 CEOP
                0040       DO j=jMin,jMax
                0041        DO i=iMin,iMax
                0042         KappaRU(i,j,k) = max ( KappaRU(i,j,k),
                0043      &         KappaRU(i,j,k) - viscArNr(k) + _maskW(i,j,k,bi,bj) *
                0044      &         0.5*(KLviscAr(i,j,k,bi,bj)+KLviscAr(i-1,j,k,bi,bj)) )
                0045        ENDDO
                0046       ENDDO
                0047
                0048       DO j=jMin,jMax
                0049        DO i=iMin,iMax
                0050         KappaRV(i,j,k) = max ( KappaRV(i,j,k),
                0051      &         KappaRV(i,j,k) - viscArNr(k) + _maskS(i,j,k,bi,bj) *
                0052      &         0.5*(KLviscAr(i,j,k,bi,bj)+KLviscAr(i,j-1,k,bi,bj)) )
                0053        ENDDO
                0054       ENDDO
                0055
                0056 #endif /* ALLOW_KL10 */
                0057
                0058       RETURN
                0059       END