rectangular mesh of cells of width ox and height Oy one might define the vertical height, h, of the free boundary above the bottom of the mesh in each column of cells. This would approximate a curve h = f(x, t) by assigning values of h to discrete values of x. This method does not work well when the boundary slope, dh/dx,