PCS_run function

Simulate the run of a PCS model