Oyster_phase_run function

Runs the entire trace element model