anchor.diff:      read(iunt,fmt=*) dk
       read(iunt,fmt=*) band
       read(iunt,fmt=*) prg
  -    if(prg .eq. 'C') then
  -        read(iunt,'(3e20.12)') ((f(i,j),j=1,3),i=1,3)
  -    else if(prg .eq. 'V') then
  -        read(iunt,*) ((f(i,j),j=1,3),i=1,3)
  -    end if
  +    read(iunt,*) ((f(i,j),j=1,3),i=1,3)
   close(iunt)
  +    
  +if(prg .eq. 'C') then
  +    f=f/(2.d0*pi)
  +end if
   
   write(ilog,"(A,I5,F10.6)") "band, dk: ", band, dk
   
  
positive.code:     i1 = A(i)
      do j=1, size(A)
          j1 = A(j)
          call set_next_eigeval(iunt, kpr, 0, j1, i1, w, dk)
      end do
  end do
  
  
negative.code: ! y
  do i = 1, size(A)
      i1 = A(i)
      call set_next_eigeval(iunt, g, kpr, 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)
  end do
  
  ! xy
  