subroutine GetForces(x,a,N,size,periodic,m,neighbourlist,rpotcut, &
     morseDe,morsealpha,morseRe,Epot,virial)
  
  USE GLOBAL
  implicit none

  type(vector)  :: x(*),a(*)
  integer       :: N,neighbourlist(*)
  real(double)  :: m,rpotcut
  type(vector)  :: size
  type(ivector) :: periodic
  real(double)  :: morseDe,morsealpha,morseRe,Epot(*)
  real(double)  :: virial

  ! Stuff removed

end subroutine GetForces