>>>  vr = halo[i].virial_radius(virial_overdensity=200.0, bins=100)

will return the virial radius of the halo in code units.

What is code units? How to convert it into kps or mps?

--