//+
SetFactory("OpenCASCADE");
Circle(1) = {0, -0, 0, 1, 0, 2*Pi};
//+
Circle(2) = {-0, -0, 0, 0.8, 0, 2*Pi};
//+
Circle(3) = {0, 0.9, 0, 0.05, 0, 2*Pi};
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/16} {
  Duplicata { Curve{3}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{4}; Curve{3}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{5}; Curve{6}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{7}; Curve{8}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{10}; Curve{9}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{11}; Curve{12}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{13}; Curve{14}; }
}
//+
Rotate {{0, 0, 1}, {0, 0, 0}, Pi/4} {
  Duplicata { Curve{15}; Curve{16}; }
}
//+
Circle(19) = {-0, -0, 0, 5, 0, 2*Pi};
//+
Circle(20) = {0, -0, 0, 4, 0, 2*Pi};
//+
//+
Physical Curve("dielectric_inner") = {2};
//+
Physical Curve("dielectric_outer") = {1};
//+
Physical Curve("pml_outer") = {20};
//+
Physical Curve("pml_inner") = {19};
//+
Physical Curve("airholerings") = {4, 3, 18, 17, 16, 15, 14, 13, 12, 11, 9, 10, 7, 8, 5, 6};