Allows options of running the simulations in parallel on multi-core machines
Added option to generate density grid from a gam object.
Added option to generate cluster sizes for observations.
Note: generating a population based on D (rather than N) does not work correctly.