lmd_vmix_mod Module Reference

Functions/Subroutines
subroutine, public	lmd_vmix (ng, tile)
subroutine	lmd_vmix_tile (ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, nstp, hz, rho, u, v, t, alfaobeta, bvf, akt, akv)
subroutine	lmd_finish (ng, tile)
subroutine	lmd_finish_tile (ng, tile, lbi, ubi, lbj, ubj, nstp, hz, bvf, akt, akv)

Function/Subroutine Documentation

lmd_finish()

subroutine lmd_vmix_mod::lmd_finish ( integer, intent(in) ng, integer, intent(in) tile )

private

Definition at line 437 of file lmd_vmix.F.

!***********************************************************************
!
      USE mod_param
      USE mod_grid
      USE mod_mixing
      USE mod_ocean
      USE mod_stepping
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
!
!  Local variable declarations.
!
# include "tile.h"
!
      CALL lmd_finish_tile (ng, tile,                                   &
     &                      lbi, ubi, lbj, ubj,                         &
     &                      nstp(ng),                                   &
     &                      grid(ng) % Hz,                              &
     &                      mixing(ng) % bvf,                           &
     &                      mixing(ng) % Akt,                           &
     &                      mixing(ng) % Akv)
      RETURN

References mod_grid::grid, lmd_finish_tile(), mod_mixing::mixing, and mod_stepping::nstp.

Referenced by lmd_vmix().

Here is the call graph for this function:

Here is the caller graph for this function:

lmd_finish_tile()

subroutine lmd_vmix_mod::lmd_finish_tile ( integer, intent(in) ng, integer, intent(in) tile, integer, intent(in) lbi, integer, intent(in) ubi, integer, intent(in) lbj, integer, intent(in) ubj, integer, intent(in) nstp, real(r8), dimension(lbi:,lbj:,:), intent(in) hz, real(r8), dimension(lbi:,lbj:,0:), intent(in) bvf, real(r8), dimension(lbi:,lbj:,0:,:), intent(inout) akt, real(r8), dimension(lbi:,lbj:,0:), intent(inout) akv )

private

Definition at line 465 of file lmd_vmix.F.

!***********************************************************************
!
      USE mod_param
      USE mod_scalars
!
      USE bc_3d_mod, ONLY : bc_w3d_tile
# ifdef DISTRIBUTE
      USE mp_exchange_mod, ONLY : mp_exchange3d, mp_exchange4d
# endif
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
      integer, intent(in) :: LBi, UBi, LBj, UBj
      Integer, intent(in) :: nstp
!
# ifdef ASSUMED_SHAPE
      real(r8), intent(in) :: Hz(LBi:,LBj:,:)
      real(r8), intent(in) :: bvf(LBi:,LBj:,0:)
 
      real(r8), intent(inout) :: Akt(LBi:,LBj:,0:,:)
      real(r8), intent(inout) :: Akv(LBi:,LBj:,0:)
# else
      real(r8), intent(in) :: Hz(LBi:UBi,LBj:UBj,N(ng))
      real(r8), intent(in) :: bvf(LBi:UBi,LBj:UBj,0:N(ng))
 
      real(r8), intent(inout) :: Akt(LBi:UBi,LBj:UBj,0:N(ng),NAT)
      real(r8), intent(inout) :: Akv(LBi:UBi,LBj:UBj,0:N(ng))
# endif
!
!  Local variable declarations.
!
      integer :: i, itrc, j, k
 
      real(r8), parameter :: eps = 1.0e-14_r8
 
      real(r8) :: cff, lmd_iwm, lmd_iws, nu_sx, nu_sxc, shear2
 
# include "set_bounds.h"
!
!-----------------------------------------------------------------------
!  Compute "interior" viscosities and diffusivities everywhere as
!  the superposition of three processes: local Richardson number
!  instability due to resolved vertical shear, internal wave
!  breaking, and double diffusion.
!-----------------------------------------------------------------------
!
      DO k=1,n(ng)-1
        DO j=jstr,jend
          DO i=istr,iend
!
!  Compute interior convective diffusivity due to static instability
!  mixing.
!
# ifdef LMD_CONVEC
            cff=max(bvf(i,j,k),lmd_bvfcon)
            cff=min(1.0_r8,(lmd_bvfcon-cff)/lmd_bvfcon)
            nu_sxc=1.0_r8-cff*cff
            nu_sxc=nu_sxc*nu_sxc*nu_sxc
