integer*8 peph
integer res
double precision AU

res = f90calceph_open(peph, "example1.dat")
if (res.eq.1) then
   ! print the value of AU
   if (f90calceph_getconstant(peph, "AU", AU).eq.1) then
        write (*,*) "AU=", AU
    endif

  call f90calceph_close(peph)
endif