integer*8 peph
integer res
double precision jd0
double precision dt1
double precision P(3)

jd0 = 2442457
dt1 = 0.5D0
res = f90calceph_open(peph, "example1.dat")
if (res.eq.1) then

  res = f90calceph_orient_order(peph,jd0, dt1, NAIFID_MOON,
&                              CALCEPH_USE_NAIFID+CALCEPH_UNIT_RAD+CALCEPH_UNIT_SEC,
&                              0, P)
  write(*,*) P

  call f90calceph_close(peph)
endif