/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis_c4.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis_c4.F

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

revision 1.5 by jmc, Mon Jun 20 23:23:00 2005 UTC revision 1.6 by jmc, Sat Oct 1 19:17:17 2005 UTC
# Line 88  C-      using SadournyCoriolis discretiz Line 88  C-      using SadournyCoriolis discretiz
88  C-      not using SadournyCoriolis discretization:  C-      not using SadournyCoriolis discretization:
89    
90           vBarXY=0.25*(           vBarXY=0.25*(
91       &       vFld( i , j )*dxG( i , j ,bi,bj)*hFacS( i , j ,k,bi,bj)       &      (vFld( i , j )*dxG( i , j ,bi,bj)*hFacS( i , j ,k,bi,bj)
92       &      +vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj))
93       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj)       &     +(vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)
94       &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacS(i-1,j+1,k,bi,bj))       &      +vFld(i-1,j+1)*dxG(i-1,j+1,bi,bj)*hFacS(i-1,j+1,k,bi,bj))
95         &               )
96           IF (upwindVort3) THEN           IF (upwindVort3) THEN
97            IF (vBarXY.GT.0.) THEN            IF (vBarXY.GT.0.) THEN
98             vort3u=omega3(i,j)*r_hFacZ(i,j)             vort3u=omega3(i,j)*r_hFacZ(i,j)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22