subroutine kenrgyp(u,wkdx,wkdy,wkin,wkix,wkiy,wki)
c
c....  volume integral of (grad u)**2 = 2.0*kinetic energy
c.... wkdx=volume integral of (d u/dx)**2
c.... wkdy=volume integral of (d u/dy)**2
c
      implicit none
      include 'itg.par'
      include 'itg.cmn'
c     include 'post.cmn'

      real u(lz,mz,nz,nspecz),wki(mz,nz),wkix(mz,nz),wkiy(mz,nz)
      real wkdx,wkdy,wkin

c Local variables:
      real coef,wx,wy,du1(lz,mz,nz,nspecz),du2(lz,mz,nz,nspecz)
      integer l,m,n
c
      call wgrad(u,du1,du2)
c
c     ... note that dr(l)=2.0*dx
c     ... normalize by delta=pi*y0 in the x-direction
c
      coef=0.25/(pi*y0)
      wkin=0.0
      wkdx=0.0
      wkdy=0.0
c
      do 10 n=1,nd
      do 10 m=1,md
      wki(m,n)=0.0
      wx=0.0
      wy=0.0
      wx=wx+dr(1)*du1(1,m,n,1)*du1(1,m,n,1)*0.5
      do 20 l=2,ldb
      wx=wx+dr(l)*du1(l,m,n,1)*du1(l,m,n,1)
      wy=wy+dr(l)*du2(l,m,n,1)*du2(l,m,n,1)
   20 continue
      wx=wx+dr(ld)*du1(ld,m,n,1)*du1(ld,m,n,1)*0.5
      wkix(m,n)=coef*wx
      wkiy(m,n)=coef*wy
      wki(m,n)=coef*(wx+wy)
      wkin=wkin+wki(m,n)
      wkdx=wkdx+coef*wx
      wkdy=wkdy+coef*wy
   10 continue
c
      return
      end