Completed PD calculation code for colony subset.

It is working but I would like to figure out how
to make it wait for processes to end, and then
continue with the rest of the analysis.