def read(self, debug, interpolate):
    """
    Read finite-element mesh and store in Sieve mesh object.

    @returns PETSc mesh object containing finite-element mesh
    """
    from pylith.mpi.Communicator import mpi_comm_world
    comm = mpi_comm_world()
    if 0 == comm.rank:
      self._info.log("Reading finite-element mesh")

    # Set flags
    self.debug(debug)
    self.interpolate(interpolate)

    # Initialize coordinate system
    if self.coordsys is None:
      raise ValueError, "Coordinate system for mesh is unknown."

    from pylith.mpi.Communicator import petsc_comm_world
    from pylith.topology.Mesh import Mesh    
    mesh = Mesh(dim=self.coordsys.spaceDim(), comm=petsc_comm_world())
    mesh.coordsys(self.coordsys)

    # Read mesh
    ModuleMeshIO.read(self, mesh)
    return mesh

 def test_dimension(self):
   """
   Test debug().
   """
   mesh = Mesh()
   self.assertEqual(0, mesh.dimension())
   return

 def test_coordsys(self):
   """
   Test coordsys().
   """
   mesh = self._getMesh()
   submesh = Mesh(mesh=mesh, label="bc")
   self.assertEqual(2, submesh.coordsys().spaceDim())
   return

 def test_constructorA(self):
   """
   Test constructor.
   """
   mesh = Mesh()
   self.assertEqual(0, mesh.dimension())
   self.assertEqual(False, mesh.debug())
   return

  def test_dimension(self):
    """
    Test debug().
    """
    mesh = self._getMesh()
    submesh = Mesh(mesh=mesh, label="bc")

    self.assertEqual(1, submesh.dimension())
    return

 def test_constructorB(self):
   """
   Test constructor.
   """
   mesh = self._getMesh()
   submesh = Mesh(mesh=mesh, label="bc")
   self.assertEqual(1, submesh.dimension())
   self.assertEqual(False, mesh.debug())
   return

 def test_constructorB(self):
   """
   Test constructor.
   """
   dim = 2
   mesh = Mesh(dim=dim)
   self.assertEqual(dim, mesh.dimension())
   self.assertEqual(False, mesh.debug())
   return

 def test_constructorC(self):
   """
   Test constructor.
   """
   dim = 2
   from pylith.mpi.Communicator import mpi_comm_self
   mesh = Mesh(dim=dim, comm=mpi_comm_self())
   self.assertEqual(dim, mesh.dimension())
   self.assertEqual(False, mesh.debug())
   return

 def test_comm(self):
   """
   Test comm().
   """
   from pylith.mpi.Communicator import petsc_comm_self
   mesh = Mesh(dim=3, comm=petsc_comm_self())
   comm = mesh.comm()
   self.assertEqual(0, comm.rank)
   self.assertEqual(1, comm.size)
   return

  def test_comm(self):
    """
    Test comm().
    """
    mesh = self._getMesh()
    submesh = Mesh(mesh=mesh, label="bc")

    comm = submesh.comm()
    self.assertEqual(0, comm.rank)
    self.assertEqual(1, comm.size)
    return

  def test_debug(self):
    """
    Test debug().
    """
    mesh = Mesh()

    self.assertEqual(False, mesh.debug())

    mesh.debug(True)
    self.assertEqual(True, mesh.debug())
    return

  def test_debug(self):
    """
    Test debug().
    """
    mesh = self._getMesh()
    submesh = Mesh(mesh=mesh, label="bc")

    self.assertEqual(False, submesh.debug())

    submesh.debug(True)
    self.assertEqual(True, submesh.debug())
    return

  def test_coordsys(self):
    """
    Test coordsys().
    """
    from spatialdata.geocoords.CSCart import CSCart
    cs = CSCart()
    cs.inventory.spaceDim = 2
    cs._configure()

    mesh = Mesh()
    mesh.coordsys(cs)
    self.assertEqual(cs.spaceDim(), mesh.coordsys().spaceDim())
    return