 
 //coil
 pcoil = newp;
-Point (pcoil+1) = {9, 11, -24.6} ;
-Point (pcoil+2) = {-9, 11, -24.6} ;
-Point (pcoil+3) = {-9, -11, -24.6} ;
-Point (pcoil+4) = {9, -11, -24.6} ;
+Point (pcoil+1) = {9, 11, -24.6,lc} ;
+Point (pcoil+2) = {-9, 11, -24.6,lc} ;
+Point (pcoil+3) = {-9, -11, -24.6,lc} ;
+Point (pcoil+4) = {9, -11, -24.6,lc} ;
 
 lcoil = newl;
 Line (lcoil+1) = {pcoil+1, pcoil+2} ;

//SetFactory("OpenCASCADE");

halfWidth = 8.05;
//down part
//left
p1 = newp;
Point(p1+1) = {0, -halfWidth, -34};
Point(p1+2) = {0, -halfWidth, -12.15};
Point(p1+3) = {0.98, -halfWidth, -12.15};
Point(p1+4) = {4.825, -halfWidth, 2.2};
Point(p1+5) = {5.95, -halfWidth, 2.2};
Point(p1+6) = {5.95, -halfWidth, -28};
Point(p1+7) = {17.1, -halfWidth, -28};
Point(p1+8) = {17.1, -halfWidth, 22.25};
Point(p1+9) = {17.686, -halfWidth, 22.25};
Point(p1+10) = {23, -halfWidth, 7.65};
Point(p1+11) = {25.85, -halfWidth, 7.65};
Point(p1+12) = {25.85, -halfWidth, 5.475};
Point(p1+13) = {23, -halfWidth, 5.475};
Point(p1+14) = {23, -halfWidth, 3.3};
Point(p1+15) = {25.85, -halfWidth, 3.3};
Point(p1+16) = {25.85, -halfWidth, -2.557};
Point(p1+17) = {23.25, -halfWidth, -9.7};
Point(p1+18) = {23.25, -halfWidth, -34};

l1=newl;
//Line (l1+1) = {p1+1, p1+2} ;


//duichen
p22 = newp;
//Point(p22+1) = {0, halfWidth, -11.15,lc};
Point(p22+2) = {-0.595, halfWidth, -11.15,lc};
Point(p22+3) = {-4.44, halfWidth, 3.2,lc};
Point(p22+4) = {-5.5, halfWidth, 3.2,lc};
Point(p22+5) = {-5.5, halfWidth, 30,lc};
Point(p22+6) = {-18, halfWidth, 30,lc};
Point(p22+7) = {-18, halfWidth, 23.625,lc};
Point(p22+8) = {-23.45, halfWidth, 8.65,lc};
Point(p22+9) = {-25.35, halfWidth, 8.65,lc};
Point(p22+10) = {-25.35, halfWidth, 25,lc};
Point(p22+11) = {-23.2, halfWidth, 25,lc};
Point(p22+12) = {-23.2, halfWidth, 35,lc};
//Point(p22+13) = {-0, halfWidth, 35,lc};

l33 = newl;
//Line (l33+1) = {p22+1, p22+2} ;
