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