diff --git a/gmsh/cube_tetra3d4.geo b/gmsh/cube_tetra3d4.geo
index d7880e5..3cb109b 100644
--- a/gmsh/cube_tetra3d4.geo
+++ b/gmsh/cube_tetra3d4.geo
@@ -1,76 +1,102 @@
 //+
-Point(1) = {0, 0, 0, 0.1};
+Point(1) = {0, 0, 0, 0.05};
 //+
-Point(2) = {1, 0, 0, 0.1};
+Point(2) = {1, 0, 0, 0.05};
 //+
-Point(3) = {0, 1, 0, 0.1};
+Point(3) = {1, 1, 0, 0.05};
 //+
-Point(4) = {1, 1, 0, 0.1};
+Point(4) = {0, 1, 0, 0.05};
+//+
+Point(5) = {0, 0, 1, 0.05};
+//+
+Point(6) = {1, 0, 1, 0.03};
+//+
+Point(7) = {1, 1, 1, 0.01};
+//+
+Point(8) = {0, 1, 1, 0.03};
 //+
 Line(1) = {1, 2};
 //+
-Line(2) = {2, 4};
+Line(2) = {2, 3};
 //+
-Line(3) = {4, 3};
+Line(3) = {3, 4};
 //+
-Line(4) = {3, 1};
+Line(4) = {4, 1};
 //+
-Curve Loop(1) = {4, 1, 2, 3};
+Line(5) = {1, 5};
 //+
-Plane Surface(1) = {1};
+Line(6) = {2, 6};
 //+
-Extrude {0, 0, 1} {
-  Surface{1}; 
-}
+Line(7) = {3, 7};
 //+
-Delete {
-  Surface{26}; 
-}
+Line(8) = {4, 8};
 //+
-Delete {
-  Surface{26}; 
-}
+Point(9) = {0.5, 0.5, 1, 0.02};
 //+
-Delete {
-  Volume{1}; 
-}
+Point(10) = {1.0, 0.5, 1, 0.02};
 //+
-Delete {
-  Surface{26}; 
-}
+Point(11) = {0.5, 1.0, 1, 0.02};
 //+
-Point(15) = {0.5, 0.5, 1.0, 0.1};
+Line(9) = {5, 6};
 //+
-Line(21) = {15, 10};
+Line(10) = {6, 10};
 //+
-Line(22) = {15, 14};
+Line(11) = {10, 7};
 //+
-Line(23) = {15, 5};
+Line(12) = {7, 11};
 //+
-Line(24) = {15, 6};
+Line(13) = {11, 8};
 //+
-Curve Loop(2) = {7, -21, 24};
+Line(14) = {8, 5};
 //+
-Surface(26) = {2};
+Line(15) = {9, 10};
 //+
-Curve Loop(3) = {21, 8, -22};
+Line(16) = {9, 11};
 //+
-Surface(27) = {3};
+Curve Loop(1) = {3, 4, 1, 2};
 //+
-Curve Loop(4) = {22, 9, -23};
+Surface(1) = {1};
 //+
-Surface(28) = {4};
+Curve Loop(2) = {4, 5, -14, -8};
 //+
-Curve Loop(5) = {23, 6, -24};
+Plane Surface(2) = {2};
 //+
-Surface(29) = {5};
+Curve Loop(3) = {8, -13, -12, -7, 3};
 //+
-Surface Loop(1) = {26, 17, 1, 13, 29, 28, 27, 21, 25};
+Plane Surface(3) = {3};
+//+
+Curve Loop(4) = {9, -6, -1, 5};
+//+
+Plane Surface(4) = {4};
+//+
+Curve Loop(5) = {6, 10, 11, -7, -2};
+//+
+Plane Surface(5) = {5};
+//+
+Curve Loop(6) = {15, 11, 12, -16};
+//+
+Plane Surface(6) = {6};
+//+
+Curve Loop(7) = {13, 14, 9, 10, -15, 16};
+//+
+Plane Surface(7) = {7};
+//+
+Surface Loop(1) = {4, 7, 3, 2, 1, 5, 6};
 //+
 Volume(1) = {1};
 //+
-Physical Surface("top", 30) = {27};
+Physical Volume("cube", 17) = {1};
 //+
-Physical Surface("bottom", 31) = {1};
+Physical Surface("bottom", 18) = {1};
 //+
-Physical Volume("body", 32) = {1};
+Physical Surface("top_big", 19) = {7};
+//+
+Physical Surface("top_small", 20) = {6};
+//+
+Physical Curve("top_south", 21) = {9};
+//+
+Physical Curve("top_west", 22) = {14};
+//+
+Physical Curve("top_east", 23) = {10, 11};
+//+
+Physical Curve("top_north", 24) = {12, 13};