Add new functions to shared_ancestry_simulator.R; can run single old model as was, now with max_cons output