Four Padé–Petrov–Galerkin models, for numerical computation discontinuous flow in open channels, are proposed in this work. Time discretization, in multi-stage approach of Saint-Venant’s basic equations, is obtained by first- and second-sub-diagonal L-stable Padé approximations in...