Include "parameters_gmsh_getdp.dat";
cl__1 = lambda0/15;
cl__2 = lambda0/60;
cl__3 = lambda0/80;
Point(1) = {-2*lambda0, -2*lambda0, 0, cl__1};
Point(2) = {2*lambda0, -2*lambda0, 0, cl__1};
Point(3) = {2*lambda0, 2*lambda0, 0, cl__1};
Point(4) = {-2*lambda0, 2*lambda0, 0, cl__1};
Point(5) = {-3*lambda0, -2*lambda0, 0, cl__1};
Point(6) = {-2*lambda0, -3*lambda0, 0, cl__1};
Point(7) = {2*lambda0, -3*lambda0, 0, cl__1};
Point(8) = {-3*lambda0, -3*lambda0, 0, cl__1};
Point(9) = {3*lambda0, -3*lambda0, 0, cl__1};
Point(10) = {3*lambda0, -2*lambda0, 0, cl__1};
Point(11) = {3*lambda0, 2*lambda0, 0, cl__1};
Point(12) = {3*lambda0, 3*lambda0, 0, cl__1};
Point(13) = {2*lambda0, 3*lambda0, 0, cl__1};
Point(14) = {-2*lambda0, 3*lambda0, 0, cl__1};
Point(15) = {-3*lambda0, 3*lambda0, 0, cl__1};
Point(16) = {-3*lambda0, 2*lambda0, 0, cl__1};
Point(17) = {x_source, y_source, 0, cl__3};
Point(18) = {x_source           , y_source-lambda0/6 , 0, cl__2};
Point(19) = {x_source+lambda0/6 , y_source           , 0, cl__2};
Point(20) = {x_source           , y_source+lambda0/6 , 0, cl__2};
Point(21) = {x_source-lambda0/6 , y_source           , 0, cl__2};

Line(1) = {14, 15};
Line(2) = {15, 16};
Line(3) = {16, 4};
Line(4) = {4, 14};
Line(5) = {16, 5};
Line(6) = {5, 8};
Line(7) = {8, 6};
Line(8) = {5, 1};
Line(9) = {1, 6};
Line(10) = {1, 4};
Line(11) = {4, 3};
Line(12) = {3, 13};
Line(13) = {14, 13};
Line(14) = {13, 12};
Line(15) = {12, 11};
Line(16) = {11, 3};
Line(17) = {3, 2};
Line(18) = {2, 10};
Line(19) = {10, 11};
Line(20) = {2, 7};
Line(21) = {9, 7};
Line(22) = {9, 10};
Line(23) = {1, 2};
Line(24) = {6, 7};

Circle(41) = {18, 17, 19};
Circle(42) = {19, 17, 20};
Circle(43) = {20, 17, 21};
Circle(44) = {21, 17, 18};

Line Loop(26) = {2, 3, 4, 1};
Plane Surface(26) = {26};
Line Loop(28) = {4, 13, -12, -11};
Plane Surface(28) = {28};
Line Loop(30) = {12, 14, 15, 16};
Plane Surface(30) = {30};
Line Loop(32) = {17, 18, 19, 16};
Plane Surface(32) = {32};
Line Loop(34) = {22, -18, 20, -21};
Plane Surface(34) = {34};
Line Loop(36) = {23, 20, -24, -9};
Plane Surface(36) = {36};
Line Loop(38) = {7, -9, -8, 6};
Plane Surface(38) = {38};
Line Loop(40) = {10, -3, 5, 8};
Plane Surface(40) = {40};

Line Loop(45) = {43, 44, 41, 42};
Plane Surface(46) = {45};
Line Loop(47) = {10, 11, 17, -23};
Plane Surface(48) = {45, 47};



Point{17} In Surface {46};


Physical Surface(1) = {26,30,34,38}; // PML corners
Physical Surface(2) = {28,36};       // PML top bot
Physical Surface(3) = {32,40};	// PML left right
Physical Surface(4) = {46};          // air around source
Physical Surface(5) = {48};          // air2
Physical Point(9)   = {17}; 			  // source point
Physical Point(10)   = {1}; 				 // print point
