c=======================================================================
c     this routine outputs the solution from twsplt.f at specified
c     times
c=======================================================================
      subroutine out(t, y)

      implicit none

c-----------------------------------------------------------------------
c     common communication with main program
c-----------------------------------------------------------------------
      include 'fcn.inc'

c-----------------------------------------------------------------------
c     n: number of odes
c     y: solution (y(1) = x, y(2) = p)
c     t: time independent variable
c-----------------------------------------------------------------------
      integer n
      parameter (n = 2)
      
      real t
      real y(n)

c-----------------------------------------------------------------------
c     write solution to std. out and increment output time
c-----------------------------------------------------------------------
      write (*,*) y(1), y(2)

      t = 0.0e0 + real(i)*per
      i = i + 1

      return

      end