raio = 300;
lesfera = 3;


lx = 1368;
ly = 1368;
lz = 4446;
lc = 150;


//esfera no centro
Point(1) = {0, 0, 0, lesfera}; //Centro
Point(8) = {raio, 0, 0, lesfera}; Point(9) = {-raio, 0, 0, lesfera};
Point(10) = {0, raio, 0, lesfera}; Point(11) = {0, -raio, 0, lesfera};
Point(12) = {0, 0, raio, lesfera}; Point(13) = {0, 0, -raio, lesfera};


// esfera interna
Circle(30) = {12, 1, 11}; Circle(31) = {11, 1, 13};
Circle(32) = {12, 1, 8}; Circle(33) = {8, 1, 13};
Circle(34) = {12, 1, 10}; Circle(35) = {10, 1, 13};
Circle(36) = {12, 1, 9}; Circle(37) = {9, 1, 13};
Circle(38) = {11, 1, 8}; Circle(39) = {8, 1, 10};
Circle(40) = {10, 1, 9}; Circle(41) = {9, 1, 11};

Line Loop(42) = {30, 38, -32}; Ruled Surface(43) = {42};
Line Loop(44) = {32, 39, -34}; Ruled Surface(45) = {44};
Line Loop(46) = {33, -31, 38}; Ruled Surface(47) = {46};
Line Loop(48) = {39, 35, -33}; Ruled Surface(49) = {48};
Line Loop(50) = {35, -37, -40}; Ruled Surface(51) = {50};
Line Loop(52) = {40, -36, 34}; Ruled Surface(53) = {52};
Line Loop(54) = {41, -30, 36}; Ruled Surface(55) = {54};
Line Loop(56) = {41, 31, -37}; Ruled Surface(57) = {56};

Surface Loop(57) = {57, 55, 43, 47, 49, 45, 53, 51};
Physical Surface(58) = {57, 55, 43, 47, 49, 45, 53, 51};

Volume(58) = {57};
Physical Volume(58) = {58};



Point(10001) = { - lx, - ly, - lz, lc};
Point(10002) = {   lx, - ly, - lz, lc};
Point(10003) = {   lx,   ly, - lz, lc};
Point(10004) = { - lx,   ly, - lz, lc};


Line(10001) = {10004, 10003};
Line(10002) = {10003, 10002};
Line(10003) = {10002, 10001};
Line(10004) = {10001, 10004};
Line Loop(10005) = {10002, 10003, 10004, 10001};
Plane Surface(10006) = {10005};
tmp[] = Extrude {0, 0, 2*lz} {
  Surface{10006};
};
Delete { Volume{tmp[1]}; }

Surface Loop(10029) = {10027,10006,10015,10019,10023,10028};
Volume(10031) = {10029,57};

Physical Surface(10033) = {10027}; // Cima
Physical Surface(10034) = {10019}; // Baixo

Physical Surface(10035) = {10023}; // Esquerda
Physical Surface(10036) = {10015}; // Direita

Physical Surface(10037) = {10006}; // Atrás
Physical Surface(10038) = {10028}; // Frente


Physical Volume(10039) = {10031};

Characteristic Length {10001:10014} = lc;










