// Radius
r=0.5;
// Height
h=1.5;
// Length
l=14*r;
// Spacing
s=1.1;

// Lower body
Point(1) = {0, 0, 0, 1.0};
Point(2) = {l, 0, 0, 1.0};
Point(3) = {l, h, 0, 1.0};
Point(4) = {0, h, 0, 1.0};

// Upper body
Point(5) = {0, h+s, 0, 1.0};
Point(6) = {l, h+s, 0, 1.0};
Point(7) = {l, 2*h+s, 0, 1.0};
Point(8) = {0, 2*h+s, 0, 1.0};

// Midpoints lower body
Point(9) = {r, h, 0, 1.0};
Point(10) = {2*r, h, 0, 1.0};
Point(11) = {3*r, h, 0, 1.0};
Point(12) = {4*r, h, 0, 1.0};
Point(13) = {5*r, h, 0, 1.0};
Point(14) = {6*r, h, 0, 1.0};
Point(15) = {7*r, h, 0, 1.0};
Point(16) = {8*r, h, 0, 1.0};
Point(17) = {9*r, h, 0, 1.0};
Point(18) = {10*r, h, 0, 1.0};
Point(19) = {11*r, h, 0, 1.0};
Point(20) = {12*r, h, 0, 1.0};
Point(21) = {13*r, h, 0, 1.0};

// Midpoints upper body
Point(22) = {r, h+s, 0, 1.0};
Point(23) = {2*r, h+s, 0, 1.0};
Point(24) = {3*r, h+s, 0, 1.0};
Point(25) = {4*r, h+s, 0, 1.0};
Point(26) = {5*r, h+s, 0, 1.0};
Point(27) = {6*r, h+s, 0, 1.0};
Point(28) = {7*r, h+s, 0, 1.0};
Point(29) = {8*r, h+s, 0, 1.0};
Point(30) = {9*r, h+s, 0, 1.0};
Point(31) = {10*r, h+s, 0, 1.0};
Point(32) = {11*r, h+s, 0, 1.0};
Point(33) = {12*r, h+s, 0, 1.0};
Point(34) = {13*r, h+s, 0, 1.0};//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {6, 7};
//+
Line(4) = {7, 8};
//+
Line(5) = {8, 5};
//+
Line(6) = {4, 1};
//+
Circle(7) = {10, 9, 4};
//+
Circle(8) = {14, 13, 12};
//+
Circle(9) = {18, 17, 16};
//+
Circle(10) = {3, 21, 20};
//+
Circle(11) = {10, 11, 12};
//+
Circle(12) = {14, 15, 16};
//+
Circle(13) = {18, 19, 20};
//+
Circle(14) = {5, 22, 23};
//+
Circle(15) = {25, 26, 27};
//+
Circle(16) = {29, 30, 31};
//+
Circle(17) = {33, 34, 6};
//+
Circle(18) = {33, 32, 31};
//+
Circle(19) = {29, 28, 27};
//+
Circle(20) = {25, 24, 23};
//+
Curve Loop(1) = {1, 2, 10, -13, 9, -12, 8, -11, 7, 6};
//+
Plane Surface(1) = {1};
//+
Curve Loop(2) = {19, -15, 20, -14, -5, -4, -3, -17, 18, -16};
//+
Plane Surface(2) = {2};
//+
Physical Surface("lower_body", 21) = {1};
//+
Physical Surface("upper_body", 22) = {2};
//+
Physical Curve("top_edge_lower_body", 23) = {7, 11, 8, 12, 9, 13, 10};
//+
Physical Curve("bottom_edge_upper_body", 24) = {14, 20, 15, 19, 16, 18, 17};
//+
Physical Curve("bottom_edge_lower_body", 25) = {1};
//+
Physical Curve("right_edge_lower_body", 26) = {2};
//+
Physical Curve("left_edge_lower_body", 27) = {6};
//+
Physical Curve("top_edge_upper_body", 28) = {4};
//+
Physical Curve("right_edge_upper_body", 29) = {3};
//+
Physical Curve("left_edge_upper_body", 30) = {5};