include bhs.h routine to obtaine Mdot and distance of a black hole binary