update relatedness normalization
Add an option -r to choose the time window to calculate relatedness. Default -1 means do not calculate relatedness.
Move systematic calculation of relatedness in a separate function.
Output relatedness in a file instead of stdout.
Avoid division by 0 when population is invaded from a single ancestor
(relatedness undefined, we take 0).