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