The LBM method is a relatively recent method. Contrary to the classic methods of finite elements, finite differences and finite volumes which are based on the macroscopic description of the variables, LBM is based on statistical physics to propose a mesoscopic description [1]. The LBM algorithm is of a massively parallelizable nature, especially with GPU technology. This parallelization property makes LBM more efficient, more attractive and the most preferred of classical methods. Several researchers around the world have applied LBM to their field of application (subterranean flows, hydrodynamics, aerodynamics, multiphysics flows, acoustics, and many others). Similarly LBM has been applied to coastal flows governed by 2D shallow water type equations [2]. The aim of this work is to extend LBM to 3D coastal flows in order to be able to simulate the characteristics of the flow according to the third dimension. The knowledge of the flow according to the vertical is imposed by the phenomena of the turbulence inseparable from the coastal flows. The 3D version that we propose will also allow simulations of other closely related phenomena