!
  program kgen
  !
  ! generate a uniform mesh of kpoints in the first BZ [0,1)
  !
  implicit none
  integer, parameter :: n = 100
  integer :: i, j, k
  real(kind=8) :: x,y,z,w
  !
  write(*,*) n*n*n
  w = 2.d0/float(n*n*n)
  do k=1,n
  do i=1,n
  do j=1,n
     x = (float(i)-1.d0)/float(n)
     y = (float(j)-1.d0)/float(n)
     z = (float(k)-1.d0)/float(n)
     write(*,'4f20.15') x,y,z,w
  enddo
  enddo
  enddo
  !
  end
  !