richness_change function

Calculate change in richness resulting from a percent reduction in flow