lc = 25.0;
Point(11) = {-100, 0, 0, lc};
Point(12) = {-100, 100, 0, lc};
Point(1) = {0, 0, 0, lc};
Point(2) = {3000, 0, 0, lc};
Point(3) = {3000, 100, 0, lc};
Point(4) = {0, 100, 0, lc};
Point(5) = {3100, 0, 0, lc};
Point(6) = {3100, 100, 0, lc};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
Line Loop(1) = {4, 1, 2, 3};
Plane Surface(1) = {1};
Transfinite Surface {1};
Recombine Surface {1};//+
Line(5) = {2, 5};
Line(6) = {5, 6};
Line(7) = {6, 3};
Line Loop(2) = {2, -7, -6, -5};
Plane Surface(2) = {2};
Transfinite Surface {2};
Recombine Surface {2};

//+
Line(8) = {4, 12};
//+
Line(9) = {12, 11};
//+
Line(10) = {11, 11};
//+
Line(11) = {11, 1};
//+
Line Loop(3) = {8, 9, 11, -4};
//+
Plane Surface(3) = {3};
Transfinite Surface {3};
Recombine Surface {3};
//+
Physical Point("all_points") = {11, 12, 4, 1, 3, 2, 5, 6};
//+
Physical Line("all_lines") = {9, 8, 4, 11, 3, 1, 7, 2, 5, 6};
//+
Physical Surface("pmlL") = {3};
//+
Physical Surface("waveguide") = {1};
//+
Physical Surface("pmlR") = {2};
