745 |
#endif |
#endif |
746 |
ENDIF |
ENDIF |
747 |
|
|
748 |
IF (nonHydrostatic.OR.quasiHydrostatic) THEN |
C-- 3.D Coriolis term (horizontal momentum, Eastward component: -f'*w) |
749 |
CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid) |
IF ( nonHydrostatic.OR.quasiHydrostatic ) THEN |
750 |
DO j=jMin,jMax |
CALL MOM_U_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid) |
751 |
DO i=iMin,iMax |
DO j=jMin,jMax |
752 |
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j) |
DO i=iMin,iMax |
753 |
|
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+fuFac*cf(i,j) |
754 |
|
ENDDO |
755 |
ENDDO |
ENDDO |
756 |
ENDDO |
IF ( usingCurvilinearGrid ) THEN |
757 |
|
C- presently, non zero angleSinC array only supported with Curvilinear-Grid |
758 |
|
CALL MOM_V_CORIOLIS_NH(bi,bj,k,wVel,cf,myThid) |
759 |
|
DO j=jMin,jMax |
760 |
|
DO i=iMin,iMax |
761 |
|
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+fvFac*cf(i,j) |
762 |
|
ENDDO |
763 |
|
ENDDO |
764 |
|
ENDIF |
765 |
ENDIF |
ENDIF |
766 |
|
|
767 |
C-- Set du/dt & dv/dt on boundaries to zero |
C-- Set du/dt & dv/dt on boundaries to zero |