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 |