88 |
C- not using SadournyCoriolis discretization: |
C- not using SadournyCoriolis discretization: |
89 |
|
|
90 |
uBarXY=0.25*( |
uBarXY=0.25*( |
91 |
& uFld( i , j )*dyG( i , j ,bi,bj)*hFacW( i , j ,k,bi,bj) |
& (uFld( i , j )*dyG( i , j ,bi,bj)*hFacW( i , j ,k,bi,bj) |
92 |
& +uFld(i+1, j )*dyG(i+1, j ,bi,bj)*hFacW(i+1, j ,k,bi,bj) |
& +uFld( i ,j-1)*dyG( i ,j-1,bi,bj)*hFacW( i ,j-1,k,bi,bj)) |
93 |
& +uFld( i ,j-1)*dyG( i ,j-1,bi,bj)*hFacW( i ,j-1,k,bi,bj) |
& +(uFld(i+1, j )*dyG(i+1, j ,bi,bj)*hFacW(i+1, j ,k,bi,bj) |
94 |
& +uFld(i+1,j-1)*dyG(i+1,j-1,bi,bj)*hFacW(i+1,j-1,k,bi,bj)) |
& +uFld(i+1,j-1)*dyG(i+1,j-1,bi,bj)*hFacW(i+1,j-1,k,bi,bj)) |
95 |
|
& ) |
96 |
IF (upwindVort3) THEN |
IF (upwindVort3) THEN |
97 |
IF (uBarXY.GT.0.) THEN |
IF (uBarXY.GT.0.) THEN |
98 |
vort3v=omega3(i,j)*r_hFacZ(i,j) |
vort3v=omega3(i,j)*r_hFacZ(i,j) |