fix coordinate system
added function to compute Read tidal radius
don't always compute tail_bit