program xephem_2 ! ! Exercise driver for subroutine ephem--call it many times ! in a loop (for the Sun) to test performance. ! ! B. Knapp, 2002-03-15 ! implicit none integer*4 i, n parameter (n=10000) 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=1) !Spherical coordinates real*8 jd, jd0, pv(6) parameter (jd0=2447892.5d0) !1990 Jan 1.0 integer*4 ierr ! Externals external ephem do i=0,n-1 jd = jd0 + i call ephem(jd, p, posnid, framid, coorid, pv, ierr) enddo end