| 17 |
C Calculates the area integrated vertical viscous fluxes of U |
C Calculates the area integrated vertical viscous fluxes of U |
| 18 |
C at vertical interface k (between level k & k-1): |
C at vertical interface k (between level k & k-1): |
| 19 |
C \begin{equation*} |
C \begin{equation*} |
| 20 |
C F^r = - \frac{ {\cal A}_w }{\Delta r_c} A_r \delta_k u |
C F^r = - \frac{ {\cal A}_w }{\Delta r_c} A_r \delta_k u |
| 21 |
C \end{equation*} |
C \end{equation*} |
| 22 |
|
|
| 23 |
C !USES: =============================================================== |
C !USES: =============================================================== |
| 59 |
DO i=1-Olx,sNx+Olx-1 |
DO i=1-Olx,sNx+Olx-1 |
| 60 |
rViscFluxU(i,j) = |
rViscFluxU(i,j) = |
| 61 |
& -KappaRU(i,j,k) |
& -KappaRU(i,j,k) |
| 62 |
& *rAw(i,j,bi,bj) |
& *rAw(i,j,bi,bj)*deepFac2F(k)*rhoFacF(k) |
| 63 |
& *( uFld(i,j,k,bi,bj)-uFld(i,j,k-1,bi,bj) |
& *( uFld(i,j,k,bi,bj)-uFld(i,j,k-1,bi,bj) |
| 64 |
& )*rkSign*recip_drC(k) |
& )*rkSign*recip_drC(k) |
| 65 |
& *_maskW(i,j,k,bi,bj) |
& *_maskW(i,j,k,bi,bj) |