#include "dims.h"
      SUBROUTINE DIVB(DVB,J)
      use cons_module,only: imax,imaxp2,cs,dlamda,dphi,re
      implicit none
C     ****
C     ****      EVALUATES THE DIVERGENCE OF B, THE UNIT MAGNETIC FIELD
C     ****      VECTOR
C     ****
#include "params.h"
      real :: bxm,bx,bxp,by,byp,bz,bzp,bmod,bmodp
      COMMON/MAGFLD/BXM(zimxp,2),BX(zimxp,zjmx),BXP(zimxp,4),
     1  BY(zimxp,zjmx),BYP(zimxp,4),BZ(zimxp,zjmx),BZP(zimxp,4),
     2  BMOD(zimxp,zjmx),BMODP(zimxp,2)
!
! Args:
      real,intent(out) :: DVB(ZIMXP)
      integer,intent(in) :: j
!
! Local:
      integer :: i
!
      DO 1 I=3,IMAXP2
        DVB(I) = (((BX(I+1,J)-BX(I-1,J))/(2.*dlamda)+(CS(J+1)*BY(I,J+1)-
     A   CS(J-1)*BY(I,J-1))/(2.*dphi))/CS(J)+2.*BZ(I,J))/re
    1 CONTINUE
C     ****
C     ****      PERIODIC POINTS
C     ****
      DO 2 I=1,2
        DVB(I) = DVB(I+IMAX)
        DVB(I+IMAXP2) = DVB(I+2)
    2 CONTINUE
      RETURN
      END
C