len = 1e-2;
dia = 2e-3;
lc = 1e-4;
odia = 2.5e-3;
Point(2) = {dia/2,0,0,lc};
Point(3) = {dia/2,len,0,lc};
Point(4) = {0,len,0,lc};
Point(10) = {odia/2,0,0,lc};
Point(11) = {odia/2,len,0,lc};
Point(12) = {10*odia,0,0,lc*8};
Point(14) = {0,10*odia,0,lc*8};
Point(15) = {0,0,0,lc/2};
Line(5) = {3,11};
Line(6) = {11,10};
Line(7) = {10,2};
Line(8) = {2,3};



//Line(11) = {15,2};

Circle(12) = {12,15,14};

Point(16) = {0,6*odia,0,lc*4};
Point(17) = {6*odia,0,0,lc*4};

Circle(14) = {17,15,16};

Line(18) = {10,17};
Line(19) = {17,12};
Line(20) = {16,14};
//Line(40) = {15,4};
Line(41) = {4,16};
Line(42) = {4,3};


Point(103) = {dia/2,-len,0,lc};
Point(104) = {0,-len,0,lc};
Point(105) = {odia/2,-len,0,lc};
Point(106) = {0,-6*odia,0,lc*4};
Point(107) = {0,-10*odia,0,lc*4};

Line Loop(43) = {-14,19,12,-20};
Plane Surface(44) = {43};
//Line Loop(45) = {18,14,-41,42,5,6};
//Plane Surface(46) = {45};
//Line Loop(47) = {-8,-11,40,42};
//Plane Surface(48) = {47};
Line Loop(49) = {5,6,7,8};
Plane Surface(50) = {49};

//Line(56) = {15,104};
Line(57) = {104,103};
Line(58) = {103,105};
Line(59) = {105,10};
Line(60) = {104,106};
Line(61) = {106,107};
Circle(62) = {106,15,17};
Circle(63) = {107,15,12};
Line(64) = {103,2};

//Line Loop(65) = {-11,56,57,64};
//Plane Surface(66) = {65};
Line Loop(67) = {7,-64,58,59};
Plane Surface(68) = {67};
Line Loop(69) = {62,-18,-59,-58,-57,60};
Plane Surface(70) = {69};
Line Loop(71) = {-19,-62,61,63};
Plane Surface(72) = {71};

c_dia = 0.0002;
Point(108) = {0.005,0,0,lc};
Point(109) = {c_dia/2,0,0,lc/4};
Point(110) = {0,c_dia/2,0,lc/4};
Point(111) = {0,-c_dia/2,0,lc/4};
Circle(73) = {111,15,109};
Circle(74) = {109,15,110};
Line(75) = {104,111};
Line(76) = {15,110};
Line(77) = {111,15};
Line(78) = {110,4};
Line(79) = {109,2};
Line Loop(81) = {64,-79,-73,-75,57};
Plane Surface(82) = {81};
Line Loop(83) = {8,-42,-78,-74,79};
Plane Surface(84) = {83};
Line Loop(85) = {-76,-77,73,74};
Plane Surface(86) = {85};
Line Loop(87) = {-14,-18,-6,-5,-42,41};
Plane Surface(88) = {87};

Physical Surface(51) = {44,72};   //outer sphere air infinity
Physical Surface(52) = {88,70};   //inner sphere air
Physical Surface(53) = {84,82};   //core
Physical Surface(54) = {50,68};   //coil
Physical Surface(57) = {86};   //central sphere
Physical Line(55) = {12,63};   // outer boundary
Physical Line(56) = {20,41,75,76,77,78,60,61};   //central axis
//Point(1305) = {0,len/5,0,lc};
//Point(1306) = {0,-len/5,0,lc};
