// #1, isotropic material
// Material constants are arranged as:
// temperature density 
// E v
Physical Point("1 0 1
0 0
1900 0.38
0 0
") = {};
// Material name: 1 -- PCABS

cl = 0.5;
Point(1) = {0, 1, -0, cl};
Point(2) = {0, -1, -0, cl};
Point(3) = {2, 0, -0, cl};
Point(4) = {-2, 0, -0, cl};
Point(5) = {-3, 1.5, -0, cl};
Point(6) = {3, 1.5, -0, cl};
Point(7) = {3, -1.5, -0, cl};
Point(8) = {-3, -1.5, -0, cl};
Line(2) = {8, 7};
Line(3) = {7, 6};
Line(4) = {6, 5};
Line(5) = {5, 8};
Ellipse(6) = {2, 7, 2, 3};
Ellipse(7) = {3, 6, 1, 1};
Ellipse(8) = {1, 5, 1, 4};
Ellipse(9) = {4, 8, 2, 2};
Line Loop(12) = {4, 5, 2, 3, -7, -6, -9, -8};
Plane Surface(12) = {12};
Physical Surface(1) = {12};

Periodic Line {5}={-3};
Periodic Line {2}={-4};

Mesh.Points = 1;
Mesh.SurfaceFaces=1;
Mesh.SurfaceEdges=1;
Mesh.VolumeEdges=1;
Mesh.ColorCarousel=2;


Mesh.Points = 1;
Mesh.SurfaceFaces=1;
Mesh.SurfaceEdges=1;
Mesh.VolumeEdges=1;
Mesh.ColorCarousel=2;