This code takes a clade ID as identified from the extant phylogeny (sub.phylo) at the end of the simulation, and identifies all of the clade members both extinct and extant. It is not fully automated (requires some inspection by the investigator), but uses the following steps. 1) Find the node number in phylo.out (which includes extinct species) that corresponds to the clade ID in sub.phylo. We do this by identifying all of the extant descendants, finding the pair with the greatest phylogenetic distance between them, and then calculating the most recent common ancestor of those species within phylo.out. Seems like we may potentially miss species that went extinct that were phylogenetically peripheral, and therefore the true node representing the clade might be a split or two deeper, but this can't be that common and should only be a difference of an extinct species or two.