integer res
integer timescale
TYPE(C_PTR) :: peph

peph = calceph_open("example1.dat"//C_NULL_CHAR)
if (C_ASSOCIATED(peph)) then
   ! print the time scale
   timescale = calceph_gettimescale(peph)
   write (*,*) "timescale=", timescale

   call calceph_close(peph)
 endif