#include "dims.h"
      SUBROUTINE DIVB(DVB,J)
      use cons_module,only: imax,imaxp2,cs,re,dphi,dlamda
      implicit none
C     ****
C     ****      EVALUATES THE DIVERGENCE OF B, THE UNIT MAGNETIC FIELD
C     ****      VECTOR
C     ****
#include "params.h"
#include "magfld.h"
!
! Args:
      integer,intent(in) :: j
      real,intent(out) :: DVB(ZIMXP)
!
! Local:
      integer :: i

      DO I=3,IMAXP2
        DVB(I) = (((BX(I+1,J)-BX(I-1,J))/(2.*dlamda)+(CS(J+1)*BY(I,J+1)-
     |   CS(J-1)*BY(I,J-1))/(2.*dphi))/CS(J)+2.*BZ(I,J))/re
      enddo
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