| 295 |
C---- Zonal momentum equation starts here |
C---- Zonal momentum equation starts here |
| 296 |
|
|
| 297 |
C Bi-harmonic term del^2 U -> v4F |
C Bi-harmonic term del^2 U -> v4F |
| 298 |
IF (momViscosity) |
IF (momViscosity .AND. viscA4.NE.0. ) |
| 299 |
& CALL MOM_U_DEL2U(bi,bj,k,uFld,hFacZ,v4f,myThid) |
& CALL MOM_U_DEL2U(bi,bj,k,uFld,hFacZ,v4f,myThid) |
| 300 |
|
|
| 301 |
C--- Calculate mean and eddy fluxes between cells for zonal flow. |
C--- Calculate mean and eddy fluxes between cells for zonal flow. |
| 328 |
& CALL MOM_U_YVISCFLUX(bi,bj,k,uFld,v4F,hFacZ,vF,myThid) |
& CALL MOM_U_YVISCFLUX(bi,bj,k,uFld,v4F,hFacZ,vF,myThid) |
| 329 |
|
|
| 330 |
C Combine fluxes -> fMer |
C Combine fluxes -> fMer |
| 331 |
DO j=jMin,jMax |
DO j=jMin,jMax+1 |
| 332 |
DO i=iMin,iMax |
DO i=iMin,iMax |
| 333 |
fMer(i,j) = vDudyFac*aF(i,j) + AhDudyFac*vF(i,j) |
fMer(i,j) = vDudyFac*aF(i,j) + AhDudyFac*vF(i,j) |
| 334 |
ENDDO |
ENDDO |
| 447 |
C---- Meridional momentum equation starts here |
C---- Meridional momentum equation starts here |
| 448 |
|
|
| 449 |
C Bi-harmonic term del^2 V -> v4F |
C Bi-harmonic term del^2 V -> v4F |
| 450 |
IF (momViscosity) |
IF (momViscosity .AND. viscA4.NE.0. ) |
| 451 |
& CALL MOM_V_DEL2V(bi,bj,k,vFld,hFacZ,v4f,myThid) |
& CALL MOM_V_DEL2V(bi,bj,k,vFld,hFacZ,v4f,myThid) |
| 452 |
|
|
| 453 |
C--- Calculate mean and eddy fluxes between cells for meridional flow. |
C--- Calculate mean and eddy fluxes between cells for meridional flow. |
| 464 |
|
|
| 465 |
C Combine fluxes -> fZon |
C Combine fluxes -> fZon |
| 466 |
DO j=jMin,jMax |
DO j=jMin,jMax |
| 467 |
DO i=iMin,iMax |
DO i=iMin,iMax+1 |
| 468 |
fZon(i,j) = uDvdxFac*aF(i,j) + AhDvdxFac*vF(i,j) |
fZon(i,j) = uDvdxFac*aF(i,j) + AhDvdxFac*vF(i,j) |
| 469 |
ENDDO |
ENDDO |
| 470 |
ENDDO |
ENDDO |