phint_length function

Compute the length of a phinterval in seconds