program test_ph
implicit none
! variables needed for ephemeris
integer, parameter 			:: ncmax=100
integer, parameter 			:: nimax=3000
integer 	   			:: nc,ni
real(kind=8)       			:: a,b
real(kind=8), dimension(ncmax,nimax,6)	:: cf
real(kind=8), dimension(nimax,2) 	:: dat
real(kind=8)				:: t
real(kind=8), dimension(6)		:: x
!
integer :: i
!
call header(nc,ni,a,b)		! Read the header
call load(nc,ni,cf,dat)		! Load Chebychev coefficients
do i=0,10
	t=2450000.5d0+i*10d0
	call pos_PH(a,b,nc,ni,cf,dat,t,x)
	write(*,20) 'DJ = ',t
	write(*,21) 'x  = ',x(1:3)
	write(*,21) 'v  = ',x(4:6)
	write(*,*) ' '
enddo
20 format(a5,f12.3) 
21 format(a5,3f20.13) 
end program
include 'r-ph.f90'