Added set_population() function.
The simulate_* functions can now handle population sizes that vary across regions.