55 |
& +omega3(i,j+1)*r_hFacZ(i,j+1)) |
& +omega3(i,j+1)*r_hFacZ(i,j+1)) |
56 |
ENDIF |
ENDIF |
57 |
uCoriolisTerm(i,j)= |
uCoriolisTerm(i,j)= |
|
C high order vorticity advection term |
|
58 |
& +vort3u*vBarXY*recip_dxc(i,j,bi,bj) |
& +vort3u*vBarXY*recip_dxc(i,j,bi,bj) |
59 |
|
& *_maskW(I,J,K,bi,bj) |
60 |
|
cph *note* put these comments after end of continued line |
61 |
|
cph to ensure TAMC compatibility |
62 |
|
C high order vorticity advection term |
63 |
|
c & ... |
64 |
C linear Coriolis term |
C linear Coriolis term |
65 |
c & +0.5*(fCoriG(I,J,bi,bj)+fCoriG(I,J+1,bi,bj))*vBarXY |
c & +0.5*(fCoriG(I,J,bi,bj)+fCoriG(I,J+1,bi,bj))*vBarXY |
66 |
C full nonlinear Coriolis term |
C full nonlinear Coriolis term |
70 |
c & fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj) ) |
c & fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj) ) |
71 |
C original form of Coriolis term (copied from calc_mom_rhs) |
C original form of Coriolis term (copied from calc_mom_rhs) |
72 |
c & +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY |
c & +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY |
|
& *_maskW(I,J,K,bi,bj) |
|
73 |
ENDDO |
ENDDO |
74 |
ENDDO |
ENDDO |
75 |
|
|