Hello,
Actually the Boussinesq equations are Saint-Venant equations with extra terms in the momentum to take into account vertical accelerations. The celerity of long waves is sqrt(gravity * depth) and it is what give the Saint-Venant equations, it is not exact for non linear waves. In our application we looked at the waves arriving at a dam. There were two waves with Saint-Venant equations and only one with Boussinesq equations.
In terms of wave run-up there is no difference, there is a reference depth with Boussinesq, this reference depth is zero on dry land, so the extra Boussinesq terms vanish.
In Telemac-2D the Boussinesq equations (you just need the keyword EQUATIONS : BOUSSINESQ in your parameter file) can only be solved with primitive equations : TREATMENT OF THE LINEAR SYSTEM : 1, with SOLVER : 7 (GMRES), so it will be more expensive than Saint-Venant with wave equation. It could even be that Telemac-3D with 2 or 3 planes is less expensive. Note that Telemac-3D with non hydrostatic option and two planes has a vertical velocity linear between bottom and free surface, which is the starting assumption of Serre equations, another "Boussinesq like" set of equations.
You can set the wave run-up zones as tidal flat area.
With best regards, have a good week-end,
Jean-Michel Hervouet