/[MITgcm]/MITgcm/pkg/kpp/kpp_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_calc_visc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by jmc, Thu May 3 21:31:35 2007 UTC revision 1.4 by jmc, Thu Oct 8 20:06:22 2009 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "KPP_OPTIONS.h"  #include "KPP_OPTIONS.h"
5    
6        SUBROUTINE KPP_CALC_VISC(        SUBROUTINE KPP_CALC_VISC(
7       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
8       U        KappaRU,KappaRV,       U        KappaRU,KappaRV,
9       I        myThid)       I        myThid)
10    
11  C     /==========================================================\  C     *==========================================================*
12  C     | SUBROUTINE KPP_CALC_VISC                                 |  C     | SUBROUTINE KPP_CALC_VISC                                 |
13  C     | o Add contrubution to net viscosity from KPP mixing      |  C     | o Add contrubution to net viscosity from KPP mixing      |
14  C     \==========================================================/  C     *==========================================================*
15        IMPLICIT NONE        IMPLICIT NONE
16    
17  C     == GLobal variables ==  C     == GLobal variables ==
# Line 40  C     I, J, K - Loop counters Line 40  C     I, J, K - Loop counters
40        DO j=jMin,jMax        DO j=jMin,jMax
41         DO i=iMin,iMax         DO i=iMin,iMax
42          KappaRU(i,j,k) = max ( KappaRU(i,j,k),          KappaRU(i,j,k) = max ( KappaRU(i,j,k),
43       &         KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *       &         KappaRU(i,j,k) - viscArNr(1) + _maskW(i,j,k,bi,bj) *
44       &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )       &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )
45         ENDDO         ENDDO
46        ENDDO        ENDDO
# Line 48  C     I, J, K - Loop counters Line 48  C     I, J, K - Loop counters
48        DO j=jMin,jMax        DO j=jMin,jMax
49         DO i=iMin,iMax         DO i=iMin,iMax
50          KappaRV(i,j,k) = max ( KappaRV(i,j,k),          KappaRV(i,j,k) = max ( KappaRV(i,j,k),
51       &         KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *       &         KappaRV(i,j,k) - viscArNr(1) + _maskS(i,j,k,bi,bj) *
52       &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )       &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )
53         ENDDO         ENDDO
54        ENDDO        ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22