General purpose data structures for the motion of completely faceted, two-dimensional surfaces as described in [2] were used. To adapt that program to this problem involved changing the formula by which the velocity of a facet is computed. These crystalline methods work for initial data consistent with definitions in section 2.1 and can automatically handle topological changes in the surface.