Include "wave_pml_param.pro";
Mesh.CharacteristicLengthFactor = 1;

// Volume de propagation
Point(100) = {a,0,0,lc};
Point(101) = {0,a,0,lc};
Point(1) = {0,0,0,lc};
Point(2) = {Lx,0,0,lc};
Point(3) = {Lx,Ly,0,lc};
Point(4) = {0,Ly,0,lc};
Line(1) = {100,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,101};
Circle(102) = {101,1,100};
Line Loop(5) = {1,2,3,4,102};
Plane Surface(1) = {5};

// Volume PML
Point(5) = {Lx+PmlThickness,0,0,lc};
Point(6) = {Lx+PmlThickness,Ly+PmlThickness,0,lc};
Point(7) = {0,Ly+PmlThickness,0,lc};
Line(6) = {2,5};
Line(7) = {5,6};
Line(8) = {6,7};
Line(9) = {7,4};
Line Loop(12) = {6,7,8,9,-3,-2};
Plane Surface(2) = {12};


Physical Surface(VolumePropaIndex) = {1};
Physical Surface(VolumePmlIndex) = {2};
Physical Line(SurfaceGaucheIndex) = {9,4};
Physical Line(SurfaceHautIndex) = {8};
Physical Line(SurfaceBasIndex) = {1,6};
Physical Line(SurfaceDroiteIndex) = {7};
Physical Line(SurfaceSourceIndex) = {102};
