// Parameters
l1=1.0;
l2=0.01;

// Points lower body
Point(1) = {-2, 0, 0, l1};
Point(2) = {2, 0, 0, l1};
Point(3) = {2, 1, 0, l1};
Point(4) = {-2, 1, 0, l1};

// Points upper body
Point(5) = {-1, 1.05, 0, l2};
Point(6) = {1, 1.05, 0, l2};
Point(7) = {1, 2.0, 0, l2};
Point(8) = {-1, 2.0, 0, l2};

// Lines
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};

// Loops
Curve Loop(1) = {3, 4, 1, 2};
Curve Loop(2) = {8, 5, 6, 7};

// Surfaces
Plane Surface(1) = {1};
Plane Surface(2) = {2};

// Physical groups
Physical Surface("lower_body", 9) = {1};
Physical Surface("upper_body", 10) = {2};
Physical Curve("lower_body_bottom_edge", 11) = {1};
Physical Curve("lower_body_right_edge", 12) = {2};
Physical Curve("lower_body_top_edge", 13) = {3};
Physical Curve("lower_body_left_edge", 14) = {4};
Physical Curve("top_body_bottom_edge", 15) = {5};
Physical Curve("top_body_right_edge", 16) = {6};
Physical Curve("top_body_top_edge", 17) = {7};
Physical Curve("top_body_left_edge", 18) = {8};