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