Mesh.MshFileVersion = 1;
R1 =    0.00000;
R2 =    0.58400;
H =    0.26300;
Fin =   16.00000;
h = 0.02;
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//////////


// domaine 1
Point(1) = {R1,0,-h/2,Fin};
Point(2) = {R2/58,0,-h/2,Fin};
Point(3) = {R2/58,H/26,-h/2,Fin};
Point(4) = {R1,H/26,-h/2,Fin};

Point(5) = {R1,0,h/2,Fin};
Point(6) = {R2/58,0,h/2,Fin};
Point(7) = {R2/58,H/26,h/2,Fin};
Point(8) = {R1,H/26,h/2,Fin};


// contour
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,5};

Line(9) = {1,5};
Line(10) = {2,6};
Line(11) = {3,7};
Line(12) = {4,8};


Line Loop(9) = {1,2,3,4};
Plane Surface(10) = {9};

Line Loop(11) = {5,6,7,8};
Plane Surface(12) = {11};

// s t r u c t u r a t i o n du m aill a g e
Transfinite Line {1,3,5,7}=4;
Transfinite Line {2,4,6,8}=3;
Transfinite Surface "*";

//maillage quadrangulaire
Recombine Surface "*";
Transfinite Line {9,10,11,12}=4;

Geometry.CopyMeshingMethod = 1;
Translate {R2/58*1,H/26*0,0}{Duplicata { Line{9,10,11,12}; Surface{10,12};}}
Translate {R2/58*2,H/26*0,0}{Duplicata { Line{9,10,11,12}; Surface{10,12};}}
DOMAIN_1=1000;

// mesh everything at order 1
Mesh 2;

// display only surfaces, and set the mesh order of visible entities to 2
Mesh.MeshOnlyVisible = 1;
Mesh.ElementOrder = 2;
Mesh.SecondOrderIncomplete = 1;
Hide "*";
Show { Surface { Surface "*" } ; }
SetOrder 2;
Mesh.MeshOnlyVisible = 0;

// show everything again
Show "*";
