--- MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis.F 2001/05/29 14:01:39 1.2 +++ MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis.F 2001/09/05 17:46:03 1.3 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis.F,v 1.2 2001/05/29 14:01:39 adcroft Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vi_u_coriolis.F,v 1.3 2001/09/05 17:46:03 heimbach Exp $ C $Name: $ #include "CPP_OPTIONS.h" @@ -55,8 +55,12 @@ & +omega3(i,j+1)*r_hFacZ(i,j+1)) ENDIF uCoriolisTerm(i,j)= -C high order vorticity advection term & +vort3u*vBarXY*recip_dxc(i,j,bi,bj) + & *_maskW(I,J,K,bi,bj) +cph *note* put these comments after end of continued line +cph to ensure TAMC compatibility +C high order vorticity advection term +c & ... C linear Coriolis term c & +0.5*(fCoriG(I,J,bi,bj)+fCoriG(I,J+1,bi,bj))*vBarXY C full nonlinear Coriolis term @@ -66,7 +70,6 @@ c & fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj) ) C original form of Coriolis term (copied from calc_mom_rhs) c & +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY - & *_maskW(I,J,K,bi,bj) ENDDO ENDDO