# else
            nu_sxc=0.0_r8
# endif
!
! Sum contributions due to internal wave breaking, shear instability
! and convective diffusivity due to shear instability.
!
            akv(i,j,k)=akv(i,j,k)+lmd_nu0c*nu_sxc
            akt(i,j,k,itemp)=akt(i,j,k,itemp)+lmd_nu0c*nu_sxc
# ifdef SALINITY
            akt(i,j,k,isalt)=akt(i,j,k,isalt)+lmd_nu0c*nu_sxc
# endif
# if defined LIMIT_VDIFF || defined LIMIT_VVISC
!
!  Limit vertical mixing coefficients with the upper threshold value.
!  This is an engineering fix but it can be based on the fact that
!  vertical mixing in the ocean from indirect observations is not
!  higher than the threshold value.
!
#  ifdef LIMIT_VVISC
            akv(i,j,k)=min(akv_limit(ng), akv(i,j,k))
#  endif
#  ifdef LIMIT_VDIFF
            akt(i,j,k,itemp)=min(akt_limit(itemp,ng), akt(i,j,k,itemp))
#   ifdef SALINITY
            akt(i,j,k,isalt)=min(akt_limit(isalt,ng), akt(i,j,k,isalt))
#   endif
#  endif
# endif
          END DO
        END DO
      END DO
!
!  Apply boundary conditions.
!
      DO k=0,n(ng)
        IF (domain(ng)%Western_Edge(tile)) THEN
          DO j=jstr,jend
            DO itrc=1,nat
              akt(istr-1,j,k,itrc)=akt(istr,j,k,itrc)
            END DO
            akv(istr-1,j,k)=akv(istr,j,k)
          END DO
        END IF
        IF (domain(ng)%Eastern_Edge(tile)) THEN
          DO j=jstr,jend
            DO itrc=1,nat
              akt(iend+1,j,k,itrc)=akt(iend,j,k,itrc)
            END DO
            akv(iend+1,j,k)=akv(iend,j,k)
          END DO
 
        END IF
        IF (domain(ng)%Southern_Edge(tile)) THEN
          DO i=istr,iend
            DO itrc=1,nat
              akt(i,jstr-1,k,itrc)=akt(i,jstr,k,itrc)
            END DO
            akv(i,jstr-1,k)=akv(i,jstr,k)
          END DO
        END IF
        IF (domain(ng)%Northern_Edge(tile)) THEN
          DO i=istr,iend
            DO itrc=1,nat
              akt(i,jend+1,k,itrc)=akt(i,jend,k,itrc)
            END DO
            akv(i,jend+1,k)=akv(i,jend,k)
          END DO
        END IF
        IF (domain(ng)%SouthWest_Corner(tile)) THEN
          DO itrc=1,nat
            akt(istr-1,jstr-1,k,itrc)=0.5_r8*                           &
     &                                (akt(istr  ,jstr-1,k,itrc)+       &
     &                                 akt(istr-1,jstr  ,k,itrc))
          END DO
          akv(istr-1,jstr-1,k)=0.5_r8*                                  &
     &                         (akv(istr  ,jstr-1,k)+                   &
     &                          akv(istr-1,jstr  ,k))
        END IF
        IF (domain(ng)%SouthEast_Corner(tile)) THEN
          DO itrc=1,nat
            akt(iend+1,jstr-1,k,itrc)=0.5_r8*                           &
     &                                (akt(iend  ,jstr-1,k,itrc)+       &
     &                                 akt(iend+1,jstr  ,k,itrc))
          END DO
          akv(iend+1,jstr-1,k)=0.5_r8*                                  &
     &                         (akv(iend  ,jstr-1,k)+                   &
     &                          akv(iend+1,jstr  ,k))
        END IF
        IF (domain(ng)%NorthWest_Corner(tile)) THEN
          DO itrc=1,nat
            akt(istr-1,jend+1,k,itrc)=0.5_r8*                           &
     &                                (akt(istr  ,jend+1,k,itrc)+       &
     &                                 akt(istr-1,jend  ,k,itrc))
          END DO
          akv(istr-1,jend+1,k)=0.5_r8*                                  &
     &                         (akv(istr  ,jend+1,k)+                   &
     &                          akv(istr-1,jend  ,k))
        END IF
        IF (domain(ng)%NorthEast_Corner(tile)) THEN
          DO itrc=1,nat
            akt(iend+1,jend+1,k,itrc)=0.5_r8*                           &
     &                                (akt(iend  ,jend+1,k,itrc)+       &
     &                                 akt(iend+1,jend  ,k,itrc))
          END DO
          akv(iend+1,jend+1,k)=0.5_r8*                                  &
     &                         (akv(iend  ,jend+1,k)+                   &
     &                          akv(iend+1,jend  ,k))
        END IF
      END DO
