// geometry for 3D example b-conforming

d=3.e-4;

Point(1) = {0.0, 0, 0, d*4} ;
Point(2) = {0.01, 0, 0, d*4} ;
Point(5) = {0.01, 0.02, 0, d*4} ;
Point(6) = {0.0, 0.02, 0, d*4} ;
Point(3) = {0.0, 0.004, 0, d*2} ;
Point(4) = {0.01, 0.004, 0, d} ;//
Point(9) = {0.0, 2*0.004, 0, d*4} ;
Point(10) = {0.01-0.004, 2*0.004, 0, d*2} ;
Point(11) = {0.01-0.004, 0.004+0.001, 0, d} ;//8
Point(12) = {0.01, 0.004+0.001, 0, d} ;//*
Point(13) = {0.0, 0.02-0.004, 0, d*2} ;
Point(14) = {0.01-0.004, 0.02-0.004, 0, d*4} ;

Line(201) = {1,2};
Line(202) = {2,4};
Line(203) = {4,3};
Line(204) = {3,1};
Line(207) = {3,9};
Line(208) = {9,10};
Line(209) = {10,11};
Line(210) = {11,12};
Line(211) = {12,4};
Line(212) = {12,5};
Line(213) = {5,6};
Line(214) = {6,13};
Line(215) = {10,14};
Line(216) = {14,13};
Line(217) = {13,9};

Line Loop(218) = {213,214,-216,-215,209,210,212};
Plane Surface(219) = {218};
Line Loop(220) = {217,208,215,216};
Plane Surface(221) = {220};
Line Loop(222) = {207,208,209,210,211,203};
Plane Surface(223) = {222};
Line Loop(224) = {204,201,202,203};
Plane Surface(225) = {224};

Extrude Surface {221, {0.0,0.0,0.02}};
Extrude Surface {219, {0.0,0.0,0.02}};
Extrude Surface {225, {0.0,0.0,0.02}};
Extrude Surface {223, {0.0,0.0,0.02}};

Physical Surface(1010) = {259,283,284,247,338,306,297,301,333,219,221,225,223}; //Dirichlet

Physical Volume(801) = {2}; //Inductor
Physical Volume(802) = {1}; //Core
Physical Volume(803) = {4}; //Armature
Physical Volume(804) = {3}; //Air