Oyster_subincr_phases function

Calculates proportion of phases in each sub-increment