/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_fill_state.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_fill_state.F

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

revision 1.16 by jmc, Thu Feb 10 15:01:34 2005 UTC revision 1.17 by jmc, Thu Feb 10 16:30:39 2005 UTC
# Line 170  C---+----1----+----2----+----3----+----4 Line 170  C---+----1----+----2----+----3----+----4
170              DO j = 1,sNy              DO j = 1,sNy
171               DO i = 1,sNx+1               DO i = 1,sNx+1
172                tmpMk(i,j,K,bi,bj) = 0.25 _d 0                tmpMk(i,j,K,bi,bj) = 0.25 _d 0
173       &           *(uVel(i,km1,K,bi,bj)+uVel(i,j,K,bi,bj))       &           *(uVel(i,j,km1,bi,bj)+uVel(i,j,K,bi,bj))
174       &           *(wVel(i-1,j,K,bi,bj)*rA(i-1,j,bi,bj)       &           *(wVel(i-1,j,K,bi,bj)*rA(i-1,j,bi,bj)
175       &            +wVel( i ,j,K,bi,bj)*rA( i ,j,bi,bj)       &            +wVel( i ,j,K,bi,bj)*rA( i ,j,bi,bj)
176       &            )*recip_rAw(i,j,bi,bj)       &            )*recip_rAw(i,j,bi,bj)
# Line 190  C---+----1----+----2----+----3----+----4 Line 190  C---+----1----+----2----+----3----+----4
190              DO j = 1,sNy+1              DO j = 1,sNy+1
191               DO i = 1,sNx               DO i = 1,sNx
192                tmpMk(i,j,K,bi,bj) = 0.25 _d 0                tmpMk(i,j,K,bi,bj) = 0.25 _d 0
193       &           *(vVel(i,km1,K,bi,bj)+vVel(i,j,K,bi,bj))       &           *(vVel(i,j,km1,bi,bj)+vVel(i,j,K,bi,bj))
194       &           *(wVel(i,j-1,K,bi,bj)*rA(i,j-1,bi,bj)       &           *(wVel(i,j-1,K,bi,bj)*rA(i,j-1,bi,bj)
195       &            +wVel(i, j ,K,bi,bj)*rA(i, j ,bi,bj)       &            +wVel(i, j ,K,bi,bj)*rA(i, j ,bi,bj)
196       &            )*recip_rAs(i,j,bi,bj)       &            )*recip_rAs(i,j,bi,bj)
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22