The role of bed-load thickness closure in a two-phase morphodynamical model