integer*8 peph
integer res
double precision AU

res = f90calceph_open(peph, "example1.dat")
if (res.eq.1) then
  ! print the time scale
  timescale = calceph_gettimescale(peph)
  write (*,*) "timescale=", timescale

  call f90calceph_close(peph)
endif