mrgsim_q function

Simulate from a model object with quicker turnaround