My branch of Stefany's code to test whether sims fit with analytical stuff