inimy.F90

      subroutine inimy
      use mod_xc         ! HYCOM communication interface
      use mod_cb_arrays  ! HYCOM saved arrays
!
! --- hycom version 2.1
      implicit none
!
! -------------------------------------------------------------
! --- initialize mellor-yamada level 2.5 vertical mixing scheme
! -------------------------------------------------------------
!
      integer i,j,k
!
# include "stmt_fns.h"
!
      a1my = 0.92
      b1my = 16.6
      a2my = 0.74
      b2my = 10.1
      c1my = 0.08
!
      e1my = 1.8
      e2my = 1.33
      e3my = 1.0
      sef  = 1.0
      smll = 1.e-8
!
! --- 'vonk'           = von karman constant
! --- 'ghc'            = constant for calculating tke production
! --- 'coef4','coef5'  = coefficients for estimating viscosity/diffusivity
! --- 'const1'         = constant for estimating surface and bottom bc's
!
      vonk   =  0.4
      ghc    = -6.0
      coef4  =  18.*a1my*a1my+9.*a1my*a2my
      coef5  =  9.*a1my*a2my
      const1 =  16.6**.6666667*sef
!
      do j=1,jdm
        do i=1,idm
          do k=0,kdm+1
            q2(i,j,k,1)=smll
            q2(i,j,k,2)=smll
            q2l(i,j,k,1)=smll
            q2l(i,j,k,2)=smll
            vctymy(i,j,k)=diwm(i,j)
            diftmy(i,j,k)=diws(i,j)
            difqmy(i,j,k)=diws(i,j)
          enddo
          do k=1,kdm+1
            vcty(i,j,k)=diwm(i,j)
            dift(i,j,k)=diws(i,j)
            difs(i,j,k)=diws(i,j)
! --- no nonlocal forcing
            ghats(i,j,k)=0.0
          enddo
        enddo
      enddo
      return
      end
!
!
!> Revision history:
!>