diff --git a/src/bvpy/domains/custom_domain.py b/src/bvpy/domains/custom_domain.py
index 553a49fda9b1e2475c0a2f37814569c7f5167d45..a3a0ae0438b608230d2591c07e8ae09175fca666 100644
--- a/src/bvpy/domains/custom_domain.py
+++ b/src/bvpy/domains/custom_domain.py
@@ -83,14 +83,12 @@ class CustomDomain(AbstractDomain, BuiltInModel):
 
         points = np.array(points)
         cells = np.array(cells)
-        # self.mesh = fe.Mesh()
-        self.mesh = None
 
         self.dim = points.shape[1]
 
-        print(f'self.dim = {self.dim}')
-
         # editor = fe.MeshEditor()
+
+        # OA: What is the purpose of the following block ?
         if self._cell_type == 'line':
             topo_dim = 1 # OA: changed variable name top -> topo_dim
                          # OA: what is the use of this variable ?
@@ -98,6 +96,7 @@ class CustomDomain(AbstractDomain, BuiltInModel):
             topo_dim = 2
         elif self._cell_type == 'tetra':
             topo_dim = 3
+            
         # editor.open(self.mesh, self._cell_type, topo_dim, self.dim)
         # editor.init_vertices(len(points))
         # editor.init_cells(len(cells))
@@ -107,7 +106,6 @@ class CustomDomain(AbstractDomain, BuiltInModel):
 
         # editor.close()
 
-        print(f'pts used to create mesh = {points[:, :self.dim]}')
 
         _cell = ufl.Cell(self._cell_type, geometric_dimension=self.dim)
         _domain = ufl.Mesh(ufl.VectorElement("CG", _cell, 1))