Modified to calculate mean and standard deviations of each metric across simulation replicates. Plot now includes envelope +/- 1 SD around mean for just 4 metrics: global richness, latitude-richness corr, gamma, and MRD-richness corr