11 |
I bi,bj,k, |
I bi,bj,k, |
12 |
I uFld,wFld,rTrans, |
I uFld,wFld,rTrans, |
13 |
O advectiveFluxWU, |
O advectiveFluxWU, |
14 |
I myThid) |
I myThid ) |
15 |
|
|
16 |
C !DESCRIPTION: |
C !DESCRIPTION: |
17 |
C Calculates the vertical advective flux of zonal momentum: |
C Calculates the vertical advective flux of zonal momentum: |
52 |
& useRealFreshWaterFlux .AND. usingPCoords ) THEN |
& useRealFreshWaterFlux .AND. usingPCoords ) THEN |
53 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
54 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
55 |
advectiveFluxWU(i,j) = rTrans(i,j)*uFld(i,j,k-1,bi,bj) |
advectiveFluxWU(i,j) = rTrans(i,j)*uFld(i,j,k-1,bi,bj) |
56 |
ENDDO |
ENDDO |
57 |
ENDDO |
ENDDO |
58 |
|
|
70 |
|
|
71 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
72 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
73 |
advectiveFluxWU(i,j) = rTrans(i,j)*uFld(i,j,k,bi,bj) |
advectiveFluxWU(i,j) = rTrans(i,j)*uFld(i,j,k,bi,bj) |
74 |
ENDDO |
ENDDO |
75 |
ENDDO |
ENDDO |
76 |
|
|
77 |
ELSE |
ELSE |
78 |
|
|
79 |
C Vertical advection - interior ; assume uFld & wFld are masked |
C Vertical advection - interior ; assume uFld & wFld are masked |
80 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
81 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
82 |
advectiveFluxWU(i,j) = rTrans(i,j)* |
advectiveFluxWU(i,j) = rTrans(i,j)* |
85 |
ENDDO |
ENDDO |
86 |
|
|
87 |
IF ( select_rStar.EQ.0 .AND. .NOT.rigidLid ) THEN |
IF ( select_rStar.EQ.0 .AND. .NOT.rigidLid ) THEN |
88 |
|
c & .AND. usingPCoords ) THEN |
89 |
C (linear) Free-surface correction at k>1 |
C (linear) Free-surface correction at k>1 |
90 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
91 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
92 |
advectiveFluxWU(i,j) = advectiveFluxWU(i,j) |
advectiveFluxWU(i,j) = advectiveFluxWU(i,j) |
93 |
& +0.25*( |
& +0.25*( |
94 |
& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)* |
& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)* |
95 |
& (maskC(i,j,k,bi,bj) - maskC(i,j,k-1,bi,bj)) |
& (maskC( i ,j,k,bi,bj)-maskC( i ,j,k-1,bi,bj)) |
96 |
& +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)* |
& +wFld(i-1,j,k,bi,bj)*rA(i-1,j,bi,bj)* |
97 |
& (maskC(i-1,j,k,bi,bj)-maskC(i-1,j,k-1,bi,bj)) |
& (maskC(i-1,j,k,bi,bj)-maskC(i-1,j,k-1,bi,bj)) |
98 |
& )*uFld(i,j,k,bi,bj) |
& )*deepFac2F(k)*rhoFacF(k) |
99 |
|
& *uFld(i,j,k,bi,bj) |
100 |
ENDDO |
ENDDO |
101 |
ENDDO |
ENDDO |
102 |
C- endif NOT rigidLid |
C- endif NOT rigidLid |