The boundary element method (BEM) with Lagrangian formulation is a conceptually simple and efficient for the simulation of nonlinear wave shoaling, or without impermeable coastal structures, up to breaking. However, in post-breaking flows, domain no longer simply connected, BEM not generation new free surface. Volumes fluid (VOF) methods are made track fluid-free surfaces after breaking, but th...