let mod2pi t =
  mod_float t twopi