def METHOD_NAME(self, faces=None, color=None, join_faces=False):
    """Draw a selection of faces.
    Parameters
    ----------
    faces : list[list[int]], optional
        A list of faces to draw.
        The default is None, in which case all faces are drawn.
    color : :class:`~compas.colors.Color` | dict[int, :class:`~compas.colors.Color`], optional
        The color specification for the faces.
        The default color is :attr:`VolMeshArtist.default_facecolor`.
    join_faces : bool, optional
        If True, join the faces into one mesh.
    Returns
    -------
    list[:rhino:`Rhino.Geometry.Mesh`]
    """
    self.face_color = color
    faces = faces or self.faces
    vertex_xyz = self.vertex_xyz
    facets = []
    for face in faces:
        facets.append(
            {
                "points": [vertex_xyz[vertex] for vertex in self.volmesh.halfface_vertices(face)],
                "name": "{}.face.{}".format(self.volmesh.name, face),
                "color": self.face_color[face].rgb255,
            }
        )
    return compas_ghpython.METHOD_NAME(facets)