!
      CALL bc_w3d_tile (ng, tile,                                       &
     &                  lbi, ubi, lbj, ubj, 0, n(ng),                   &
     &                  akv)
      DO itrc=1,nat
        CALL bc_w3d_tile (ng, tile,                                     &
     &                    lbi, ubi, lbj, ubj, 0, n(ng),                 &
     &                    akt(:,:,:,itrc))
      END DO
# ifdef DISTRIBUTE
      CALL mp_exchange3d (ng, tile, inlm, 1,                            &
     &                    lbi, ubi, lbj, ubj, 0, n(ng),                 &
     &                    nghostpoints,                                 &
     &                    ewperiodic(ng), nsperiodic(ng),               &
     &                    akv)
      CALL mp_exchange4d (ng, tile, inlm, 1,                            &
     &                    lbi, ubi, lbj, ubj, 0, n(ng), 1, nat,         &
     &                    nghostpoints,                                 &
     &                    ewperiodic(ng), nsperiodic(ng),               &
     &                    akt)
# endif
!
      RETURN

References mod_scalars::akt_limit, mod_scalars::akv_limit, bc_3d_mod::bc_w3d_tile(), mod_param::domain, mod_scalars::ewperiodic, mod_param::inlm, mod_scalars::isalt, mod_scalars::itemp, mod_scalars::lmd_bvfcon, mod_scalars::lmd_nu0c, mp_exchange_mod::mp_exchange3d(), mp_exchange_mod::mp_exchange4d(), mod_param::nghostpoints, and mod_scalars::nsperiodic.

Referenced by lmd_finish().

Here is the call graph for this function:

Here is the caller graph for this function:

lmd_vmix()

subroutine, public lmd_vmix_mod::lmd_vmix	(	integer, intent(in)	ng,
		integer, intent(in)	tile )

Definition at line 33 of file lmd_vmix.F.

!***********************************************************************
!
      USE mod_param
      USE mod_grid
      USE mod_mixing
      USE mod_ocean
      USE mod_scalars
      USE mod_stepping
!
# ifdef LMD_SKPP
      USE lmd_skpp_mod, ONLY : lmd_skpp
# endif
# ifdef LMD_BKPP
      USE lmd_bkpp_mod, ONLY : lmd_bkpp
# endif
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
!
!  Local variable declarations.
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__
!
# include "tile.h"
!
# ifdef PROFILE
      CALL wclock_on (ng, inlm, 18, __line__, myfile)
# endif
      IF (.not.perfectrst(ng).or.iic(ng).ne.ntstart(ng)) THEN
        CALL lmd_vmix_tile (ng, tile,                                   &
     &                      lbi, ubi, lbj, ubj,                         &
     &                      imins, imaxs, jmins, jmaxs,                 &
     &                      nstp(ng),                                   &
     &                      grid(ng) % Hz,                              &
# ifndef RI_SPLINES
     &                      grid(ng) % z_r,                             &
# endif
     &                      ocean(ng) % rho,                            &
     &                      ocean(ng) % u,                              &
     &                      ocean(ng) % v,                              &
# ifdef LMD_DDMIX
     &                      ocean(ng) % t,                              &
     &                      mixing(ng) % alfaobeta,                     &
# endif
     &                      mixing(ng) % bvf,                           &
     &                      mixing(ng) % Akt,                           &
     &                      mixing(ng) % Akv)
# ifdef LMD_SKPP
        CALL lmd_skpp (ng, tile)
# endif
# ifdef LMD_BKPP
        CALL lmd_bkpp (ng, tile)
