program xephem_table_2
!
!     Exercise driver for subroutine ephem_table--call it many times
!     in a loop (for the Sun) to test performance.
!
!     B. Knapp, 2002-03-19
!
      implicit none

      integer*4 i, n
      parameter (n=100000)

      integer*4 p
      parameter (p=0) !the Sun

      integer*4 posnid, framid, coorid
      parameter (posnid=2) !Apparent position
      parameter (framid=1) !Equatorial frame
      parameter (coorid=2) !Rectangular coordinates

      real*8 jd, jd0, pv(6), dpv(6)
      parameter (jd0=2447892.5d0) !1990 Jan 1.0
      integer*4 ierr
      real*4 dtime, tarray(2)

!     Externals
      external ephem_table

!     Call once to get the initialization done
      write(*,*) dtime(tarray)
      call ephem_table(jd0, p, posnid, framid, coorid, pv, dpv, ierr)

      write(*,*) dtime(tarray)
      do i=0,n-1
         jd = jd0 + dble(i)/dble(n)
         call ephem_table(jd, p, posnid, framid, coorid, pv, dpv, ierr)
      enddo
      write(*,*) dtime(tarray)
      end