subroutine Solve1(x,v,a,N,size,periodic,deltat)
  USE GLOBAL
  implicit none

  type(vector)  :: x(*),v(*),a(*)
  integer       :: N
  real(double)  :: deltat
  type(vector)  :: size
  type(ivector) :: periodic

  ! stuff removed
  
end subroutine Solve1


subroutine Solve2(v,a,N,deltat,btctau,T,desiredT)
  USE GLOBAL
  implicit none

  type(vector) :: v(*),a(*)
  integer      :: N
  real(double) :: deltat
  real(double) :: T,desiredT,btctau

  ! stuff removed
  
end subroutine Solve2