// definition des différents points utilisés   
lcar3 = 0.00411;  
// the bottom edge of the square coord2   
p1 = newp; Point(p1) = { 0  , 0.116623  , 0  , 0.00411  }; 
p2 = newp; Point(p2) = { 0.0697396  , 0.0449548  , 0  , 0.00411  }; 
p3 = newp; Point(p3) = { 0.159065  , 0  , 0  , 0.00411  }; 
p4 = newp; Point(p4) = { 0.239659  , 0  , 0  , 0.00411  }; 
p5 = newp; Point(p5) = { 0.296476  , 0.082291  , 0  , 0.00411  }; 
p6 = newp; Point(p6) = { 0.296476  , 0.182291  , 0  , 0.00411  }; 
p7 = newp; Point(p7) = { 0.353293  , 0  , 0  , 0.00411  }; 
p8 = newp; Point(p8) = { 0.550038  , 0  , 0  , 0.00411  }; 
p9 = newp; Point(p9) = { 0.633544  , -0.0550167  , 0  , 0.00411  }; 
p10 = newp; Point(p10) = { 0.633544  , 0.0449833  , 0  , 0.00411  }; 
p11 = newp; Point(p11) = { 0.717049  , 0  , 0  , 0.00411  }; 
// the right side  edge of the square coord3  
p12 = newp; Point(p12) = { 0.980414  , 0  , 0  , 0.00411  }; 
p13 = newp; Point(p13) = { 1.06974  , 0.0449548  , 0  , 0.00411  }; 
p14 = newp; Point(p14) = { 1  , 0.116623  , 0  , 0.00411  }; 
p15 = newp; Point(p15) = { 1  , 0.437302  , 0  , 0.00411  }; 
p16 = newp; Point(p16) = { 1.0818  , 0.494826  , 0  , 0.00411  }; 
p17 = newp; Point(p17) = { 0.981799  , 0.494826  , 0  , 0.00411  }; 
p18 = newp; Point(p18) = { 1  , 0.552349  , 0  , 0.00411  }; 
p19 = newp; Point(p19) = { 1  , 0.716771  , 0  , 0.00411  }; 
p20 = newp; Point(p20) = { 1.00587  , 0.816599  , 0  , 0.00411  }; 
p21 = newp; Point(p21) = { 0.905871  , 0.816599  , 0  , 0.00411  }; 
p22 = newp; Point(p22) = { 1  , 0.916426  , 0  , 0.00411  }; 
//  the top edge of the square coord1 
p23 = newp; Point(p23) = { 1  , 0.973286  , 0  , 0.00411  }; 
p24 = newp; Point(p24) = { 1.06974  , 1.04495  , 0  , 0.00411  }; 
p25 = newp; Point(p25) = { 0.980414  , 1  , 0  , 0.00411  }; 
p26 = newp; Point(p26) = { 0.717049  , 1  , 0  , 0.00411  }; 
p27 = newp; Point(p27) = { 0.633544  , 0.844983  , 0  , 0.00411  }; 
p28 = newp; Point(p28) = { 0.633544  , 0.944983  , 0  , 0.00411  }; 
p29 = newp; Point(p29) = { 0.550038  , 1  , 0  , 0.00411  }; 
p30 = newp; Point(p30) = { 0.353293  , 1  , 0  , 0.00411  }; 
p31 = newp; Point(p31) = { 0.296476  , 0.982291  , 0  , 0.00411  }; 
p32 = newp; Point(p32) = { 0.296476  , 1.08229  , 0  , 0.00411  }; 
p33 = newp; Point(p33) = { 0.239659  , 1  , 0  , 0.00411  }; 
// the top edge of the square coord4 
p34 = newp; Point(p34) = { 0.159065  , 1  , 0  , 0.00411  }; 
p35 = newp; Point(p35) = { 0.0697396  , 1.04495  , 0  , 0.00411  }; 
p36 = newp; Point(p36) = { 0  , 0.973286  , 0  , 0.00411  }; 
p37 = newp; Point(p37) = { 0  , 0.916426  , 0  , 0.00411  }; 
p38 = newp; Point(p38) = { 0.105871  , 0.816599  , 0  , 0.00411  }; 
p39 = newp; Point(p39) = { 0.00587081  , 0.816599  , 0  , 0.00411  }; 
p40 = newp; Point(p40) = { 0  , 0.716771  , 0  , 0.00411  }; 
p41 = newp; Point(p41) = { 0  , 0.552349  , 0  , 0.00411  }; 
p42 = newp; Point(p42) = { 0.181799  , 0.494826  , 0  , 0.00411  }; 
p43 = newp; Point(p43) = { 0.0817989  , 0.494826  , 0  , 0.00411  }; 
p44 = newp; Point(p44) = { 0  , 0.437302  , 0  , 0.00411  }; 
p45 = newp; Point(p45) = { 0 ,1.166233e-001 ,0 ,lcar3 };
// the body of the square coord9 
p46 = newp; Point(p46) = { 0.398677  , 0.301152  , 0  , 0.0041  }; 
p47 = newp; Point(p47) = { 0.498677  , 0.301152  , 0  , 0.0041  }; 
p48 = newp; Point(p48) = { 0.398677  , 0.401152  , 0  , 0.0041  }; 
p49 = newp; Point(p49) = { 0.298677  , 0.301152  , 0  , 0.0041  }; 
p50 = newp; Point(p50) = { 0.398677  , 0.201152  , 0  , 0.0041  }; 
p51 = newp; Point(p51) = { 0.260841  , 0.704013  , 0  , 0.0041  }; 
p52 = newp; Point(p52) = { 0.360841  , 0.704013  , 0  , 0.0041  }; 
p53 = newp; Point(p53) = { 0.260841  , 0.804013  , 0  , 0.0041  }; 
p54 = newp; Point(p54) = { 0.160841  , 0.704013  , 0  , 0.0041  }; 
p55 = newp; Point(p55) = { 0.260841  , 0.604013  , 0  , 0.0041  }; 
p56 = newp; Point(p56) = { 0.636348  , 0.532735  , 0  , 0.0041  }; 
p57 = newp; Point(p57) = { 0.736348  , 0.532735  , 0  , 0.0041  }; 
p58 = newp; Point(p58) = { 0.636348  , 0.632735  , 0  , 0.0041  }; 
p59 = newp; Point(p59) = { 0.536348  , 0.532735  , 0  , 0.0041  }; 
p60 = newp; Point(p60) = { 0.636348  , 0.432735  , 0  , 0.0041  }; 
p61 = newp; Point(p61) = { 0.822221  , 0.411498  , 0  , 0.0041  }; 
p62 = newp; Point(p62) = { 0.922221  , 0.411498  , 0  , 0.0041  }; 
p63 = newp; Point(p63) = { 0.822221  , 0.511498  , 0  , 0.0041  }; 
p64 = newp; Point(p64) = { 0.722221  , 0.411498  , 0  , 0.0041  }; 
p65 = newp; Point(p65) = { 0.822221  , 0.311498  , 0  , 0.0041  }; 
p66 = newp; Point(p66) = { 0.163245  , 0.276342  , 0  , 0.0041  }; 
p67 = newp; Point(p67) = { 0.263245  , 0.276342  , 0  , 0.0041  }; 
p68 = newp; Point(p68) = { 0.163245  , 0.376342  , 0  , 0.0041  }; 
p69 = newp; Point(p69) = { 0.0632451  , 0.276342  , 0  , 0.0041  }; 
p70 = newp; Point(p70) = { 0.163245  , 0.176342  , 0  , 0.0041  }; 
p71 = newp; Point(p71) = { 0.377614  , 0.51972  , 0  , 0.0041  }; 
p72 = newp; Point(p72) = { 0.477614  , 0.51972  , 0  , 0.0041  }; 
p73 = newp; Point(p73) = { 0.377614  , 0.61972  , 0  , 0.0041  }; 
p74 = newp; Point(p74) = { 0.277614  , 0.51972  , 0  , 0.0041  }; 
p75 = newp; Point(p75) = { 0.377614  , 0.41972  , 0  , 0.0041  }; 
p76 = newp; Point(p76) = { 0.602238  , 0.182839  , 0  , 0.0041  }; 
p77 = newp; Point(p77) = { 0.702238  , 0.182839  , 0  , 0.0041  }; 
p78 = newp; Point(p78) = { 0.602238  , 0.282839  , 0  , 0.0041  }; 
p79 = newp; Point(p79) = { 0.502238  , 0.182839  , 0  , 0.0041  }; 
p80 = newp; Point(p80) = { 0.602238  , 0.0828386  , 0  , 0.0041  }; 
// tracé des cercles de la zone 2 
C1=newreg;Circle(C1)={p1,p2,p3};
l2=newreg;Line(l2)={p3,p4};
C3=newreg;Circle(C3)={p4,p5,p6};
C4=newreg;Circle(C4)={p6,p5,p7};
l5=newreg;Line(l5)={p7,p8};
C6=newreg;Circle(C6)={p8,p9,p10};
C7=newreg;Circle(C7)={p10,p9,p11};
l8=newreg;Line(l8)={p11,p12};
// tracé des cercles de la zone 3 
C9=newreg;Circle(C9)={p12,p13,p14};
l10=newreg;Line(l10)={p14,p15};
C11=newreg;Circle(C11)={p15,p16,p17};
C12=newreg;Circle(C12)={p17,p16,p18};
l13=newreg;Line(l13)={p18,p19};
C14=newreg;Circle(C14)={p19,p20,p21};
C15=newreg;Circle(C15)={p21,p20,p22};
l16=newreg;Line(l16)={p22,p23};
// tracé des cercles de la zone 1 
C17=newreg;Circle(C17)={p23,p24,p25};
l18=newreg;Line(l18)={p25,p26};
C20=newreg;Circle(C20)={p26,p28,p27};
C21=newreg;Circle(C21)={p27,p28,p29};
l22=newreg;Line(l22)={p29,p30};
C23=newreg;Circle(C23)={p30,p32,p31};
C24=newreg;Circle(C24)={p31,p32,p33};
l25=newreg;Line(l25)={p33,p34};
// tracé des cercles de la zone 4 
C26=newreg;Circle(C26)={p34,p35,p36};
l27=newreg;Line(l27)={p36,p37};
C28=newreg;Circle(C28)={p37,p39,p38};
C29=newreg;Circle(C29)={p38,p39,p40};
l30=newreg;Line(l30)={p40,p41};
C31=newreg;Circle(C31)={p41,p43,p42};
C32=newreg;Circle(C32)={p42,p43,p44};
l33=newreg;Line(l33)={p44,p45};
l34=newreg;Line(l34)={p44,p1};
// tracé des cercles de la zone 9 
C34=newreg;Circle(C34)={p47,p46,p48};
C35=newreg;Circle(C35)={p48,p46,p49};
C36=newreg;Circle(C36)={p49,p46,p50};
C37=newreg;Circle(C37)={p50,p46,p47};
Line Loop(1)={C34,C35,C36,C37};
//Plane Surface(2)={1};
C38=newreg;Circle(C38)={p52,p51,p53};
C39=newreg;Circle(C39)={p53,p51,p54};
C40=newreg;Circle(C40)={p54,p51,p55};
C41=newreg;Circle(C41)={p55,p51,p52};
Line Loop(3)={C38,C39,C40,C41};
//Plane Surface(4)={3};
C42=newreg;Circle(C42)={p57,p56,p58};
C43=newreg;Circle(C43)={p58,p56,p59};
C44=newreg;Circle(C44)={p59,p56,p60};
C45=newreg;Circle(C45)={p60,p56,p57};
Line Loop(5)={C42,C43,C44,C45};
//Plane Surface(6)={5};
C46=newreg;Circle(C46)={p62,p61,p63};
C47=newreg;Circle(C47)={p63,p61,p64};
C48=newreg;Circle(C48)={p64,p61,p65};
C49=newreg;Circle(C49)={p65,p61,p62};
Line Loop(7)={C46,C47,C48,C49};
//Plane Surface(8)={7};
C50=newreg;Circle(C50)={p67,p66,p68};
C51=newreg;Circle(C51)={p68,p66,p69};
C52=newreg;Circle(C52)={p69,p66,p70};
C53=newreg;Circle(C53)={p70,p66,p67};
Line Loop(9)={C50,C51,C52,C53};
//Plane Surface(10)={9};
C54=newreg;Circle(C54)={p72,p71,p73};
C55=newreg;Circle(C55)={p73,p71,p74};
C56=newreg;Circle(C56)={p74,p71,p75};
C57=newreg;Circle(C57)={p75,p71,p72};
Line Loop(11)={C54,C55,C56,C57};
//Plane Surface(12)={11};
C58=newreg;Circle(C58)={p77,p76,p78};
C59=newreg;Circle(C59)={p78,p76,p79};
C60=newreg;Circle(C60)={p79,p76,p80};
C61=newreg;Circle(C61)={p80,p76,p77};
Line Loop(13)={C58,C59,C60,C61};
//Plane Surface(14)={13};
Line Loop(400)={C1,l2,C3,C4,l5,C6,C7,l8,C9,l10,C11,C12,l13,C14,C15,l16,C17,l18,C20,C21,l22,C23,C24,l25,C26,l27,C28,C29,l30,C31,C32,l34};
Plane Surface(401)={400,1,3,5,7,9,11,13};

h = 0.1;

Extrude Surface { 401, {0,0,h} } { 
//  Layers { {8,4,2,1}, {9000:9003}, {0.25,0.5,0.75,1} }; 
  Recombine; Layers { 1, 9000, 1 }; 
};