We consider a two steps solution procedure composed by: a first step where the non hydrostatic source term is recovered by inverting the elliptic coercive operator associated to the dispersive effects; a second step which involves the solution of the hyperbolic shallow water system with the source term, computed in the previous phase, which accounts for the non-hydrostatic effects. Appropriate numerical methods that can be also generalized on arbitrary unstructured meshes are used to discretize the two stages: the standard $C^0$ Galerkin finite element method for the elliptic phase; either third order Finite Volume of third order stabilized Finite Element methods for the hyperbolic phase. The discrete dispersion properties of the fully coupled schemes obtained are studied, showing accuracy close or better to that of a fourth order finite difference method. The hybrid approach of locally reverting to the nonlinear shallow water equations is used to recover energy dissipation in breaking regions. To this scope we evaluate two strategies~: simply neglecting the non-hydrostatic contribution in the hyperbolic phase~; imposing a tighter coupling of the two phases, with a wave breaking indicator embedded in the elliptic phase to smoothly turn off the dispersive effects. The discrete models obtained are thoroughly tested on benchmarks involving wave dispersion, breaking and runup, showing a very promising potential for the simulation of complex near shore wave physics in terms of accuracy and robustness.< Réduire