now read in extant phylos instead of trying to recreate them
adding in sim4065 @ time=30k
removed simID argument from bamm.plot; add 4-panel bamm.plot