Mesh.MshFileVersion = 1; 
L =   10.00000;
B =   10.00000;
H =   10.00000;
b =    9.00000;
c =    2.50000;
Fin =    3.00000;
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
////////// 







Mesh.ElementOrder = 2;
//+
Point(1) = {0, 0, 0, Fin};
//+
Point(2) = {L, 0, 0, Fin};
//+
Point(3) = {L, B, 0, Fin};
//+
Point(4) = {0, B, 0, Fin};
//+
Point(5) = {(L-c)/2, (L-b)/2, 0, Fin};
//+
Point(6) = {(L-c)/2+c, (L-b)/2, 0, Fin};
//+
Point(7) = {(L-c)/2+c, (B-c)/2, 0, Fin};
//+
Point(8) = {L-(L-b)/2, (B-c)/2, 0, Fin};
//+
Point(9) = {L-(L-b)/2, (B-c)/2+c, 0, Fin};
//+
Point(10) = {(L-c)/2+c, (B-c)/2+c, 0, Fin};
//+
Point(11) = {(L-c)/2+c, B-(B-b)/2, 0, Fin};
//+
Point(12) = {(L-c)/2, B-(B-b)/2, 0, Fin};
//+
Point(13) = {(L-c)/2, (B-c)/2+c, 0, Fin};
//+
Point(14) = {(L-b)/2, (B-c)/2+c, 0, 1.0};
//+
Point(15) = {(L-b)/2, (B-c)/2, 0, 1.0};
//+
Point(16) = {(L-c)/2, (B-c)/2, 0, 1.0};
//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {3, 4};
//+
Line(4) = {4, 1};
//+
Line(5) = {5, 6};
//+
Line(6) = {6, 7};
//+
Line(7) = {7, 8};
//+
Line(8) = {8, 9};
//+
Line(9) = {9, 10};
//+
Line(10) = {10, 11};
//+
Line(11) = {11, 12};
//+
Line(12) = {12, 13};
//+
Line(13) = {13, 14};
//+
Line(14) = {14, 15};
//+
Line(15) = {15, 16};
//+
Line(16) = {16, 5};
//+
Line Loop(1) = {1, 2, 3, 4};
//+
Line Loop(2) = {7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 5, 6};
//+
Plane Surface(1) = {1, 2};

out[] = Extrude {0, 0, 10} { Surface{1};};

//For t In {0:17}
//Printf("%g",out[t]);
//EndFor

DOMAIN_1=1000;
Physical Volume(DOMAIN_1)={out[1]};

BOUNDARY_1001 = 1001;
BOUNDARY_1002 = 1002;
BOUNDARY_1003 = 1003; 
BOUNDARY_1004 = 1004;
BOUNDARY_1005 = 1005;
BOUNDARY_1006 = 1006;
BOUNDARY_1007 = 1007;

Physical Surface(BOUNDARY_1001) = {-1};
Physical Surface(BOUNDARY_1002) = {out[0]};
Physical Surface(BOUNDARY_1003) = {out[2]};
Physical Surface(BOUNDARY_1004) = {out[3]};
Physical Surface(BOUNDARY_1005) = {out[4]};
Physical Surface(BOUNDARY_1006) = {out[5]};
Physical Surface(BOUNDARY_1007) = {out[6],out[7],out[8],out[9],out[10],out[11],out[12],out[13],out[14],out[15],out[16],out[17]};

Periodic Surface{out[3]} = {out[5]} Translate{L,0,0};
Periodic Surface{out[4]} = {out[2]} Translate{0,B,0};