Running a simple simulation with quasi-steady model approach (input.mdv.flowtype=2) : constant discharge upstream, constant water level downstream, still doesn't work. It gives following Warning: the power of the excess shield should be the same in the sediment transport relation than in the entrainement formulation. Otherwise tt may create a discontinuity in Dk because it is 0 when Qbk_st is 0 but it may cancel with the term in Ek_st 
And ends with ¡¡ ERROR !! 
Supposedly you do not end up here in bed_level_update, at line 137 
Same simulation with input.mdv.flowtype=4 (onlu change in input) does start, although it gives the same warning. Any idea/tip?