# endif
        CALL lmd_finish (ng, tile)
      END IF
# ifdef PROFILE
      CALL wclock_off (ng, inlm, 18, __line__, myfile)
# endif
!
      RETURN

References mod_grid::grid, mod_scalars::iic, mod_param::inlm, lmd_bkpp(), lmd_finish(), lmd_skpp(), lmd_vmix_tile(), mod_mixing::mixing, mod_stepping::nstp, mod_scalars::ntstart, mod_ocean::ocean, mod_scalars::perfectrst, wclock_off(), and wclock_on().

Referenced by main3d().

Here is the call graph for this function:

Here is the caller graph for this function:

lmd_vmix_tile()

subroutine lmd_vmix_mod::lmd_vmix_tile ( integer, intent(in) ng, integer, intent(in) tile, integer, intent(in) lbi, integer, intent(in) ubi, integer, intent(in) lbj, integer, intent(in) ubj, integer, intent(in) imins, integer, intent(in) imaxs, integer, intent(in) jmins, integer, intent(in) jmaxs, integer, intent(in) nstp, real(r8), dimension(lbi:,lbj:,:), intent(in) hz, real(r8), dimension(lbi:,lbj:,:), intent(in) rho, real(r8), dimension(lbi:,lbj:,:,:), intent(in) u, real(r8), dimension(lbi:,lbj:,:,:), intent(in) v, real(r8), dimension(lbi:,lbj:,:,:,:), intent(in) t, real(r8), dimension(lbi:,lbj:,0:), intent(in) alfaobeta, real(r8), dimension(lbi:,lbj:,0:), intent(in) bvf, real(r8), dimension(lbi:,lbj:,0:,:), intent(inout) akt, real(r8), dimension(lbi:,lbj:,0:), intent(inout) akv )

private

Definition at line 99 of file lmd_vmix.F.

!***********************************************************************
!
      USE mod_param
      USE mod_scalars
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
      integer, intent(in) :: LBi, UBi, LBj, UBj
      integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
      Integer, intent(in) :: nstp
!
# ifdef ASSUMED_SHAPE
      real(r8), intent(in) :: Hz(LBi:,LBj:,:)
#  ifndef RI_SPLINES
      real(r8), intent(in) :: z_r(LBi:,LBj:,:)
#  endif
      real(r8), intent(in) :: rho(LBi:,LBj:,:)
      real(r8), intent(in) :: u(LBi:,LBj:,:,:)
      real(r8), intent(in) :: v(LBi:,LBj:,:,:)
#  ifdef LMD_DDMIX
      real(r8), intent(in) :: t(LBi:,LBj:,:,:,:)
      real(r8), intent(in) :: alfaobeta(LBi:,LBj:,0:)
#  endif
      real(r8), intent(in) :: bvf(LBi:,LBj:,0:)
 
      real(r8), intent(inout) :: Akt(LBi:,LBj:,0:,:)
      real(r8), intent(inout) :: Akv(LBi:,LBj:,0:)
# else
      real(r8), intent(in) :: Hz(LBi:UBi,LBj:UBj,N(ng))
#  ifndef RI_SPLINES
      real(r8), intent(in) :: z_r(LBi:UBi,LBj:UBj,N(ng))
#  endif
      real(r8), intent(in) :: rho(LBi:UBi,LBj:UBj,N(ng))
      real(r8), intent(in) :: u(LBi:UBi,LBj:UBj,N(ng),3)
      real(r8), intent(in) :: v(LBi:UBi,LBj:UBj,N(ng),3)
#  ifdef LMD_DDMIX
      real(r8), intent(in) :: t(LBi:UBi,LBj:UBj,N(ng),3,NT(ng))
      real(r8), intent(in) :: alfaobeta(LBi:UBi,LBj:UBj,0:N(ng))
#  endif
      real(r8), intent(in) :: bvf(LBi:UBi,LBj:UBj,0:N(ng))
 
      real(r8), intent(inout) :: Akt(LBi:UBi,LBj:UBj,0:N(ng),NAT)
      real(r8), intent(inout) :: Akv(LBi:UBi,LBj:UBj,0:N(ng))
