cl1 = 1;
h = 3;
r = 5;

Point(1) = {0, 0, 0, cl1};
Point(2) = {r, 0, 0, cl1};
Point(3) = {0, r, 0, cl1};
Point(4) = {-r, 0, 0, cl1};
Point(5) = {0, -r, 0, cl1};
Circle(1) = {2, 1, 3};
Circle(2) = {3, 1, 4};
Circle(3) = {4, 1, 5};
Circle(4) = {5, 1, 2};
Point(11) = {0, 0, h, cl1};
Point(12) = {r, 0, h, cl1};
Point(13) = {0, r, h, cl1};
Point(14) = {-r, 0, h, cl1};
Point(15) = {0, -r, h, cl1};
Circle(11) = {12, 11, 13};
Circle(12) = {13, 11, 14};
Circle(13) = {14, 11, 15};
Circle(14) = {15, 11, 12};
Line(15) = {14, 4};
Line(16) = {15, 5};
Line(17) = {12, 2};
Line(18) = {13, 3};
Line Loop(19) = {3, 4, 1, 2};
Plane Surface(20) = {19};
Line Loop(21) = {13, 14, 11, 12};
Plane Surface(22) = {21};
Line Loop(23) = {15, 3, -16, -13};
Ruled Surface(24) = {23};
Line Loop(25) = {4, -17, -14, 16};
Ruled Surface(26) = {25};
Line Loop(27) = {1, -18, -11, 17};
Ruled Surface(28) = {27};
Line Loop(29) = {12, 15, -2, -18};
Ruled Surface(30) = {29};
Surface Loop(31) = {26, 20, 24, 30, 22, 28};
Volume(32) = {31};

Mesh.Algorithm = 7; //Bamg
Mesh.Algorithm3D= 7; //mmg3d
Mesh.CharacteristicLengthExtendFromBoundary = 0;


Field[1] = AttractorAnisoCurve;
Field[1].EdgesList = {1,2,3,4};
Field[1].NNodesByEdge = 50;
Field[1].dMax = h*0.25;
Field[1].dMin = h*0.1;
Field[1].lMaxNormal = h / 5;
Field[1].lMaxTangent = r / 10;
Field[1].lMinNormal = h / 50;
Field[1].lMinTangent = r / 50;

Background Field = 1;
