As is discused on this post, taking some assumptions, the water surface can be simulated by a discrete aproximation of a grid of heights using this formula
Where:
HT is the new height grid
HT-1 is the grid mapped in the previous step of the simulation
N is a damping factor
The 2D simulation works pretty well. I want now to move to 3D, but I'm concerned about the bundary limit conditions.
Unless I'm missing something, the document doesn't explain how to calculate the height on the limits, where there is no x+1 or x-1.
I could mirror the point on the other site as if the boundaries were mirrors, but I'm not sure.
The following image depicts what's described above