A numerical method, based on the Zakharov equation in three dimensions, is developed to study the nonlinear dynamics of deep water gravity waves.We focus primarily on systems describing the evolution of hydrodynamic surface waves, considering a thick incompressible fluid spanning across all space available and subject to gravity. The reference height to the resting state is z = 0. In the case of hydrostatic equilibrium, the bottom of the water is located at a distance H from the water surface. This depth is considered constant. The height of the surface wave is denotedx,y,t) and the volume occupied by the fluid is restricted by –H < z