55 |
& +omega3(i+1,j)*r_hFacZ(i+1,j)) |
& +omega3(i+1,j)*r_hFacZ(i+1,j)) |
56 |
ENDIF |
ENDIF |
57 |
vCoriolisTerm(i,j)= |
vCoriolisTerm(i,j)= |
|
C high order vorticity advection term |
|
58 |
& -vort3v*uBarXY*recip_dyc(i,j,bi,bj) |
& -vort3v*uBarXY*recip_dyc(i,j,bi,bj) |
59 |
|
& *_maskS(I,J,K,bi,bj) |
60 |
|
C high order vorticity advection term |
61 |
|
c & ... |
62 |
C linear Coriolis term |
C linear Coriolis term |
63 |
c & -0.5 *(fCoriG(I,J,bi,bj)+fCoriG(I+1,J,bi,bj))*uBarXY |
c & -0.5 *(fCoriG(I,J,bi,bj)+fCoriG(I+1,J,bi,bj))*uBarXY |
64 |
C full nonlinear Coriolis term |
C full nonlinear Coriolis term |
68 |
c & fCori(I,J-1,bi,bj)*uBarX(I,J-1,K,bi,bj) ) |
c & fCori(I,J-1,bi,bj)*uBarX(I,J-1,K,bi,bj) ) |
69 |
C original form of Coriolis term (copied from calc_mom_rhs) |
C original form of Coriolis term (copied from calc_mom_rhs) |
70 |
c & -0.5*(fCori(i,j,bi,bj)+fCori(i,j-1,bi,bj))*uBarXY |
c & -0.5*(fCori(i,j,bi,bj)+fCori(i,j-1,bi,bj))*uBarXY |
|
& *_maskS(I,J,K,bi,bj) |
|
71 |
ENDDO |
ENDDO |
72 |
ENDDO |
ENDDO |
73 |
|
|