sim_distribution function

Simulate spatial and temporal distribution