def METHOD_NAME(self):
    new_corner_lons = numpy.array([-1.1, 1.1, -1.1, 1.1])
    new_corner_lats = numpy.array([0., 0., 0., 0.])
    new_corner_depths = numpy.array([0., 0., 15., 15.])
    points = [Point(lon, lat, depth) for lon, lat, depth in
              zip(new_corner_lons, new_corner_lats, new_corner_depths)]
    new_surface = PlanarSurface(self.STRIKE, self.DIP,
                                points[0], points[1], points[3], points[2])
    self.fault.modify_set_geometry(new_surface, "XYZ")
    self.assertEqual(self.fault.surface_node, "XYZ")
    numpy.testing.assert_array_almost_equal(self.fault.surface.corner_lons,
                                            new_corner_lons)
    numpy.testing.assert_array_almost_equal(self.fault.surface.corner_lats,
                                            new_corner_lats)
    numpy.testing.assert_array_almost_equal(
        self.fault.surface.corner_depths, new_corner_depths)