7 |
C !ROUTINE: MOM_W_CORIOLIS_NH |
C !ROUTINE: MOM_W_CORIOLIS_NH |
8 |
|
|
9 |
C !INTERFACE: ========================================================== |
C !INTERFACE: ========================================================== |
10 |
SUBROUTINE MOM_W_CORIOLIS_NH( |
SUBROUTINE MOM_W_CORIOLIS_NH( |
11 |
I bi,bj,k, |
I bi,bj,k, |
12 |
I uFld, vFld, |
I uFld, vFld, |
13 |
U wCoriolisTerm, |
U wCoriolisTerm, |
14 |
I myThid) |
I myThid ) |
15 |
|
|
16 |
C !DESCRIPTION: |
C !DESCRIPTION: |
17 |
C Calculates the Coriolis term in the vertical momentum equation: |
C Calculates the Coriolis term in the vertical momentum equation: |
52 |
IF ( k.GT.1 .AND. k.LE.Nr ) THEN |
IF ( k.GT.1 .AND. k.LE.Nr ) THEN |
53 |
DO j=1-Oly,sNy+Oly-1 |
DO j=1-Oly,sNy+Oly-1 |
54 |
DO i=1-Olx,sNx+Olx-1 |
DO i=1-Olx,sNx+Olx-1 |
55 |
wCoriolisTerm(i,j) = |
wCoriolisTerm(i,j) = |
56 |
& -gravitySign*fCoriCos(i,j,bi,bj)* |
& -gravitySign*fCoriCos(i,j,bi,bj)* |
57 |
& ( angleCosC(i,j,bi,bj)*( |
& ( angleCosC(i,j,bi,bj)*( |
58 |
& (uFld(i,j,k-1,bi,bj)+uFld(i+1,j,k-1,bi,bj)) |
& (uFld(i,j,k-1,bi,bj)+uFld(i+1,j,k-1,bi,bj)) |
62 |
& (vFld(i,j,k-1,bi,bj)+vFld(i,j+1,k-1,bi,bj)) |
& (vFld(i,j,k-1,bi,bj)+vFld(i,j+1,k-1,bi,bj)) |
63 |
& + (vFld(i,j, k ,bi,bj)+vFld(i,j+1, k ,bi,bj)) |
& + (vFld(i,j, k ,bi,bj)+vFld(i,j+1, k ,bi,bj)) |
64 |
& )*0.25 _d 0 |
& )*0.25 _d 0 |
65 |
& )*horiVertRatio |
& )*wUnit2rVel(k) |
66 |
ENDDO |
ENDDO |
67 |
ENDDO |
ENDDO |
68 |
ELSE |
ELSE |