As already noted, the process of embedding a discontinuous surface in a matrix of computational cells involves three separate tasks. First, it is necessary to devise a means of numerically describing the location and shape of the boundary. Second, an algorithm must be given for computing the time evolution of the boundary