pkg/kpp/kpp_calc_visc.F

77006ca295 Patr*0001 #include "KPP_OPTIONS.h"
                0002
a9d2e4c565 Jean*0003       SUBROUTINE KPP_CALC_VISC(
77006ca295 Patr*0004      I        bi,bj,iMin,iMax,jMin,jMax,K,
                0005      U        KappaRU,KappaRV,
                0006      I        myThid)
                0007
a9d2e4c565 Jean*0008 C     *==========================================================*
77006ca295 Patr*0009 C     | SUBROUTINE KPP_CALC_VISC                                 |
                0010 C     | o Add contrubution to net viscosity from KPP mixing      |
a9d2e4c565 Jean*0011 C     *==========================================================*
77006ca295 Patr*0012       IMPLICIT NONE
                0013
                0014 C     == GLobal variables ==
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
                0018 #include "DYNVARS.h"
                0019 #include "GRID.h"
ee8598ea0f Patr*0020 #include "KPP.h"
77006ca295 Patr*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 C     == Local variables ==
                0032 C     I, J, K - Loop counters
                0033       INTEGER i,j
                0034
                0035       DO j=jMin,jMax
                0036        DO i=iMin,iMax
                0037         KappaRU(i,j,k) = max ( KappaRU(i,j,k),
a9d2e4c565 Jean*0038      &         KappaRU(i,j,k) - viscArNr(1) + _maskW(i,j,k,bi,bj) *
77006ca295 Patr*0039      &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )
                0040        ENDDO
                0041       ENDDO
                0042
                0043       DO j=jMin,jMax
                0044        DO i=iMin,iMax
                0045         KappaRV(i,j,k) = max ( KappaRV(i,j,k),
a9d2e4c565 Jean*0046      &         KappaRV(i,j,k) - viscArNr(1) + _maskS(i,j,k,bi,bj) *
77006ca295 Patr*0047      &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )
                0048        ENDDO
                0049       ENDDO
                0050
                0051       RETURN
                0052       END