# endif
!
!  Local variable declarations.
!
      integer :: i, itrc, j, k
 
      real(r8), parameter :: eps = 1.0e-14_r8
 
      real(r8) :: cff, lmd_iwm, lmd_iws, nu_sx, nu_sxc, shear2
# ifdef LMD_DDMIX
      real(r8) :: Rrho, ddDS, ddDT, nu_dds, nu_ddt
# endif
 
      real(r8), dimension(IminS:ImaxS,JminS:JmaxS,0:N(ng)) :: Rig
 
      real(r8), dimension(IminS:ImaxS,0:N(ng)) :: FC
      real(r8), dimension(IminS:ImaxS,0:N(ng)) :: dR
      real(r8), dimension(IminS:ImaxS,0:N(ng)) :: dU
      real(r8), dimension(IminS:ImaxS,0:N(ng)) :: dV
 
# include "set_bounds.h"
 
# ifdef LMD_RIMIX
!
!-----------------------------------------------------------------------
! Compute gradient Richardson number.
!-----------------------------------------------------------------------
!
!  Compute gradient Richardson number at horizontal RHO-points and
!  vertical W-points.  If zero or very small velocity shear, bound
!  computation by a large negative value.
!
#  ifdef RI_SPLINES
      DO j=max(1,jstr-1),min(jend+1,mm(ng))
        DO i=max(1,istr-1),min(iend+1,lm(ng))
          fc(i,0)=0.0_r8
          dr(i,0)=0.0_r8
          du(i,0)=0.0_r8
          dv(i,0)=0.0_r8
        END DO
        DO k=1,n(ng)-1
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            cff=1.0_r8/(2.0_r8*hz(i,j,k+1)+                             &
     &                  hz(i,j,k)*(2.0_r8-fc(i,k-1)))
            fc(i,k)=cff*hz(i,j,k+1)
            dr(i,k)=cff*(6.0_r8*(rho(i,j,k+1)-rho(i,j,k))-              &
     &                   hz(i,j,k)*dr(i,k-1))
            du(i,k)=cff*(3.0_r8*(u(i  ,j,k+1,nstp)-u(i  ,j,k,nstp)+     &
     &                           u(i+1,j,k+1,nstp)-u(i+1,j,k,nstp))-    &
     &                   hz(i,j,k)*du(i,k-1))
            dv(i,k)=cff*(3.0_r8*(v(i,j  ,k+1,nstp)-v(i,j  ,k,nstp)+     &
     &                           v(i,j+1,k+1,nstp)-v(i,j+1,k,nstp))-    &
     &                   hz(i,j,k)*dv(i,k-1))
          END DO
        END DO
        DO i=max(1,istr-1),min(iend+1,lm(ng))
          dr(i,n(ng))=0.0_r8
          du(i,n(ng))=0.0_r8
          dv(i,n(ng))=0.0_r8
        END DO
        DO k=n(ng)-1,1,-1
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            dr(i,k)=dr(i,k)-fc(i,k)*dr(i,k+1)
            du(i,k)=du(i,k)-fc(i,k)*du(i,k+1)
            dv(i,k)=dv(i,k)-fc(i,k)*dv(i,k+1)
          END DO
        END DO
        DO k=1,n(ng)-1
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            shear2=du(i,k)*du(i,k)+dv(i,k)*dv(i,k)
            rig(i,j,k)=bvf(i,j,k)/(shear2+eps)
!!          Rig(i,j,k)=-gorho0*dR(i,k)/(shear2+eps)
          END DO
        END DO
      END DO
#  else
      DO k=1,n(ng)-1
        DO j=max(1,jstr-1),min(jend+1,mm(ng))
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            cff=0.5_r8/(z_r(i,j,k+1)-z_r(i,j,k))
            shear2=(cff*(u(i  ,j,k+1,nstp)-u(i  ,j,k,nstp)+             &
     &                   u(i+1,j,k+1,nstp)-u(i+1,j,k,nstp)))**2+        &
     &             (cff*(v(i,j  ,k+1,nstp)-v(i,j  ,k,nstp)+             &
     &                   v(i,j+1,k+1,nstp)-v(i,j+1,k,nstp)))**2
            rig(i,j,k)=bvf(i,j,k)/(shear2+eps)
          END DO
        END DO
      END DO
