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