diff --git a/LDDsimulation/domainSubstructuring.py b/LDDsimulation/domainSubstructuring.py
index 4048ef4957c74f0925c383a70008d4150fd5cbf7..bbedaf8ef9fe465f3645fae7fc09d6a50862820a 100644
--- a/LDDsimulation/domainSubstructuring.py
+++ b/LDDsimulation/domainSubstructuring.py
@@ -79,14 +79,14 @@ class twoSoilLayers(domainSubstructuring):
         self.__subdomain1_vertices = [
             self.__interface12_vertices[0],
             self.__interface12_vertices[1],
-            self.__sub_domain0_vertices[2],
-            self.__sub_domain0_vertices[3]
+            self.__subdomain0_vertices[2],
+            self.__subdomain0_vertices[3]
             ]
 
         # subdomain2
         self.__subdomain2_vertices = [
-            self.__sub_domain0_vertices[0],
-            self.__sub_domain0_vertices[1],
+            self.__subdomain0_vertices[0],
+            self.__subdomain0_vertices[1],
             self.__interface12_vertices[1],
             self.__interface12_vertices[0]]
 
@@ -105,15 +105,15 @@ class twoSoilLayers(domainSubstructuring):
         # dictionary entry should be 0: None
         self.__subdomain1_outer_boundary_verts = {
             0: [self.__interface12_vertices[1],
-                self.__sub_domain0_vertices[2],
-                self.__sub_domain0_vertices[3],
+                self.__subdomain0_vertices[2],
+                self.__subdomain0_vertices[3],
                 self.__interface12_vertices[0]]
         }
 
         self.__subdomain2_outer_boundary_verts = {
             0: [self.__interface12_vertices[0],
-                self.__sub_domain0_vertices[0],
-                self.__sub_domain0_vertices[1],
+                self.__subdomain0_vertices[0],
+                self.__subdomain0_vertices[1],
                 self.__interface12_vertices[1]]
         }