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