BUGFIX: compute omega_mean only in range cg%is:cg%ie range

Fixes ill defined mean on internal boundaries