Halo.
bulk_velocity
Returns the mass-weighted average velocity in cm/s.
This calculates and returns the mass-weighted average velocity of just the particles in the halo in cm/s.
Examples
>>> bv = halos[0].bulk_velocity()