/[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.6 by jmc, Sat Oct 1 19:17:17 2005 UTC revision 1.7 by heimbach, Wed Jun 7 01:55:15 2006 UTC
# Line 60  C-      using SadournyCoriolis discretiz Line 60  C-      using SadournyCoriolis discretiz
60    
61           vBarXY=1.           vBarXY=1.
62           vBarXm=0.5*(           vBarXm=0.5*(
63       &       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)
64       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj) )       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*_hFacS(i-1, j ,k,bi,bj) )
65           vBarXp=0.5*(           vBarXp=0.5*(
66       &       vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,k,bi,bj)       &       vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*_hFacS( i ,j+1,k,bi,bj)
67       &      +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) )
68           IF (upwindVorticity) THEN           IF (upwindVorticity) THEN
69            IF ( (vBarXm+vBarXp) .GT.0.) THEN            IF ( (vBarXm+vBarXp) .GT.0.) THEN
70             vort3u=vBarXm*r_hFacZ(i, j )*omega3(i, j )             vort3u=vBarXm*r_hFacZ(i, j )*omega3(i, j )
# 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-1, j )*dxG(i-1, j ,bi,bj)*hFacS(i-1, j ,k,bi,bj))       &      +vFld(i-1, j )*dxG(i-1, j ,bi,bj)*_hFacS(i-1, j ,k,bi,bj))
93       &     +(vFld( i ,j+1)*dxG( i ,j+1,bi,bj)*hFacS( i ,j+1,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
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22