/[MITgcm]/MITgcm/pkg/my82/my82_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/my82/my82_calc_visc.F

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

revision 1.3 by jmc, Tue Jan 20 00:26:04 2009 UTC revision 1.4 by jmc, Thu Oct 8 20:07:53 2009 UTC
# Line 8  C $Name$ Line 8  C $Name$
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 ==
# Line 39  CEOP Line 39  CEOP
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
# Line 47  CEOP Line 47  CEOP
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

Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4

  ViewVC Help
Powered by ViewVC 1.1.22