modes function

Identify multiple modes