#include "dims.h"
      SUBROUTINE VDOTDH(VX,VY,PHIJM1,PHIJ,PHIJP1,J,ANS)
      use cons_module,only: len2,kmax,imax,imaxp2,cs,dlamda,dphi,re
      implicit none
C     ****
C     ****      EVALUATES   ANS = (V(H).DEL(H))*PHI
C     ****
#include "params.h"
!
! Args:
      integer,intent(in) :: j
      real,intent(in) :: VX(ZIMXP,ZKMXP),VY(ZIMXP,ZKMXP),
     |  PHIJM1(ZIMXP,ZKMXP),PHIJ(ZIMXP,ZKMXP),PHIJP1(ZIMXP,ZKMXP)
      real,intent(out) :: ANS(ZIMXP,ZKMXP)
!
! Local:
      integer :: i,k
!
      DO 1 I=3,LEN2-2
        ANS(I,1) = 1./re*(VX(I,1)/(CS(J)*2.*dlamda)*(PHIJ(I+1,1)-
     A   PHIJ(I-1,1))+VY(I,1)/(2.*dphi)*(PHIJP1(I,1)-PHIJM1(I,1)))
    1 CONTINUE
C     ****
C     ****      PERIODIC POINTS
C     ****
      DO 2 I=1,2
        DO 3 K=1,KMAX
          ANS(I,K) = ANS(I+IMAX,K)
          ANS(I+IMAXP2,K) = ANS(I+2,K)
    3   CONTINUE
    2 CONTINUE
      RETURN
      END
C