split into several smaller functions
Add a function that takes two generations (gen0 and gen1) and tag each individual at gen1 with this ancestor at generation gen0. We should be able to use this to compute relatedness because it is giving us identity by descent
fix a bug: report[gen][i] is actually giving replication report creating individual 'i' at generation 'gen+1' (and not 'gen' as it was supposed in calculation of reproductive success)