Added rouitines in plot.subhalos to find subhalo functions from given halo objects. The new added routines are heavy on memory usage so the old ones have been left there as well