anchor.diff:  ! y
   do i = 1, size(A)
       i1 = A(i)
  -    call set_next_eigeval(iunt, g, kpr, 0, i1, 0, w, dk, prg)
  +    call set_next_eigeval(iunt, g, kpc, 0, i1, 0, w, dk, prg)
   end do
   
   ! z
   do i = 1, size(A)
       i1 = A(i)
  -    call set_next_eigeval(iunt, g, kpr, 0, 0, i1, w, dk, prg)
  +    call set_next_eigeval(iunt, g, kpc, 0, 0, i1, w, dk, prg)
   end do
   
   ! xy
  
positive.code:     i1 = A(i)
      do j=1, size(A)
          j1 = A(j)
          call set_next_eigeval(iunt, g, kpr, 0, j1, i1, w, dk, prg)
      end do
  end do
  
  
negative.code: integer(kind=4), parameter :: iunt = 10, ilog = 11, w = 1 ! k point weight for VASP
  character(len=3), parameter :: version_number = '1.0'
  
  real(kind=8) :: kp(3), kpc(3), dk, E(-2:2,-2:2,-2:2), A(4), f(3,3), g(3,3)
  integer(kind=4) :: i, i1, j, j1, band
  character prg
  
  