/[MITgcm]/MITgcm/pkg/ggl90/ggl90_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/ggl90/ggl90_calc_visc.F

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

revision 1.2 by ce107, Tue Jul 31 23:01:58 2007 UTC revision 1.3 by dfer, Fri Oct 10 21:18:10 2008 UTC
# Line 39  C     I, J, K - Loop counters Line 39  C     I, J, K - Loop counters
39  CEOP  CEOP
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) = KappaRU(i,j,k) + _maskW(i,j,k,bi,bj) *
43       &     KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *       &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
44       &     0.5*(GGL90viscAr(i,j,k,bi,bj)+GGL90viscAr(i-1,j,k,bi,bj)) )       &                            +GGL90viscAr(i-1,j,k,bi,bj))
45         &                    -viscAr)
46         ENDDO         ENDDO
47        ENDDO        ENDDO
48    
49        DO j=jMin,jMax        DO j=jMin,jMax
50         DO i=iMin,iMax         DO i=iMin,iMax
51          KappaRV(i,j,k) = max ( KappaRV(i,j,k),          KappaRV(i,j,k) = KappaRV(i,j,k) + _maskS(i,j,k,bi,bj) *
52       &     KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *       &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
53       &     0.5*(GGL90viscAr(i,j,k,bi,bj)+GGL90viscAr(i,j-1,k,bi,bj)) )       &                            +GGL90viscAr(i,j-1,k,bi,bj))
54         &                    -viscAr)
55         ENDDO         ENDDO
56        ENDDO        ENDDO
57    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22