#  endif
#  ifdef RI_HORAVG
      DO k=1,n(ng)-1
        IF (domain(ng)%Western_Edge(tile)) THEN
          DO j=max(1,jstr-1),min(jend+1,mm(ng))
            rig(istr-1,j,k)=rig(istr,j,k)
          END DO
        END IF
        IF (domain(ng)%Eastern_Edge(tile)) THEN
          DO j=max(1,jstr-1),min(jend+1,mm(ng))
            rig(iend+1,j,k)=rig(iend,j,k)
          END DO
        END IF
        IF (domain(ng)%Southern_Edge(tile)) THEN
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            rig(i,jstr-1,k)=rig(i,jstr,k)
          END DO
        END IF
        IF (domain(ng)%Northern_Edge(tile)) THEN
          DO i=max(1,istr-1),min(iend+1,lm(ng))
            rig(i,jend+1,k)=rig(i,jend,k)
          END DO
        END IF
        IF (domain(ng)%SouthWest_Corner(tile)) THEN
          rig(istr-1,jstr-1,k)=rig(istr,jstr,k)
        END IF
        IF (domain(ng)%NorthWest_Corner(tile)) THEN
          rig(istr-1,jend+1,k)=rig(istr,jend,k)
        END IF
        IF (domain(ng)%SouthEast_Corner(tile)) THEN
          rig(iend+1,jstr-1,k)=rig(iend,jstr,k)
        END IF
        IF (domain(ng)%NorthEast_Corner(tile)) THEN
          rig(iend+1,jend+1,k)=rig(iend,jend,k)
        END IF
!
!  Smooth gradient Richardson number horizontally.  Use Rig(:,:,0)
!  as scratch utility array.
!
        DO j=jstr-1,jend
          DO i=istr-1,iend
            rig(i,j,0)=0.25_r8*(rig(i,j  ,k)+rig(i+1,j  ,k)+            &
     &                          rig(i,j+1,k)+rig(i+1,j+1,k))
          END DO
        END DO
        DO j=jstr,jend
          DO i=istr,iend
            rig(i,j,k)=0.25_r8*(rig(i,j  ,0)+rig(i-1,j  ,0)+            &
     &                          rig(i,j-1,0)+rig(i-1,j-1,0))
          END DO
        END DO
      END DO
#  endif
#  ifdef RI_VERAVG
!
!  Smooth gradient Richardson number vertically at the interior points.
!
      DO k=n(ng)-2,2,-1
        DO j=jstr,jend
          DO i=istr,iend
            rig(i,j,k)=0.25_r8*rig(i,j,k-1)+                            &
     &                 0.50_r8*rig(i,j,k  )+                            &
     &                 0.25_r8*rig(i,j,k+1)
          END DO
        END DO
      END DO
#  endif
# endif
!
!-----------------------------------------------------------------------
!  Compute "interior" viscosities and diffusivities everywhere as
!  the superposition of three processes: local Richardson number
!  instability due to resolved vertical shear, internal wave
!  breaking, and double diffusion.
!-----------------------------------------------------------------------
!
      DO k=1,n(ng)-1
        DO j=jstr,jend
          DO i=istr,iend
!
!  Compute interior diffusivity due to shear instability mixing.
!
# ifdef LMD_RIMIX
            cff=min(1.0_r8,max(0.0_r8,rig(i,j,k))/lmd_ri0)
            nu_sx=1.0_r8-cff*cff
            nu_sx=nu_sx*nu_sx*nu_sx
!
!  The shear mixing should be also a function of the actual magnitude
!  of the shear, see Polzin (1996, JPO, 1409-1425).
!
            shear2=bvf(i,j,k)/(rig(i,j,k)+eps)
            cff=shear2*shear2/(shear2*shear2+16.0e-10_r8)
            nu_sx=cff*nu_sx
# else
            nu_sx=0.0_r8
