/[MITgcm]/MITgcm/pkg/mom_common/mom_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_calc_visc.F

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

revision 1.2 by baylor, Mon Sep 19 20:46:23 2005 UTC revision 1.3 by baylor, Tue Sep 20 13:09:10 2005 UTC
# Line 162  c This approximation is good to the same Line 162  c This approximation is good to the same
162            Asmag2=smag2fac*L2            Asmag2=smag2fac*L2
163       &           *sqrt(tension(i,j)**2       &           *sqrt(tension(i,j)**2
164       &           +0.25*(strain(i+1, j )**2+strain( i ,j+1)**2       &           +0.25*(strain(i+1, j )**2+strain( i ,j+1)**2
165       &                 +strain(i-1, j )**2+strain( i ,j-1)**2))       &                 +strain(i  , j )**2+strain(i+1,j+1)**2))
166           ELSE           ELSE
167            Asmag2=0d0            Asmag2=0d0
168           ENDIF           ENDIF
# Line 170  c This approximation is good to the same Line 170  c This approximation is good to the same
170           IF (smag4fac.NE.0.) THEN           IF (smag4fac.NE.0.) THEN
171            Asmag4=smag4fac*L4            Asmag4=smag4fac*L4
172       &             *sqrt(tension(i,j)**2       &             *sqrt(tension(i,j)**2
173       &             +0.25*(strain(i+1, j )**2+strain( i ,j+1)**2       &           +0.25*(strain(i+1, j )**2+strain( i ,j+1)**2
174       &                   +strain(i-1, j )**2+strain( i ,j-1)**2))       &                 +strain(i  , j )**2+strain(i+1,j+1)**2))
175           ELSE           ELSE
176            Asmag4=0d0            Asmag4=0d0
177           ENDIF           ENDIF
# Line 275  C This if statement is just to prevent b Line 275  C This if statement is just to prevent b
275           IF (smag2fac.NE.0.) THEN           IF (smag2fac.NE.0.) THEN
276            Asmag2=smag2fac*L2            Asmag2=smag2fac*L2
277       &           *sqrt(strain(i,j)**2       &           *sqrt(strain(i,j)**2
278       &              +0.25*(tension( i , j )**2+tension( i ,j+1)**2       &              +0.25*(tension( i , j )**2+tension( i ,j-1)**2
279       &                +tension(i+1, j )**2+tension(i+1,j+1)**2))       &                    +tension(i-1, j )**2+tension(i-1,j-1)**2))
280           ELSE           ELSE
281            Asmag2=0d0            Asmag2=0d0
282           ENDIF           ENDIF
# Line 284  C This if statement is just to prevent b Line 284  C This if statement is just to prevent b
284           IF (smag4fac.NE.0.) THEN           IF (smag4fac.NE.0.) THEN
285            Asmag4=smag4fac*L4            Asmag4=smag4fac*L4
286       &             *sqrt(strain(i,j)**2       &             *sqrt(strain(i,j)**2
287       &             +0.25*(tension( i , j )**2+tension( i ,j+1)**2       &              +0.25*(tension( i , j )**2+tension( i ,j-1)**2
288       &                +tension(i+1, j )**2+tension(i+1,j+1)**2))       &                    +tension(i-1, j )**2+tension(i-1,j-1)**2))
289           ELSE           ELSE
290            Asmag4=0d0            Asmag4=0d0
291           ENDIF           ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22