  We consider here different family of Boussinesq systems in two space
dimensions. These systems approximate the three-dimensional Euler equations and
consist of three coupled nonlinear dispersive wave equations that describe
propagation of long surface waves of small amplitude in ideal fluids over a
horizontal bottom. We present here a FreeFem++ code aimed at solving
numerically these systems where a discretization using P1 finite element for
these systems was taken in space and a second order Runge-Kutta scheme in time.
We give the detail of our code where we use a mesh adaptation technique. An
optimization of the used algorithm is done and a comparison of the solution for
different Boussinesq family is done too. The results we obtained agree with
those of the literature.
