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: |
45 |
|
|
46 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
47 |
C i,j :: loop indices |
C i,j :: loop indices |
48 |
INTEGER I,J |
INTEGER i,j |
49 |
CEOP |
CEOP |
|
_RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
|
50 |
|
|
51 |
IF ( k.EQ.Nr+1 .AND. |
IF ( k.EQ.Nr+1 .AND. |
52 |
& useRealFreshWaterFlux .AND. |
& useRealFreshWaterFlux .AND. usingPCoords ) THEN |
|
& buoyancyRelation .EQ. 'OCEANICP' ) 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)* |
83 |
|
#ifdef MOM_BOUNDARY_CONSERVE |
84 |
|
& 0.5 _d 0*( uFld(i,j,k,bi,bj)*_maskW(i,j,k-1,bi,bj) |
85 |
|
& +uFld(i,j,k-1,bi,bj)*_maskW(i,j,k,bi,bj) ) |
86 |
|
#else |
87 |
& 0.5 _d 0*( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) ) |
& 0.5 _d 0*( uFld(i,j,k,bi,bj)+uFld(i,j,k-1,bi,bj) ) |
88 |
|
#endif |
89 |
ENDDO |
ENDDO |
90 |
ENDDO |
ENDDO |
91 |
|
|
92 |
IF ( select_rStar.EQ.0 .AND. .NOT.rigidLid ) THEN |
IF ( select_rStar.EQ.0 .AND. .NOT.rigidLid ) THEN |
93 |
|
c & .AND. usingPCoords ) THEN |
94 |
C (linear) Free-surface correction at k>1 |
C (linear) Free-surface correction at k>1 |
95 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
96 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
97 |
advectiveFluxWU(i,j) = advectiveFluxWU(i,j) |
advectiveFluxWU(i,j) = advectiveFluxWU(i,j) |
98 |
& +0.25*( |
& +0.25*( |
99 |
& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)* |
& wFld(i, j ,k,bi,bj)*rA(i, j ,bi,bj)* |
100 |
& (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)) |
101 |
& +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)* |
102 |
& (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)) |
103 |
& )*uFld(i,j,k,bi,bj) |
& )*deepFac2F(k)*rhoFacF(k) |
104 |
|
& *uFld(i,j,k,bi,bj) |
105 |
ENDDO |
ENDDO |
106 |
ENDDO |
ENDDO |
107 |
C- endif NOT rigidLid |
C- endif NOT rigidLid |