# endif
!
!  Compute interior diffusivity due to wave breaking (Gargett and
!  Holloway.
!
            cff=1.0_r8/sqrt(max(bvf(i,j,k),1.0e-7_r8))
            lmd_iwm=1.0e-6_r8*cff
            lmd_iws=1.0e-7_r8*cff
!           lmd_iwm=lmd_nuwm
!           lmd_iws=lmd_nuws
!
! Sum contributions due to internal wave breaking, shear instability
! and convective diffusivity due to shear instability.
!
            akv(i,j,k)=lmd_iwm+lmd_nu0m*nu_sx
            akt(i,j,k,itemp)=lmd_iws+lmd_nu0s*nu_sx
# ifdef SALINITY
            akt(i,j,k,isalt)=akt(i,j,k,itemp)
# endif
          END DO
        END DO
# ifdef LMD_DDMIX
!
!-----------------------------------------------------------------------
!  Compute double-diffusive mixing.  It can occur when vertical
!  gradient of density is stable but the vertical gradient of
!  salinity (salt figering) or temperature (diffusive convection)
!  is unstable.
!-----------------------------------------------------------------------
!
!  Compute double-diffusive density ratio, Rrho.
!
        DO j=jstr,jend
          DO i=istr,iend
            dddt=t(i,j,k+1,nstp,itemp)-t(i,j,k,nstp,itemp)
            ddds=t(i,j,k+1,nstp,isalt)-t(i,j,k,nstp,isalt)
            ddds=sign(1.0_r8,ddds)*max(abs(ddds),1.0e-14_r8)
            rrho=alfaobeta(i,j,k)*dddt/ddds
!
!  Salt fingering case.
!
            IF ((rrho.gt.1.0_r8).and.(ddds.gt.0.0_r8)) THEN
!
!  Compute interior diffusivity for double diffusive mixing of
!  salinity.  Upper bound "Rrho" by "Rrho0"; (lmd_Rrho0=1.9,
!  lmd_nuf=0.001).
!
              rrho=min(rrho,lmd_rrho0)
              nu_dds=1.0_r8-((rrho-1.0_r8)/(lmd_rrho0-1.0_r8))**2
              nu_dds=lmd_nuf*nu_dds*nu_dds*nu_dds
!
!  Compute interior diffusivity for double diffusive mixing
!  of temperature (lmd_fdd=0.7).
!
              nu_ddt=lmd_fdd*nu_dds
!
!  Diffusive convection case.
!
            ELSE IF ((0.0_r8.lt.rrho).and.(rrho.lt.1.0_r8).and.         &
     &              (ddds.lt.0.0_r8)) THEN
!
!  Compute interior diffusivity for double diffusive mixing of
!  temperature (Marmorino and Caldwell, 1976); (lmd_nu=1.5e-6,
!  lmd_tdd1=0.909, lmd_tdd2=4.6, lmd_tdd3=0.54).
!
              nu_ddt=lmd_nu*lmd_tdd1*                                   &
     &               exp(lmd_tdd2*                                      &
     &                   exp(-lmd_tdd3*((1.0_r8/rrho)-1.0_r8)))
!
!  Compute interior diffusivity for double diffusive mixing
!  of salinity (lmd_sdd1=0.15, lmd_sdd2=1.85, lmd_sdd3=0.85).
!
              IF (rrho.lt.0.5_r8) THEN
                nu_dds=nu_ddt*lmd_sdd1*rrho
              ELSE
                nu_dds=nu_ddt*(lmd_sdd2*rrho-lmd_sdd3)
              END IF
            ELSE
              nu_ddt=0.0_r8
              nu_dds=0.0_r8
            END IF
!
!  Add double diffusion contribution to temperature and salinity
!  mixing coefficients.
!
            akt(i,j,k,itemp)=akt(i,j,k,itemp)+nu_ddt
#  ifdef SALINITY
            akt(i,j,k,isalt)=akt(i,j,k,isalt)+nu_dds
#  endif
          END DO
        END DO
# endif
      END DO
 
      RETURN

References mod_param::domain, mod_scalars::isalt, mod_scalars::itemp, mod_param::lm, mod_scalars::lmd_fdd, mod_scalars::lmd_nu, mod_scalars::lmd_nu0m, mod_scalars::lmd_nu0s, mod_scalars::lmd_nuf, mod_scalars::lmd_ri0, mod_scalars::lmd_rrho0, mod_scalars::lmd_sdd1, mod_scalars::lmd_sdd2, mod_scalars::lmd_sdd3, mod_scalars::lmd_tdd1, mod_scalars::lmd_tdd2, mod_scalars::lmd_tdd3, and mod_param::mm.

Referenced by lmd_vmix().

Here is the caller graph for this function: