A new two-phase flow model based on coupling of the depth-integrated continuum method and discrete element method. (January 2021)