lc = 0.34;
r1 = 1.0;
r2 = 2.0;
xc = 0.0;
yc = 0.0;
zc = 0.0;

//Sphere
Point(1) = {xc,yc,zc,lc};
Point(2) = {xc,yc,zc-r1,lc};
Point(3) = {xc,yc-r1,zc,lc};
Point(4) = {xc,yc,zc+r1,lc};

Point(5) = {xc,yc,zc-r2,lc};
Point(6) = {xc,yc-r2,zc,lc};
Point(7) = {xc,yc,zc+r2,lc};

Line(1) = {1,2};
Circle(2) = {2,1,3};
Circle(3) = {3,1,4};
Line(4) = {4,1};

Circle(5) = {5,1,6};
Circle(6) = {6,1,7};

Line(7) = {2, 5};
Line(8) = {4, 7};

Line Loop(1) = {2, 3, 8, -6, -5, -7};
Plane Surface(1) = {1};

Extrude {{0, 0, 1}, {0, 0, 0}, Pi/4} {  Surface{1}; }
Delete {Volume{1};}

Physical Surface(1) = {1};
Physical Surface(2) = {30};
Physical Surface(3) = {21, 18};
Physical Surface(4) = {25, 28};

Surface Loop(1) = {1, 18, 21, 30, 25, 28};
Volume(1) = {1};
Physical Volume(1) = {1};

Mesh.SecondOrderIncomplete = 0;
