pkg/mom_common/mom_v_coriolis_nh.F

1a741329d2 Jean*0001 #include "MOM_COMMON_OPTIONS.h"
                0002
                0003 CBOP
                0004 C !ROUTINE: MOM_V_CORIOLIS_NH
                0005
                0006 C !INTERFACE: ==========================================================
                0007       SUBROUTINE MOM_V_CORIOLIS_NH(
                0008      I               bi,bj,k,wFld,
                0009      O               vCoriolisTerm,
                0010      I               myThid )
                0011
                0012 C !DESCRIPTION:
                0013 C Calculates the 3.D Coriolis term in the meridional momentum equation:
                0014 C \begin{equation*}
3daafce20b Jean*0015 C - \overline{ f_prime \overline{w}^{k} }^{i}
1a741329d2 Jean*0016 C \end{equation*}
                0017 C consistent with Non-Hydrostatic (or quasi-hydrostatic) formulation
                0018
                0019 C !USES: ===============================================================
                0020       IMPLICIT NONE
                0021 #include "SIZE.h"
                0022 #include "EEPARAMS.h"
                0023 #include "PARAMS.h"
                0024 #include "GRID.h"
                0025
                0026 C !INPUT PARAMETERS: ===================================================
                0027 C  bi,bj                :: tile indices
                0028 C  k                    :: vertical level
                0029 C  wFld                 :: vertical flow
                0030 C  myThid               :: thread number
8ecb2a8ef1 Jean*0031       INTEGER bi,bj,k
1a741329d2 Jean*0032       _RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0033       INTEGER myThid
                0034
                0035 C !OUTPUT PARAMETERS: ==================================================
                0036 C  vCoriolisTerm        :: Coriolis term
                0037       _RL vCoriolisTerm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0038
                0039 C !LOCAL VARIABLES: ====================================================
                0040 C  i,j                  :: loop indices
                0041       INTEGER i,j,kp1
                0042       _RL wMsk
                0043 CEOP
                0044
                0045
                0046       kp1=min(k+1,Nr)
                0047       wMsk=1.
                0048       IF (k.EQ.Nr) wMsk=0.
                0049
                0050 C Energy conserving discretization of 2*Omega*cos(phi)*w
                0051       DO j=2-Oly,sNy+Oly
                0052        DO i=1-Olx,sNx+Olx
                0053         vCoriolisTerm(i,j) =
                0054      &   -0.5*( fCoriCos( i ,j,bi,bj)*angleSinC(i, j ,bi,bj)
8ecb2a8ef1 Jean*0055      &         *0.5*( wFld( i ,j, k ,bi,bj)*rVel2wUnit( k )
                0056      &               +wFld( i ,j,kp1,bi,bj)*rVel2wUnit(kp1)*wMsk )
1a741329d2 Jean*0057      &        + fCoriCos(i,j-1,bi,bj)*angleSinC(i,j-1,bi,bj)
8ecb2a8ef1 Jean*0058      &         *0.5*( wFld(i,j-1, k ,bi,bj)*rVel2wUnit( k )
                0059      &               +wFld(i,j-1,kp1,bi,bj)*rVel2wUnit(kp1)*wMsk )
                0060      &        )*gravitySign
1a741329d2 Jean*0061        ENDDO
                0062       ENDDO
                0063
                0064       RETURN
                0065       END