/[MITgcm]/MITgcm/verification/tutorial_deep_convection/code/mom_uv_smag_3d.F
ViewVC logotype

Diff of /MITgcm/verification/tutorial_deep_convection/code/mom_uv_smag_3d.F

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

revision 1.1 by jmc, Mon Sep 30 18:19:41 2013 UTC revision 1.2 by jmc, Tue Oct 1 14:38:56 2013 UTC
# Line 58  C     !LOCAL VARIABLES: Line 58  C     !LOCAL VARIABLES:
58        IF ( k.LE.1 )  maskM1 = 0.        IF ( k.LE.1 )  maskM1 = 0.
59        IF ( k.GE.Nr ) maskP1 = 0.        IF ( k.GE.Nr ) maskP1 = 0.
60    
61    C-    Initialise output arrays:
62          DO j=1-OLy,sNy+OLy
63           DO i=1-OLx,sNx+OLx
64             uDissip(i,j) = 0. _d 0
65             vDissip(i,j) = 0. _d 0
66           ENDDO
67          ENDDO
68    
69  C     - Laplacian and bi-harmonic terms  C     - Laplacian and bi-harmonic terms
70  c     IF (harmonic) THEN  c     IF (harmonic) THEN
71    
# Line 90  c    &         *maskW(i,j,km1,bi,bj)*mas Line 98  c    &         *maskW(i,j,km1,bi,bj)*mas
98    
99         DO j=2-OLy,sNy+OLy-1         DO j=2-OLy,sNy+OLy-1
100          DO i=2-OLx,sNx+OLx-1          DO i=2-OLx,sNx+OLx-1
101          vDissip(i,j) = recip_rAs(i,j,bi,bj)*(           vDissip(i,j) = recip_rAs(i,j,bi,bj)*(
102       &     ( viscAh3d_12(i+1, j , k )*str12(i+1, j ,k)       &     ( viscAh3d_12(i+1, j , k )*str12(i+1, j ,k)
103       &         *dyU(i+1, j ,bi,bj)       &         *dyU(i+1, j ,bi,bj)
104       &      -viscAh3d_12( i , j , k )*str12( i , j ,k)       &      -viscAh3d_12( i , j , k )*str12( i , j ,k)

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

  ViewVC Help
Powered by ViewVC 1.1.22