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

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

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

revision 1.1 by adcroft, Wed Jun 21 19:45:50 2000 UTC revision 1.2 by adcroft, Thu Jun 29 17:39:12 2000 UTC
# Line 50  C     I, J, K - Loop counters Line 50  C     I, J, K - Loop counters
50        DO j=jMin,jMax        DO j=jMin,jMax
51         DO i=iMin,iMax         DO i=iMin,iMax
52          KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*          KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*
53       &                   (KPPdiffKzS(i,j,k,bi,bj)-diffKrT)       &                   (KPPdiffKzS(i,j,k,bi,bj)-diffKrS)
54         ENDDO         ENDDO
55        ENDDO        ENDDO
56    
57        DO j=jMin,jMax        DO j=jMin,jMax
58         DO i=iMin,iMax         DO i=iMin,iMax
59          KappaRU(i,j,k) = KappaRU(i,j,k)          KappaRU(i,j,k) = max ( KappaRU(i,j,k),
60       &       + 0.5*( KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj) )       &         KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *
61         &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )
62         ENDDO         ENDDO
63        ENDDO        ENDDO
64    
65        DO j=jMin,jMax        DO j=jMin,jMax
66         DO i=iMin,iMax         DO i=iMin,iMax
67          KappaRV(i,j,k) = KappaRV(i,j,k)          KappaRV(i,j,k) = max ( KappaRV(i,j,k),
68       &       + 0.5*( KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj) )       &         KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *
69         &         0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )
70         ENDDO         ENDDO
71        ENDDO        ENDDO
72    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22