integer*8 peph
integer res
integer n

res = f90calceph_open(peph, "example1.dat")
if (res.eq.1) then
  n = f90calceph_getpositionrecordcount(peph)
  write (*,*) "number of position's record", n
  call f90calceph_close(peph)
endif