anchor.diff:        Ep = 0._xp
   
         !$OMP PARALLEL REDUCTION(+:Ec,Ep) PRIVATE(j)
  -      !$OMP DO SCHEDULE(GUIDED)
  +      !$OMP DO SCHEDULE(RUNTIME)
         do i = 1, npoints
   
            Ec = Ec + 0.5_xp * m(i) * norm2(v(:,i))**2
  
positive.code: 
           end do
  
        end do
        !$OMP END DO
        !$OMP END PARALLEL
  
negative.code:       real(kind = xp), intent(in)  :: r(:,:)
        real(kind = xp), intent(out) :: v(:,:)
  
        v(1, :) =   r(2, :)
        v(2, :) = - r(1, :)
        v(3, :) = 0._xp
  
     end subroutine initial_speeds
  