Extension of the gambin model to multimodal species abundance distributions