C $Header: /u/gcmpack/MITgcm/pkg/mom_fluxform/mom_u_coriolis.F,v 1.9 2008/11/16 17:25:26 jmc Exp $ C $Name: $ #include "MOM_FLUXFORM_OPTIONS.h" CBOP C !ROUTINE: MOM_U_CORIOLIS C !INTERFACE: ========================================================== SUBROUTINE MOM_U_CORIOLIS( I bi,bj,k,vFld, U uCoriolisTerm, I myThid) C !DESCRIPTION: C Calculates the horizontal Coriolis term in the zonal equation: C \begin{equation*} C \overline{f}^i \overline{v}^{ij} C \end{equation*} C !USES: =============================================================== IMPLICIT NONE #include "SIZE.h" #include "EEPARAMS.h" #include "PARAMS.h" #include "GRID.h" #include "SURFACE.h" C !INPUT PARAMETERS: =================================================== C bi,bj :: tile indices C k :: vertical level C vFld :: meridional flow C myThid :: thread number INTEGER bi,bj,k _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) INTEGER myThid C !OUTPUT PARAMETERS: ================================================== C uCoriolisTerm :: Coriolis term _RL uCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy) C !LOCAL VARIABLES: ==================================================== C i,j :: loop indices INTEGER i,j _RS one PARAMETER( one = 1. _d 0 ) CEOP IF (useEnergyConservingCoriolis) THEN C Energy conserving discretization DO j=1-Oly,sNy+Oly-1 DO i=1-Olx+1,sNx+Olx uCoriolisTerm(i,j) = & 0.5*( _fCori( i ,j,bi,bj) & *0.5*( vFld( i ,j)+vFld( i ,j+1) ) & + _fCori(i-1,j,bi,bj) & *0.5*( vFld(i-1,j)+vFld(i-1,j+1) ) ) ENDDO ENDDO ELSE C Original discretization DO j=1-Oly,sNy+Oly-1 DO i=1-Olx+1,sNx+Olx uCoriolisTerm(i,j) = & 0.5*( _fCori( i ,j,bi,bj) + & _fCori(i-1,j,bi,bj) ) & *0.25*( & vFld( i ,j)+vFld( i ,j+1) & +vFld(i-1,j)+vFld(i-1,j+1) & ) ENDDO ENDDO ENDIF IF (useJamartWetPoints) THEN C Scale term so that only "wet" points are used C Due to: Jamart and Ozer, 1986, JGR 91 (C9), 10,621-10,631 C "Numerical Boundary Layers and Spurious Residual Flows" DO j=1-Oly,sNy+Oly-1 DO i=1-Olx+1,sNx+Olx uCoriolisTerm(i,j) = uCoriolisTerm(i,j) & *4. _d 0/MAX( one, & maskS( i , j ,k,bi,bj)+maskS( i ,j+1,k,bi,bj) & +maskS(i-1, j ,k,bi,bj)+maskS(i-1,j+1,k,bi,bj) ) ENDDO ENDDO ENDIF RETURN END