
h = 1;
Mesh.ElementOrder = 3;


Point(1) = {0, 0, 0, h};
Point(2) = {1, 0, 0, h};
Point(3) = {1,1, 0, h};
Point(4) = {0, 1, 0, h};
Point(5) = {2, 0, 0, h};
Point(6) = {2, 1, 0, h};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
Line(5) = {2, 5};
Line(6) = {5, 6};
Line(7) = {6, 3};

// Physical domain
Physical Surface(1000) = {8};
Physical Surface(1001) = {9};
Physical Line(1003) = {1, 2, 3, 4};
Physical Line(1004) = {5, 6, 7, -2};
Line Loop(1005) = {4, 1, 2, 3};
Plane Surface(1006) = {1005};
Line Loop(1007) = {-2, 7, 6, 5};
Plane Surface(1008) = {1007};
Physical Surface(1009) = {1006};
Physical Surface(1010) = {1008};
