subroutine UpdateNeighbourlist(x,N,size,periodic,neighbourlist,rskincut)
  USE GLOBAL
  implicit none

  ! Subroutine constructs a Verlet neighbour list in the
  ! following format:
  !
  ! nngbr1          Number of neighbours for atom 1
  ! ingbr1,1        Index of neighbour 1 of atom 1
  ! ingbr2,1        Index of neighbour 2 of atom 1
  ! ....
  ! ingbrnngbr1,1   Index of neighbour nngbr1 of atom 1
  ! nngbr2          Number of neighbours for atom 2
  ! ...

  
  type(vector)  :: x(*)
  integer       :: N,neighbourlist(*)
  type(vector)  :: size
  type(ivector) :: periodic
  real(double)  :: rskincut

  integer :: nneighbourstot

  
  
  ! Stuff removed

  
  print '(A,G13.5,A)','Neighbour list update found ',    &
       1.0*nneighbourstot/N,' neighbours per atom'
  
end subroutine UpdateNeighbourlist