[Gmsh] multi-scale mesh issue

Christophe Geuzaine geuzaine at gmail.com
Thu Jul 21 00:19:18 CEST 2005


Jean-Francois Roussel wrote:
> Hi,
> 
> first of all congratulations for this nery nice software.
> 
> I am using GMSH 3D mesh-generator to mesh the "empty" space with an 
> object embedded in it (I do not mesh the object as often in FE, but its 
> complement).
> And I have a specific need, which is to make multi-scale meshing.
> For example, the computation box can be a few 100m big, while I want to 
> mesh the object at mm scale (or below if possible!).
> 
> First, if I do that brutally, simply having very small and very large 
> characteristic lengths I get:
> - a nice 2D mesh
> - but a degenerated 3D mesh: very long tetrahedra starting from 
> different small 2D triangles and ending at a common node in volume.
> I solved that by helping GMSH thanks to intermediate boxes with 
> intermediate characteristic lengths.
> They are a little long to define, but it works fine.
> (If ever you have a better idea, I'd appreciate).
> 
> The second point I am now stuck with, is that gmsh does not seem be able 
> to cope with arbitrary "multi-scale ratio".

Hello Jean-Francois - That's (unfortunately) a known problem with the 
current implementation of our 3D algorithm. I'll have a look at your 
file and see if we can improve on the current situation.

Best,

Christophe



> When the ratio of the box size to the smallest element size goes beyond 
> ~ 100 000, gmsh either crashes or enters a seemingly ever lasting loop 
> (in all cases with error messages, rather variable, you have a few 
> examples below).
> I tried some changes without much result:
> - gmsh version: 1.56 -> 1.60
> - parameters: characteristic length factor, random perturbation (NetGen 
> does not work, probably I should check my normals orientation)
> A change that brought some relative improvement was to move the region 
> of fine meshing around x~y~z~0.
> That made me think that the limitation might come from the machine 
> precision when making differences in intermediate computations... 
> (errors reports are in single precision...).
> 
> I attached the file geom8du.geo, which causes the problem.
> I made the following tests on this file
> - if I remove the largest volume (Volume(3)): it works fine
> - if I remove the smallest 2 volumes (Volume(30113) and Volume(30139)): 
> it works fine, just with a few classical warnings ("Coplanar points in 
> circum sphere computation")
> - if I remove the smallest 3 volumes (previous + Volume(30165)): it 
> works just fine
> - if I remove intermediate volumes: it does not work
> which makes me think:
> - my CAD is correct (since all my volumes can be meshed independently 
> (but not all together))
> - it is an issue of relative scales (suppressing smallest or largest 
> volumes improves, suppressing intermediate ones does not)
> 
> I also attached a slightly more complex file geom8eu.geo (the one I 
> initially worked with) in case the threshold for triggering the issue is 
> different on your machine and geom8du.geo works on yours (mine is 
> Windows - Pentium).
> 
> Thanks a lot for your help!
> 
> Jean-Francois Roussel
> 
> -------Error messages -------
> 
> A couple of messages that led to infinite loop or crash:
> 
> "Huh! Panic in CrSi"
> "Unable to add point... (will try later)"
> 
> A recovered error:
> 
> Info    : 12/13: Edge 4622->1589 => 0 division(s)
> Error   : Missing edge without any intersection 
> (0,1.08048e-06,-0.00146627) (0,-1.05818e-16,-0.0005)
> Info    : 13/13: Edge 1588->4622 => 1 division(s)
> Info    : Boundary faces recovery
> Info    : Volume = 5.488e+06
> Info    : ===================================================
>        : Number of missing faces = 14
>        : ===================================================
> Info    : 1/14: Face 1562 1563 4440
> Info    : Intersections left (the face contains 0 initial faces, divided 
> in 3 points)
> Info    : The face is divided in 3 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (0 <--> 1=2*(3-1)-3)
> Info    : 2/14: Face 1591 1592 4431
> Info    : Intersections left (the face contains 5 initial faces, divided 
> in 8 points)
> Info    : The face is divided in 8 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (5 <--> 6=2*(8-1)-8)
> Info    : 3/14: Face 1601 1602 4435
> Info    : Recoverable face (4 <--> 4=2*(6-1)-6)
> Info    : 4/14: Face 1592 1593 4431
> Info    : Recoverable face (7 <--> 7=2*(9-1)-9)
> Info    : 5/14: Face 1562 1563 4440
> Info    : Intersections left (the face contains 0 initial faces, divided 
> in 3 points)
> Info    : The face is divided in 3 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (0 <--> 1=2*(3-1)-3)
> Info    : 6/14: Face 1001 1602 4435
> Info    : Intersections left (the face contains 1 initial faces, divided 
> in 4 points)
> Info    : The face is divided in 4 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (1 <--> 2=2*(4-1)-4)
> Info    : 7/14: Face 1003 1004 4414
> Info    : Recoverable face (2 <--> 2=2*(4-1)-4)
> Info    : 8/14: Face 1592 1593 4629
> Info    : Intersections left (the face contains 23 initial faces, 
> divided in 13 points)
> Info    : The face is divided in 13 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (23 <--> 11=2*(13-1)-13)
> Info    : 9/14: Face 1001 1602 4637
> Info    : Intersections left (the face contains 1 initial faces, divided 
> in 4 points)
> Info    : The face is divided in 4 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (1 <--> 2=2*(4-1)-4)
> Info    : 10/14: Face 1588 1589 4622
> Info    : Intersections left (the face contains 1 initial faces, divided 
> in 4 points)
> Info    : The face is divided in 4 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (1 <--> 2=2*(4-1)-4)
> Info    : 11/14: Face 1588 1589 4622
> Info    : Intersections left (the face contains 1 initial faces, divided 
> in 4 points)
> Info    : The face is divided in 4 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (1 <--> 2=2*(4-1)-4)
> Info    : 12/14: Face 1003 1004 4642
> Info    : Recoverable face (2 <--> 2=2*(4-1)-4)
> Info    : 13/14: Face 1601 1602 4637
> Info    : Intersections left (the face contains 5 initial faces, divided 
> in 6 points)
> Info    : The face is divided in 6 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (5 <--> 4=2*(6-1)-6)
> Info    : 14/14: Face 1591 1592 4629
> Info    : Intersections left (the face contains 20 initial faces, 
> divided in 12 points)
> Info    : The face is divided in 12 points (simplexes added=0, deleted=0)
> Warning : *Unrecoverable* face (20 <--> 10=2*(12-1)-12)
> Info    : Volume after edge/face recovery = 5.488e+06
> Info    : Final check: 1 missing edges, 2 missing faces
> Info    : ===================================================
>        : Surface coherence results (number of intersections)
>        : 42 EV, 16 EE, 0 FV, 0 FF, 0 FE, 3 EEE, 0 EEEE
>        : ===================================================
> Error   : Could not restore all edges/faces
> Error   : Surface recovery failed (send a bug report with the geo file 
> to <gmsh at geuz.org>!)
> Info    : Mesh 3D... (final)
> Info    : Mesh 3D complete (27.66 s)
> 
> 
> ------------------------------------------------------------------------
> 
> Point(101) = {0, 0, -6, 10};
> Point(102) = {50, 0, -6, 10};
> Point(103) = {-50, 0, -6, 10};
> Point(104) = {0, 50, -6, 10};
> Point(105) = {0, -50, -6, 10};
> Point(106) = {0, 0, 94, 10};
> Point(107) = {0, 0, -106, 10};
> Point(201) = {-1, -1, -7, 0.3};
> Point(202) = {-1, -1, -5, 0.3};
> Point(203) = {-1, 1, -7, 0.3};
> Point(204) = {-1, 1, -5, 0.3};
> Point(205) = {1, -1, -7, 0.3};
> Point(206) = {1, -1, -5, 0.3};
> Point(207) = {1, 1, -7, 0.3};
> Point(208) = {1, 1, -5, 0.3};
> Point(301) = {0, -1, -4, 0.4};
> Point(302) = {0, -1, 6, 0.4};
> Point(303) = {0, 1, -4, 0.4};
> Point(304) = {0, 1, 6, 0.4};
> Point(305) = {-0.2, -1, -4, 0.4};
> Point(306) = {-0.2, -1, 6, 0.4};
> Point(307) = {-0.2, 1, -4, 0.4};
> Point(308) = {-0.2, 1, 6, 0.4};
> Point(401) = {0, -1, -8, 0.6};
> Point(402) = {0, -1, -18, 0.6};
> Point(403) = {0, 1, -8, 0.6};
> Point(404) = {0, 1, -18, 0.6};
> Point(405) = {-0.2, -1, -8, 0.6};
> Point(406) = {-0.2, -1, -18, 0.6};
> Point(407) = {-0.2, 1, -8, 0.6};
> Point(408) = {-0.2, 1, -18, 0.6};
> Point(1001) = {0, -0.01, 0.9995000000000001, 0.0007};
> Point(1002) = {0, -0.01, 1.0005, 0.0007};
> Point(1003) = {0, 0.01, 0.9995000000000001, 0.0007};
> Point(1004) = {0, 0.01, 1.0005, 0.0007};
> Point(1005) = {0, -0.02, 0.997, 0.00175};
> Point(1006) = {0, -0.02, 1.003, 0.00175};
> Point(1007) = {0, 0.02, 0.997, 0.00175};
> Point(1008) = {0, 0.02, 1.003, 0.00175};
> Point(1009) = {0.005590169943749475, -0.02, 0.997, 0.002625};
> Point(1010) = {0.005590169943749475, -0.02, 1.003, 0.002625};
> Point(1011) = {0.005590169943749475, 0.02, 0.997, 0.002625};
> Point(1012) = {0.005590169943749475, 0.02, 1.003, 0.002625};
> Point(1013) = {0, -0.04, 0.991, 0.004375};
> Point(1014) = {0, -0.04, 1.009, 0.004375};
> Point(1015) = {0, 0.04, 0.991, 0.004375};
> Point(1016) = {0, 0.04, 1.009, 0.004375};
> Point(1017) = {0.01397542485937369, -0.04, 0.991, 0.006562500000000001};
> Point(1018) = {0.01397542485937369, -0.04, 1.009, 0.006562500000000001};
> Point(1019) = {0.01397542485937369, 0.04, 0.991, 0.006562500000000001};
> Point(1020) = {0.01397542485937369, 0.04, 1.009, 0.006562500000000001};
> Point(1021) = {0, -0.08, 0.973, 0.0109375};
> Point(1022) = {0, -0.08, 1.027, 0.0109375};
> Point(1023) = {0, 0.08, 0.973, 0.0109375};
> Point(1024) = {0, 0.08, 1.027, 0.0109375};
> Point(1025) = {0.03493856214843422, -0.08, 0.973, 0.01640625};
> Point(1026) = {0.03493856214843422, -0.08, 1.027, 0.01640625};
> Point(1027) = {0.03493856214843422, 0.08, 0.973, 0.01640625};
> Point(1028) = {0.03493856214843422, 0.08, 1.027, 0.01640625};
> Point(1029) = {0, -0.16, 0.919, 0.02734375};
> Point(1030) = {0, -0.16, 1.081, 0.02734375};
> Point(1031) = {0, 0.16, 0.919, 0.02734375};
> Point(1032) = {0, 0.16, 1.081, 0.02734375};
> Point(1033) = {0.08734640537108554, -0.16, 0.919, 0.04101562500000001};
> Point(1034) = {0.08734640537108554, -0.16, 1.081, 0.04101562500000001};
> Point(1035) = {0.08734640537108554, 0.16, 0.919, 0.04101562500000001};
> Point(1036) = {0.08734640537108554, 0.16, 1.081, 0.04101562500000001};
> Point(1037) = {0, -0.32, 0.7569999999999999, 0.06835937500000001};
> Point(1038) = {0, -0.32, 1.243, 0.06835937500000001};
> Point(1039) = {0, 0.32, 0.7569999999999999, 0.06835937500000001};
> Point(1040) = {0, 0.32, 1.243, 0.06835937500000001};
> Point(1041) = {0.2183660134277138, -0.32, 0.7569999999999999, 0.1025390625};
> Point(1042) = {0.2183660134277138, -0.32, 1.243, 0.1025390625};
> Point(1043) = {0.2183660134277138, 0.32, 0.7569999999999999, 0.1025390625};
> Point(1044) = {0.2183660134277138, 0.32, 1.243, 0.1025390625};
> Point(1045) = {0, -0.64, 0.2709999999999999, 0.1708984375};
> Point(1046) = {0, -0.64, 1.729, 0.1708984375};
> Point(1047) = {0, 0.64, 0.2709999999999999, 0.1708984375};
> Point(1048) = {0, 0.64, 1.729, 0.1708984375};
> Point(1049) = {0.5459150335692846, -0.64, 0.2709999999999999, 0.2563476562500001};
> Point(1050) = {0.5459150335692846, -0.64, 1.729, 0.2563476562500001};
> Point(1051) = {0.5459150335692846, 0.64, 0.2709999999999999, 0.2563476562500001};
> Point(1052) = {0.5459150335692846, 0.64, 1.729, 0.2563476562500001};
> Point(1053) = {-2.5, -2.5, -24, 0.75};
> Point(1054) = {-2.5, -2.5, 12, 0.75};
> Point(1055) = {-2.5, 2.5, -24, 0.75};
> Point(1056) = {-2.5, 2.5, 12, 0.75};
> Point(1057) = {2.5, -2.5, -24, 0.75};
> Point(1058) = {2.5, -2.5, 12, 0.75};
> Point(1059) = {2.5, 2.5, -24, 0.75};
> Point(1060) = {2.5, 2.5, 12, 0.75};
> Point(1061) = {-6.25, -6.25, -33, 1.875};
> Point(1062) = {-6.25, -6.25, 21, 1.875};
> Point(1063) = {-6.25, 6.25, -33, 1.875};
> Point(1064) = {-6.25, 6.25, 21, 1.875};
> Point(1065) = {6.25, -6.25, -33, 1.875};
> Point(1066) = {6.25, -6.25, 21, 1.875};
> Point(1067) = {6.25, 6.25, -33, 1.875};
> Point(1068) = {6.25, 6.25, 21, 1.875};
> Point(1069) = {-15.625, -15.625, -46.5, 4.6875};
> Point(1070) = {-15.625, -15.625, 34.5, 4.6875};
> Point(1071) = {-15.625, 15.625, -46.5, 4.6875};
> Point(1072) = {-15.625, 15.625, 34.5, 4.6875};
> Point(1073) = {15.625, -15.625, -46.5, 4.6875};
> Point(1074) = {15.625, -15.625, 34.5, 4.6875};
> Point(1075) = {15.625, 15.625, -46.5, 4.6875};
> Point(1076) = {15.625, 15.625, 34.5, 4.6875};
> Ellipse (1) = {103, 101, 103, 104} Plane{0, 0, 1};
> Ellipse (2) = {106, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (3) = {102, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (4) = {107, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (5) = {105, 101, 106, 106} Plane{0, 0, 1};
> Ellipse (6) = {105, 101, 103, 103} Plane{0, 0, 1};
> Ellipse (7) = {105, 101, 107, 107} Plane{0, 0, 1};
> Ellipse (8) = {105, 101, 102, 102} Plane{0, 0, 1};
> Ellipse (9) = {103, 101, 106, 106} Plane{0, 0, 1};
> Ellipse (10) = {106, 101, 102, 102} Plane{0, 0, 1};
> Ellipse (11) = {102, 101, 107, 107} Plane{0, 0, 1};
> Ellipse (12) = {107, 101, 103, 103} Plane{0, 0, 1};
> Line (10031) = {201, 202};
> Line (10032) = {203, 204};
> Line (10033) = {207, 208};
> Line (10034) = {205, 206};
> Line (10035) = {205, 201};
> Line (10036) = {206, 202};
> Line (10037) = {208, 204};
> Line (10038) = {207, 203};
> Line (10039) = {203, 201};
> Line (10040) = {204, 202};
> Line (10041) = {208, 206};
> Line (10042) = {207, 205};
> Line (20057) = {304, 308};
> Line (20058) = {302, 306};
> Line (20059) = {301, 305};
> Line (20060) = {303, 307};
> Line (20061) = {303, 301};
> Line (20062) = {307, 305};
> Line (20063) = {304, 302};
> Line (20064) = {308, 306};
> Line (20065) = {304, 303};
> Line (20066) = {308, 307};
> Line (20067) = {302, 301};
> Line (20068) = {306, 305};
> Line (30057) = {404, 408};
> Line (30058) = {402, 406};
> Line (30059) = {401, 405};
> Line (30060) = {403, 407};
> Line (30061) = {403, 401};
> Line (30062) = {407, 405};
> Line (30063) = {404, 402};
> Line (30064) = {408, 406};
> Line (30065) = {404, 403};
> Line (30066) = {408, 407};
> Line (30067) = {402, 401};
> Line (30068) = {406, 405};
> Line (30082) = {1004, 1003};
> Line (30083) = {1002, 1001};
> Line (30084) = {1004, 1002};
> Line (30085) = {1003, 1001};
> Line (30088) = {1005, 1006};
> Line (30089) = {1007, 1008};
> Line (30090) = {1009, 1010};
> Line (30091) = {1011, 1012};
> Line (30092) = {1005, 1007};
> Line (30093) = {1006, 1008};
> Line (30094) = {1009, 1011};
> Line (30095) = {1010, 1012};
> Line (30096) = {1005, 1009};
> Line (30097) = {1006, 1010};
> Line (30098) = {1007, 1011};
> Line (30099) = {1008, 1012};
> Line (30114) = {1013, 1014};
> Line (30115) = {1015, 1016};
> Line (30116) = {1017, 1018};
> Line (30117) = {1019, 1020};
> Line (30118) = {1013, 1015};
> Line (30119) = {1014, 1016};
> Line (30120) = {1017, 1019};
> Line (30121) = {1018, 1020};
> Line (30122) = {1013, 1017};
> Line (30123) = {1014, 1018};
> Line (30124) = {1015, 1019};
> Line (30125) = {1016, 1020};
> Line (30140) = {1021, 1022};
> Line (30141) = {1023, 1024};
> Line (30142) = {1025, 1026};
> Line (30143) = {1027, 1028};
> Line (30144) = {1021, 1023};
> Line (30145) = {1022, 1024};
> Line (30146) = {1025, 1027};
> Line (30147) = {1026, 1028};
> Line (30148) = {1021, 1025};
> Line (30149) = {1022, 1026};
> Line (30150) = {1023, 1027};
> Line (30151) = {1024, 1028};
> Line (30166) = {1029, 1030};
> Line (30167) = {1031, 1032};
> Line (30168) = {1033, 1034};
> Line (30169) = {1035, 1036};
> Line (30170) = {1029, 1031};
> Line (30171) = {1030, 1032};
> Line (30172) = {1033, 1035};
> Line (30173) = {1034, 1036};
> Line (30174) = {1029, 1033};
> Line (30175) = {1030, 1034};
> Line (30176) = {1031, 1035};
> Line (30177) = {1032, 1036};
> Line (30192) = {1037, 1038};
> Line (30193) = {1039, 1040};
> Line (30194) = {1041, 1042};
> Line (30195) = {1043, 1044};
> Line (30196) = {1037, 1039};
> Line (30197) = {1038, 1040};
> Line (30198) = {1041, 1043};
> Line (30199) = {1042, 1044};
> Line (30200) = {1037, 1041};
> Line (30201) = {1038, 1042};
> Line (30202) = {1039, 1043};
> Line (30203) = {1040, 1044};
> Line (30218) = {1045, 1046};
> Line (30219) = {1047, 1048};
> Line (30220) = {1049, 1050};
> Line (30221) = {1051, 1052};
> Line (30222) = {1045, 1047};
> Line (30223) = {1046, 1048};
> Line (30224) = {1049, 1051};
> Line (30225) = {1050, 1052};
> Line (30226) = {1045, 1049};
> Line (30227) = {1046, 1050};
> Line (30228) = {1047, 1051};
> Line (30229) = {1048, 1052};
> Line (30244) = {1053, 1054};
> Line (30245) = {1055, 1056};
> Line (30246) = {1057, 1058};
> Line (30247) = {1059, 1060};
> Line (30248) = {1053, 1055};
> Line (30249) = {1054, 1056};
> Line (30250) = {1057, 1059};
> Line (30251) = {1058, 1060};
> Line (30252) = {1053, 1057};
> Line (30253) = {1054, 1058};
> Line (30254) = {1055, 1059};
> Line (30255) = {1056, 1060};
> Line (30269) = {1061, 1062};
> Line (30270) = {1063, 1064};
> Line (30271) = {1065, 1066};
> Line (30272) = {1067, 1068};
> Line (30273) = {1061, 1063};
> Line (30274) = {1062, 1064};
> Line (30275) = {1065, 1067};
> Line (30276) = {1066, 1068};
> Line (30277) = {1061, 1065};
> Line (30278) = {1062, 1066};
> Line (30279) = {1063, 1067};
> Line (30280) = {1064, 1068};
> Line (30294) = {1069, 1070};
> Line (30295) = {1071, 1072};
> Line (30296) = {1073, 1074};
> Line (30297) = {1075, 1076};
> Line (30298) = {1069, 1071};
> Line (30299) = {1070, 1072};
> Line (30300) = {1073, 1075};
> Line (30301) = {1074, 1076};
> Line (30302) = {1069, 1073};
> Line (30303) = {1070, 1074};
> Line (30304) = {1071, 1075};
> Line (30305) = {1072, 1076};
> Line Loop (1000014) = {2, -1, 9};
> Ruled Surface (14) = {1000014};
> Line Loop (1000016) = {4, -1, -12};
> Ruled Surface (16) = {1000016};
> Line Loop (1000018) = {2, -3, -10};
> Ruled Surface (18) = {1000018};
> Line Loop (1000020) = {3, -4, -11};
> Ruled Surface (20) = {1000020};
> Line Loop (1000022) = {5, -9, -6};
> Ruled Surface (22) = {1000022};
> Line Loop (1000024) = {7, 12, -6};
> Ruled Surface (24) = {1000024};
> Line Loop (1000026) = {7, -11, -8};
> Ruled Surface (26) = {1000026};
> Line Loop (1000028) = {8, -10, -5};
> Ruled Surface (28) = {1000028};
> Line Loop (1010044) = {10031, -10040, -10032, 10039};
> Plane Surface (10044) = {1010044};
> Line Loop (1010046) = {10042, 10034, -10041, -10033};
> Plane Surface (10046) = {1010046};
> Line Loop (1010048) = {10035, -10039, -10038, 10042};
> Plane Surface (10048) = {1010048};
> Line Loop (1010050) = {10036, -10040, -10037, 10041};
> Plane Surface (10050) = {1010050};
> Line Loop (1010052) = {10034, 10036, -10031, -10035};
> Plane Surface (10052) = {1010052};
> Line Loop (1010054) = {10032, -10037, -10033, 10038};
> Plane Surface (10054) = {1010054};
> Line Loop (1020070) = {20063, 20058, -20064, -20057};
> Plane Surface (20070) = {1020070};
> Line Loop (1020072) = {20061, 20059, -20062, -20060};
> Plane Surface (20072) = {1020072};
> Line Loop (1020074) = {20065, 20060, -20066, -20057};
> Plane Surface (20074) = {1020074};
> Line Loop (1020076) = {20067, 20059, -20068, -20058};
> Plane Surface (20076) = {1020076};
> Line Loop (1020078) = {20064, 20068, -20062, -20066};
> Plane Surface (20078) = {1020078};
> Line Loop (1030070) = {30063, 30058, -30064, -30057};
> Plane Surface (30070) = {1030070};
> Line Loop (1030072) = {30061, 30059, -30062, -30060};
> Plane Surface (30072) = {1030072};
> Line Loop (1030074) = {30065, 30060, -30066, -30057};
> Plane Surface (30074) = {1030074};
> Line Loop (1030076) = {30067, 30059, -30068, -30058};
> Plane Surface (30076) = {1030076};
> Line Loop (1030078) = {30064, 30068, -30062, -30066};
> Plane Surface (30078) = {1030078};
> Line Loop (1030080) = {30063, 30067, -30061, -30065};
> Plane Surface (30080) = {1030080};
> Line Loop (1030087) = {30084, 30083, -30085, -30082};
> Plane Surface (30087) = {1030087};
> Line Loop (1030090) = {20063, 20067, -20061, -20065, 30222, 30219, -30223, -30218};
> Plane Surface (30090) = {1030090};
> Line Loop (1030106) = {30088, 30093, -30089, -30092, 30082, 30085, -30083, -30084};
> Plane Surface (30106) = {1030106};
> Line Loop (1030107) = {-30090, 30094, 30091, -30095};
> Plane Surface (30107) = {1030107};
> Line Loop (1030108) = {-30088, 30096, 30090, -30097};
> Plane Surface (30108) = {1030108};
> Line Loop (1030109) = {30089, 30099, -30091, -30098};
> Plane Surface (30109) = {1030109};
> Line Loop (1030110) = {30092, 30098, -30094, -30096};
> Plane Surface (30110) = {1030110};
> Line Loop (1030111) = {-30093, 30097, 30095, -30099};
> Plane Surface (30111) = {1030111};
> Line Loop (1030132) = {30114, 30119, -30115, -30118, 30092, 30089, -30093, -30088};
> Plane Surface (30132) = {1030132};
> Line Loop (1030133) = {-30116, 30120, 30117, -30121};
> Plane Surface (30133) = {1030133};
> Line Loop (1030134) = {-30114, 30122, 30116, -30123};
> Plane Surface (30134) = {1030134};
> Line Loop (1030135) = {30115, 30125, -30117, -30124};
> Plane Surface (30135) = {1030135};
> Line Loop (1030136) = {30118, 30124, -30120, -30122};
> Plane Surface (30136) = {1030136};
> Line Loop (1030137) = {-30119, 30123, 30121, -30125};
> Plane Surface (30137) = {1030137};
> Line Loop (1030158) = {30140, 30145, -30141, -30144, 30118, 30115, -30119, -30114};
> Plane Surface (30158) = {1030158};
> Line Loop (1030159) = {-30142, 30146, 30143, -30147};
> Plane Surface (30159) = {1030159};
> Line Loop (1030160) = {-30140, 30148, 30142, -30149};
> Plane Surface (30160) = {1030160};
> Line Loop (1030161) = {30141, 30151, -30143, -30150};
> Plane Surface (30161) = {1030161};
> Line Loop (1030162) = {30144, 30150, -30146, -30148};
> Plane Surface (30162) = {1030162};
> Line Loop (1030163) = {-30145, 30149, 30147, -30151};
> Plane Surface (30163) = {1030163};
> Line Loop (1030184) = {30166, 30171, -30167, -30170, 30144, 30141, -30145, -30140};
> Plane Surface (30184) = {1030184};
> Line Loop (1030185) = {-30168, 30172, 30169, -30173};
> Plane Surface (30185) = {1030185};
> Line Loop (1030186) = {-30166, 30174, 30168, -30175};
> Plane Surface (30186) = {1030186};
> Line Loop (1030187) = {30167, 30177, -30169, -30176};
> Plane Surface (30187) = {1030187};
> Line Loop (1030188) = {30170, 30176, -30172, -30174};
> Plane Surface (30188) = {1030188};
> Line Loop (1030189) = {-30171, 30175, 30173, -30177};
> Plane Surface (30189) = {1030189};
> Line Loop (1030210) = {30192, 30197, -30193, -30196, 30170, 30167, -30171, -30166};
> Plane Surface (30210) = {1030210};
> Line Loop (1030211) = {-30194, 30198, 30195, -30199};
> Plane Surface (30211) = {1030211};
> Line Loop (1030212) = {-30192, 30200, 30194, -30201};
> Plane Surface (30212) = {1030212};
> Line Loop (1030213) = {30193, 30203, -30195, -30202};
> Plane Surface (30213) = {1030213};
> Line Loop (1030214) = {30196, 30202, -30198, -30200};
> Plane Surface (30214) = {1030214};
> Line Loop (1030215) = {-30197, 30201, 30199, -30203};
> Plane Surface (30215) = {1030215};
> Line Loop (1030236) = {30218, 30223, -30219, -30222, 30196, 30193, -30197, -30192};
> Plane Surface (30236) = {1030236};
> Line Loop (1030237) = {-30220, 30224, 30221, -30225};
> Plane Surface (30237) = {1030237};
> Line Loop (1030238) = {-30218, 30226, 30220, -30227};
> Plane Surface (30238) = {1030238};
> Line Loop (1030239) = {30219, 30229, -30221, -30228};
> Plane Surface (30239) = {1030239};
> Line Loop (1030240) = {30222, 30228, -30224, -30226};
> Plane Surface (30240) = {1030240};
> Line Loop (1030241) = {-30223, 30227, 30225, -30229};
> Plane Surface (30241) = {1030241};
> Line Loop (1030257) = {30244, 30249, -30245, -30248};
> Plane Surface (30257) = {1030257};
> Line Loop (1030259) = {-30246, 30250, 30247, -30251};
> Plane Surface (30259) = {1030259};
> Line Loop (1030261) = {-30244, 30252, 30246, -30253};
> Plane Surface (30261) = {1030261};
> Line Loop (1030263) = {30245, 30255, -30247, -30254};
> Plane Surface (30263) = {1030263};
> Line Loop (1030265) = {30248, 30254, -30250, -30252};
> Plane Surface (30265) = {1030265};
> Line Loop (1030267) = {-30249, 30253, 30251, -30255};
> Plane Surface (30267) = {1030267};
> Line Loop (1030282) = {30269, 30274, -30270, -30273};
> Plane Surface (30282) = {1030282};
> Line Loop (1030284) = {-30271, 30275, 30272, -30276};
> Plane Surface (30284) = {1030284};
> Line Loop (1030286) = {-30269, 30277, 30271, -30278};
> Plane Surface (30286) = {1030286};
> Line Loop (1030288) = {30270, 30280, -30272, -30279};
> Plane Surface (30288) = {1030288};
> Line Loop (1030290) = {30273, 30279, -30275, -30277};
> Plane Surface (30290) = {1030290};
> Line Loop (1030292) = {-30274, 30278, 30276, -30280};
> Plane Surface (30292) = {1030292};
> Line Loop (1030307) = {30294, 30299, -30295, -30298};
> Plane Surface (30307) = {1030307};
> Line Loop (1030309) = {-30296, 30300, 30297, -30301};
> Plane Surface (30309) = {1030309};
> Line Loop (1030311) = {-30294, 30302, 30296, -30303};
> Plane Surface (30311) = {1030311};
> Line Loop (1030313) = {30295, 30305, -30297, -30304};
> Plane Surface (30313) = {1030313};
> Line Loop (1030315) = {30298, 30304, -30300, -30302};
> Plane Surface (30315) = {1030315};
> Line Loop (1030317) = {-30299, 30303, 30301, -30305};
> Plane Surface (30317) = {1030317};
> Surface Loop (1000001) = {30282, 30284, 30286, 30288, 30290, 30292, 30257, 30259, 30261, 30263, 30265, 30267};
> Volume (1) = {1000001};
> Surface Loop (1000002) = {30307, 30309, 30311, 30313, 30315, 30317, 30282, 30284, 30286, 30288, 30290, 30292};
> Volume (2) = {1000002};
> Surface Loop (1000003) = {22, 28, 26, 24, 16, 20, 18, 14, 30307, 30309, 30311, 30313, 30315, 30317};
> Volume (3) = {1000003};
> Surface Loop (1001000) = {30257, 30259, 30261, 30263, 30265, 30267, 10052, 10046, 10048, 10044, 10050, 10054, 30090, 30237, 30238, 30239, 30240, 30241, 20070, 20076, 20072, 20078, 20074, 30080, 30070, 30076, 30072, 30078, 30074};
> Volume (1000) = {1001000};
> Surface Loop (1030113) = {30106, 30107, 30108, 30109, 30110, 30111, 30087};
> Volume (30113) = {1030113};
> Surface Loop (1030139) = {30132, 30133, 30134, 30135, 30136, 30137, 30107, 30108, 30109, 30110, 30111};
> Volume (30139) = {1030139};
> Surface Loop (1030165) = {30158, 30159, 30160, 30161, 30162, 30163, 30133, 30134, 30135, 30136, 30137};
> Volume (30165) = {1030165};
> Surface Loop (1030191) = {30184, 30185, 30186, 30187, 30188, 30189, 30159, 30160, 30161, 30162, 30163};
> Volume (30191) = {1030191};
> Surface Loop (1030217) = {30210, 30211, 30212, 30213, 30214, 30215, 30185, 30186, 30187, 30188, 30189};
> Volume (30217) = {1030217};
> Surface Loop (1030243) = {30236, 30237, 30238, 30239, 30240, 30241, 30211, 30212, 30213, 30214, 30215};
> Volume (30243) = {1030243};
> Physical Surface (1) = {14, 16, 18, 20, 22, 24, 26, 28};
> Physical Surface (2) = {10044, 10046, 10048, 10050, 10052, 10054};
> Physical Surface (3) = {30087};
> Physical Surface (5) = {30070, 30072, 30074, 30076, 30078, 30080};
> 
> 
> ------------------------------------------------------------------------
> 
> Point(101) = {0, 0, -6, 20};
> Point(102) = {100, 0, -6, 20};
> Point(103) = {-100, 0, -6, 20};
> Point(104) = {0, 100, -6, 20};
> Point(105) = {0, -100, -6, 20};
> Point(106) = {0, 0, 194, 20};
> Point(107) = {0, 0, -206, 20};
> Point(201) = {-1, -1, -7, 0.3};
> Point(202) = {-1, -1, -5, 0.3};
> Point(203) = {-1, 1, -7, 0.3};
> Point(204) = {-1, 1, -5, 0.3};
> Point(205) = {1, -1, -7, 0.3};
> Point(206) = {1, -1, -5, 0.3};
> Point(207) = {1, 1, -7, 0.3};
> Point(208) = {1, 1, -5, 0.3};
> Point(301) = {0, -1, -4, 0.4};
> Point(302) = {0, -1, 6, 0.4};
> Point(303) = {0, 1, -4, 0.4};
> Point(304) = {0, 1, 6, 0.4};
> Point(305) = {-0.2, -1, -4, 0.4};
> Point(306) = {-0.2, -1, 6, 0.4};
> Point(307) = {-0.2, 1, -4, 0.4};
> Point(308) = {-0.2, 1, 6, 0.4};
> Point(401) = {0, -1, -8, 0.6};
> Point(402) = {0, -1, -18, 0.6};
> Point(403) = {0, 1, -8, 0.6};
> Point(404) = {0, 1, -18, 0.6};
> Point(405) = {-0.2, -1, -8, 0.6};
> Point(406) = {-0.2, -1, -18, 0.6};
> Point(407) = {-0.2, 1, -8, 0.6};
> Point(408) = {-0.2, 1, -18, 0.6};
> Point(1001) = {0, -0.01, 0.9995000000000001, 0.0007};
> Point(1002) = {0, -0.01, 1.0005, 0.0007};
> Point(1003) = {0, 0.01, 0.9995000000000001, 0.0007};
> Point(1004) = {0, 0.01, 1.0005, 0.0007};
> Point(1005) = {0, -0.02, 0.997, 0.00175};
> Point(1006) = {0, -0.02, 1.003, 0.00175};
> Point(1007) = {0, 0.02, 0.997, 0.00175};
> Point(1008) = {0, 0.02, 1.003, 0.00175};
> Point(1009) = {0.005590169943749475, -0.02, 0.997, 0.002625};
> Point(1010) = {0.005590169943749475, -0.02, 1.003, 0.002625};
> Point(1011) = {0.005590169943749475, 0.02, 0.997, 0.002625};
> Point(1012) = {0.005590169943749475, 0.02, 1.003, 0.002625};
> Point(1013) = {0, -0.04, 0.991, 0.004375};
> Point(1014) = {0, -0.04, 1.009, 0.004375};
> Point(1015) = {0, 0.04, 0.991, 0.004375};
> Point(1016) = {0, 0.04, 1.009, 0.004375};
> Point(1017) = {0.01397542485937369, -0.04, 0.991, 0.006562500000000001};
> Point(1018) = {0.01397542485937369, -0.04, 1.009, 0.006562500000000001};
> Point(1019) = {0.01397542485937369, 0.04, 0.991, 0.006562500000000001};
> Point(1020) = {0.01397542485937369, 0.04, 1.009, 0.006562500000000001};
> Point(1021) = {0, -0.08, 0.973, 0.0109375};
> Point(1022) = {0, -0.08, 1.027, 0.0109375};
> Point(1023) = {0, 0.08, 0.973, 0.0109375};
> Point(1024) = {0, 0.08, 1.027, 0.0109375};
> Point(1025) = {0.03493856214843422, -0.08, 0.973, 0.01640625};
> Point(1026) = {0.03493856214843422, -0.08, 1.027, 0.01640625};
> Point(1027) = {0.03493856214843422, 0.08, 0.973, 0.01640625};
> Point(1028) = {0.03493856214843422, 0.08, 1.027, 0.01640625};
> Point(1029) = {0, -0.16, 0.919, 0.02734375};
> Point(1030) = {0, -0.16, 1.081, 0.02734375};
> Point(1031) = {0, 0.16, 0.919, 0.02734375};
> Point(1032) = {0, 0.16, 1.081, 0.02734375};
> Point(1033) = {0.08734640537108554, -0.16, 0.919, 0.04101562500000001};
> Point(1034) = {0.08734640537108554, -0.16, 1.081, 0.04101562500000001};
> Point(1035) = {0.08734640537108554, 0.16, 0.919, 0.04101562500000001};
> Point(1036) = {0.08734640537108554, 0.16, 1.081, 0.04101562500000001};
> Point(1037) = {0, -0.32, 0.7569999999999999, 0.06835937500000001};
> Point(1038) = {0, -0.32, 1.243, 0.06835937500000001};
> Point(1039) = {0, 0.32, 0.7569999999999999, 0.06835937500000001};
> Point(1040) = {0, 0.32, 1.243, 0.06835937500000001};
> Point(1041) = {0.2183660134277138, -0.32, 0.7569999999999999, 0.1025390625};
> Point(1042) = {0.2183660134277138, -0.32, 1.243, 0.1025390625};
> Point(1043) = {0.2183660134277138, 0.32, 0.7569999999999999, 0.1025390625};
> Point(1044) = {0.2183660134277138, 0.32, 1.243, 0.1025390625};
> Point(1045) = {0, -0.64, 0.2709999999999999, 0.1708984375};
> Point(1046) = {0, -0.64, 1.729, 0.1708984375};
> Point(1047) = {0, 0.64, 0.2709999999999999, 0.1708984375};
> Point(1048) = {0, 0.64, 1.729, 0.1708984375};
> Point(1049) = {0.5459150335692846, -0.64, 0.2709999999999999, 0.2563476562500001};
> Point(1050) = {0.5459150335692846, -0.64, 1.729, 0.2563476562500001};
> Point(1051) = {0.5459150335692846, 0.64, 0.2709999999999999, 0.2563476562500001};
> Point(1052) = {0.5459150335692846, 0.64, 1.729, 0.2563476562500001};
> Point(1053) = {-2.5, -2.5, -24, 0.75};
> Point(1054) = {-2.5, -2.5, 12, 0.75};
> Point(1055) = {-2.5, 2.5, -24, 0.75};
> Point(1056) = {-2.5, 2.5, 12, 0.75};
> Point(1057) = {2.5, -2.5, -24, 0.75};
> Point(1058) = {2.5, -2.5, 12, 0.75};
> Point(1059) = {2.5, 2.5, -24, 0.75};
> Point(1060) = {2.5, 2.5, 12, 0.75};
> Point(1061) = {-6.25, -6.25, -33, 1.875};
> Point(1062) = {-6.25, -6.25, 21, 1.875};
> Point(1063) = {-6.25, 6.25, -33, 1.875};
> Point(1064) = {-6.25, 6.25, 21, 1.875};
> Point(1065) = {6.25, -6.25, -33, 1.875};
> Point(1066) = {6.25, -6.25, 21, 1.875};
> Point(1067) = {6.25, 6.25, -33, 1.875};
> Point(1068) = {6.25, 6.25, 21, 1.875};
> Point(1069) = {-15.625, -15.625, -46.5, 4.6875};
> Point(1070) = {-15.625, -15.625, 34.5, 4.6875};
> Point(1071) = {-15.625, 15.625, -46.5, 4.6875};
> Point(1072) = {-15.625, 15.625, 34.5, 4.6875};
> Point(1073) = {15.625, -15.625, -46.5, 4.6875};
> Point(1074) = {15.625, -15.625, 34.5, 4.6875};
> Point(1075) = {15.625, 15.625, -46.5, 4.6875};
> Point(1076) = {15.625, 15.625, 34.5, 4.6875};
> Point(1077) = {-39.0625, -39.0625, -66.75, 11.71875};
> Point(1078) = {-39.0625, -39.0625, 54.75, 11.71875};
> Point(1079) = {-39.0625, 39.0625, -66.75, 11.71875};
> Point(1080) = {-39.0625, 39.0625, 54.75, 11.71875};
> Point(1081) = {39.0625, -39.0625, -66.75, 11.71875};
> Point(1082) = {39.0625, -39.0625, 54.75, 11.71875};
> Point(1083) = {39.0625, 39.0625, -66.75, 11.71875};
> Point(1084) = {39.0625, 39.0625, 54.75, 11.71875};
> Ellipse (1) = {103, 101, 103, 104} Plane{0, 0, 1};
> Ellipse (2) = {106, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (3) = {102, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (4) = {107, 101, 104, 104} Plane{0, 0, 1};
> Ellipse (5) = {105, 101, 106, 106} Plane{0, 0, 1};
> Ellipse (6) = {105, 101, 103, 103} Plane{0, 0, 1};
> Ellipse (7) = {105, 101, 107, 107} Plane{0, 0, 1};
> Ellipse (8) = {105, 101, 102, 102} Plane{0, 0, 1};
> Ellipse (9) = {103, 101, 106, 106} Plane{0, 0, 1};
> Ellipse (10) = {106, 101, 102, 102} Plane{0, 0, 1};
> Ellipse (11) = {102, 101, 107, 107} Plane{0, 0, 1};
> Ellipse (12) = {107, 101, 103, 103} Plane{0, 0, 1};
> Line (10031) = {201, 202};
> Line (10032) = {203, 204};
> Line (10033) = {207, 208};
> Line (10034) = {205, 206};
> Line (10035) = {205, 201};
> Line (10036) = {206, 202};
> Line (10037) = {208, 204};
> Line (10038) = {207, 203};
> Line (10039) = {203, 201};
> Line (10040) = {204, 202};
> Line (10041) = {208, 206};
> Line (10042) = {207, 205};
> Line (20057) = {304, 308};
> Line (20058) = {302, 306};
> Line (20059) = {301, 305};
> Line (20060) = {303, 307};
> Line (20061) = {303, 301};
> Line (20062) = {307, 305};
> Line (20063) = {304, 302};
> Line (20064) = {308, 306};
> Line (20065) = {304, 303};
> Line (20066) = {308, 307};
> Line (20067) = {302, 301};
> Line (20068) = {306, 305};
> Line (30057) = {404, 408};
> Line (30058) = {402, 406};
> Line (30059) = {401, 405};
> Line (30060) = {403, 407};
> Line (30061) = {403, 401};
> Line (30062) = {407, 405};
> Line (30063) = {404, 402};
> Line (30064) = {408, 406};
> Line (30065) = {404, 403};
> Line (30066) = {408, 407};
> Line (30067) = {402, 401};
> Line (30068) = {406, 405};
> Line (30082) = {1004, 1003};
> Line (30083) = {1002, 1001};
> Line (30084) = {1004, 1002};
> Line (30085) = {1003, 1001};
> Line (30088) = {1005, 1006};
> Line (30089) = {1007, 1008};
> Line (30090) = {1009, 1010};
> Line (30091) = {1011, 1012};
> Line (30092) = {1005, 1007};
> Line (30093) = {1006, 1008};
> Line (30094) = {1009, 1011};
> Line (30095) = {1010, 1012};
> Line (30096) = {1005, 1009};
> Line (30097) = {1006, 1010};
> Line (30098) = {1007, 1011};
> Line (30099) = {1008, 1012};
> Line (30114) = {1013, 1014};
> Line (30115) = {1015, 1016};
> Line (30116) = {1017, 1018};
> Line (30117) = {1019, 1020};
> Line (30118) = {1013, 1015};
> Line (30119) = {1014, 1016};
> Line (30120) = {1017, 1019};
> Line (30121) = {1018, 1020};
> Line (30122) = {1013, 1017};
> Line (30123) = {1014, 1018};
> Line (30124) = {1015, 1019};
> Line (30125) = {1016, 1020};
> Line (30140) = {1021, 1022};
> Line (30141) = {1023, 1024};
> Line (30142) = {1025, 1026};
> Line (30143) = {1027, 1028};
> Line (30144) = {1021, 1023};
> Line (30145) = {1022, 1024};
> Line (30146) = {1025, 1027};
> Line (30147) = {1026, 1028};
> Line (30148) = {1021, 1025};
> Line (30149) = {1022, 1026};
> Line (30150) = {1023, 1027};
> Line (30151) = {1024, 1028};
> Line (30166) = {1029, 1030};
> Line (30167) = {1031, 1032};
> Line (30168) = {1033, 1034};
> Line (30169) = {1035, 1036};
> Line (30170) = {1029, 1031};
> Line (30171) = {1030, 1032};
> Line (30172) = {1033, 1035};
> Line (30173) = {1034, 1036};
> Line (30174) = {1029, 1033};
> Line (30175) = {1030, 1034};
> Line (30176) = {1031, 1035};
> Line (30177) = {1032, 1036};
> Line (30192) = {1037, 1038};
> Line (30193) = {1039, 1040};
> Line (30194) = {1041, 1042};
> Line (30195) = {1043, 1044};
> Line (30196) = {1037, 1039};
> Line (30197) = {1038, 1040};
> Line (30198) = {1041, 1043};
> Line (30199) = {1042, 1044};
> Line (30200) = {1037, 1041};
> Line (30201) = {1038, 1042};
> Line (30202) = {1039, 1043};
> Line (30203) = {1040, 1044};
> Line (30218) = {1045, 1046};
> Line (30219) = {1047, 1048};
> Line (30220) = {1049, 1050};
> Line (30221) = {1051, 1052};
> Line (30222) = {1045, 1047};
> Line (30223) = {1046, 1048};
> Line (30224) = {1049, 1051};
> Line (30225) = {1050, 1052};
> Line (30226) = {1045, 1049};
> Line (30227) = {1046, 1050};
> Line (30228) = {1047, 1051};
> Line (30229) = {1048, 1052};
> Line (30244) = {1053, 1054};
> Line (30245) = {1055, 1056};
> Line (30246) = {1057, 1058};
> Line (30247) = {1059, 1060};
> Line (30248) = {1053, 1055};
> Line (30249) = {1054, 1056};
> Line (30250) = {1057, 1059};
> Line (30251) = {1058, 1060};
> Line (30252) = {1053, 1057};
> Line (30253) = {1054, 1058};
> Line (30254) = {1055, 1059};
> Line (30255) = {1056, 1060};
> Line (30269) = {1061, 1062};
> Line (30270) = {1063, 1064};
> Line (30271) = {1065, 1066};
> Line (30272) = {1067, 1068};
> Line (30273) = {1061, 1063};
> Line (30274) = {1062, 1064};
> Line (30275) = {1065, 1067};
> Line (30276) = {1066, 1068};
> Line (30277) = {1061, 1065};
> Line (30278) = {1062, 1066};
> Line (30279) = {1063, 1067};
> Line (30280) = {1064, 1068};
> Line (30294) = {1069, 1070};
> Line (30295) = {1071, 1072};
> Line (30296) = {1073, 1074};
> Line (30297) = {1075, 1076};
> Line (30298) = {1069, 1071};
> Line (30299) = {1070, 1072};
> Line (30300) = {1073, 1075};
> Line (30301) = {1074, 1076};
> Line (30302) = {1069, 1073};
> Line (30303) = {1070, 1074};
> Line (30304) = {1071, 1075};
> Line (30305) = {1072, 1076};
> Line (30319) = {1077, 1078};
> Line (30320) = {1079, 1080};
> Line (30321) = {1081, 1082};
> Line (30322) = {1083, 1084};
> Line (30323) = {1077, 1079};
> Line (30324) = {1078, 1080};
> Line (30325) = {1081, 1083};
> Line (30326) = {1082, 1084};
> Line (30327) = {1077, 1081};
> Line (30328) = {1078, 1082};
> Line (30329) = {1079, 1083};
> Line (30330) = {1080, 1084};
> Line Loop (1000014) = {2, -1, 9};
> Ruled Surface (14) = {1000014};
> Line Loop (1000016) = {4, -1, -12};
> Ruled Surface (16) = {1000016};
> Line Loop (1000018) = {2, -3, -10};
> Ruled Surface (18) = {1000018};
> Line Loop (1000020) = {3, -4, -11};
> Ruled Surface (20) = {1000020};
> Line Loop (1000022) = {5, -9, -6};
> Ruled Surface (22) = {1000022};
> Line Loop (1000024) = {7, 12, -6};
> Ruled Surface (24) = {1000024};
> Line Loop (1000026) = {7, -11, -8};
> Ruled Surface (26) = {1000026};
> Line Loop (1000028) = {8, -10, -5};
> Ruled Surface (28) = {1000028};
> Line Loop (1010044) = {10031, -10040, -10032, 10039};
> Plane Surface (10044) = {1010044};
> Line Loop (1010046) = {10042, 10034, -10041, -10033};
> Plane Surface (10046) = {1010046};
> Line Loop (1010048) = {10035, -10039, -10038, 10042};
> Plane Surface (10048) = {1010048};
> Line Loop (1010050) = {10036, -10040, -10037, 10041};
> Plane Surface (10050) = {1010050};
> Line Loop (1010052) = {10034, 10036, -10031, -10035};
> Plane Surface (10052) = {1010052};
> Line Loop (1010054) = {10032, -10037, -10033, 10038};
> Plane Surface (10054) = {1010054};
> Line Loop (1020070) = {20063, 20058, -20064, -20057};
> Plane Surface (20070) = {1020070};
> Line Loop (1020072) = {20061, 20059, -20062, -20060};
> Plane Surface (20072) = {1020072};
> Line Loop (1020074) = {20065, 20060, -20066, -20057};
> Plane Surface (20074) = {1020074};
> Line Loop (1020076) = {20067, 20059, -20068, -20058};
> Plane Surface (20076) = {1020076};
> Line Loop (1020078) = {20064, 20068, -20062, -20066};
> Plane Surface (20078) = {1020078};
> Line Loop (1030070) = {30063, 30058, -30064, -30057};
> Plane Surface (30070) = {1030070};
> Line Loop (1030072) = {30061, 30059, -30062, -30060};
> Plane Surface (30072) = {1030072};
> Line Loop (1030074) = {30065, 30060, -30066, -30057};
> Plane Surface (30074) = {1030074};
> Line Loop (1030076) = {30067, 30059, -30068, -30058};
> Plane Surface (30076) = {1030076};
> Line Loop (1030078) = {30064, 30068, -30062, -30066};
> Plane Surface (30078) = {1030078};
> Line Loop (1030080) = {30063, 30067, -30061, -30065};
> Plane Surface (30080) = {1030080};
> Line Loop (1030087) = {30084, 30083, -30085, -30082};
> Plane Surface (30087) = {1030087};
> Line Loop (1030090) = {20063, 20067, -20061, -20065, 30222, 30219, -30223, -30218};
> Plane Surface (30090) = {1030090};
> Line Loop (1030106) = {30088, 30093, -30089, -30092, 30082, 30085, -30083, -30084};
> Plane Surface (30106) = {1030106};
> Line Loop (1030107) = {-30090, 30094, 30091, -30095};
> Plane Surface (30107) = {1030107};
> Line Loop (1030108) = {-30088, 30096, 30090, -30097};
> Plane Surface (30108) = {1030108};
> Line Loop (1030109) = {30089, 30099, -30091, -30098};
> Plane Surface (30109) = {1030109};
> Line Loop (1030110) = {30092, 30098, -30094, -30096};
> Plane Surface (30110) = {1030110};
> Line Loop (1030111) = {-30093, 30097, 30095, -30099};
> Plane Surface (30111) = {1030111};
> Line Loop (1030132) = {30114, 30119, -30115, -30118, 30092, 30089, -30093, -30088};
> Plane Surface (30132) = {1030132};
> Line Loop (1030133) = {-30116, 30120, 30117, -30121};
> Plane Surface (30133) = {1030133};
> Line Loop (1030134) = {-30114, 30122, 30116, -30123};
> Plane Surface (30134) = {1030134};
> Line Loop (1030135) = {30115, 30125, -30117, -30124};
> Plane Surface (30135) = {1030135};
> Line Loop (1030136) = {30118, 30124, -30120, -30122};
> Plane Surface (30136) = {1030136};
> Line Loop (1030137) = {-30119, 30123, 30121, -30125};
> Plane Surface (30137) = {1030137};
> Line Loop (1030158) = {30140, 30145, -30141, -30144, 30118, 30115, -30119, -30114};
> Plane Surface (30158) = {1030158};
> Line Loop (1030159) = {-30142, 30146, 30143, -30147};
> Plane Surface (30159) = {1030159};
> Line Loop (1030160) = {-30140, 30148, 30142, -30149};
> Plane Surface (30160) = {1030160};
> Line Loop (1030161) = {30141, 30151, -30143, -30150};
> Plane Surface (30161) = {1030161};
> Line Loop (1030162) = {30144, 30150, -30146, -30148};
> Plane Surface (30162) = {1030162};
> Line Loop (1030163) = {-30145, 30149, 30147, -30151};
> Plane Surface (30163) = {1030163};
> Line Loop (1030184) = {30166, 30171, -30167, -30170, 30144, 30141, -30145, -30140};
> Plane Surface (30184) = {1030184};
> Line Loop (1030185) = {-30168, 30172, 30169, -30173};
> Plane Surface (30185) = {1030185};
> Line Loop (1030186) = {-30166, 30174, 30168, -30175};
> Plane Surface (30186) = {1030186};
> Line Loop (1030187) = {30167, 30177, -30169, -30176};
> Plane Surface (30187) = {1030187};
> Line Loop (1030188) = {30170, 30176, -30172, -30174};
> Plane Surface (30188) = {1030188};
> Line Loop (1030189) = {-30171, 30175, 30173, -30177};
> Plane Surface (30189) = {1030189};
> Line Loop (1030210) = {30192, 30197, -30193, -30196, 30170, 30167, -30171, -30166};
> Plane Surface (30210) = {1030210};
> Line Loop (1030211) = {-30194, 30198, 30195, -30199};
> Plane Surface (30211) = {1030211};
> Line Loop (1030212) = {-30192, 30200, 30194, -30201};
> Plane Surface (30212) = {1030212};
> Line Loop (1030213) = {30193, 30203, -30195, -30202};
> Plane Surface (30213) = {1030213};
> Line Loop (1030214) = {30196, 30202, -30198, -30200};
> Plane Surface (30214) = {1030214};
> Line Loop (1030215) = {-30197, 30201, 30199, -30203};
> Plane Surface (30215) = {1030215};
> Line Loop (1030236) = {30218, 30223, -30219, -30222, 30196, 30193, -30197, -30192};
> Plane Surface (30236) = {1030236};
> Line Loop (1030237) = {-30220, 30224, 30221, -30225};
> Plane Surface (30237) = {1030237};
> Line Loop (1030238) = {-30218, 30226, 30220, -30227};
> Plane Surface (30238) = {1030238};
> Line Loop (1030239) = {30219, 30229, -30221, -30228};
> Plane Surface (30239) = {1030239};
> Line Loop (1030240) = {30222, 30228, -30224, -30226};
> Plane Surface (30240) = {1030240};
> Line Loop (1030241) = {-30223, 30227, 30225, -30229};
> Plane Surface (30241) = {1030241};
> Line Loop (1030257) = {30244, 30249, -30245, -30248};
> Plane Surface (30257) = {1030257};
> Line Loop (1030259) = {-30246, 30250, 30247, -30251};
> Plane Surface (30259) = {1030259};
> Line Loop (1030261) = {-30244, 30252, 30246, -30253};
> Plane Surface (30261) = {1030261};
> Line Loop (1030263) = {30245, 30255, -30247, -30254};
> Plane Surface (30263) = {1030263};
> Line Loop (1030265) = {30248, 30254, -30250, -30252};
> Plane Surface (30265) = {1030265};
> Line Loop (1030267) = {-30249, 30253, 30251, -30255};
> Plane Surface (30267) = {1030267};
> Line Loop (1030282) = {30269, 30274, -30270, -30273};
> Plane Surface (30282) = {1030282};
> Line Loop (1030284) = {-30271, 30275, 30272, -30276};
> Plane Surface (30284) = {1030284};
> Line Loop (1030286) = {-30269, 30277, 30271, -30278};
> Plane Surface (30286) = {1030286};
> Line Loop (1030288) = {30270, 30280, -30272, -30279};
> Plane Surface (30288) = {1030288};
> Line Loop (1030290) = {30273, 30279, -30275, -30277};
> Plane Surface (30290) = {1030290};
> Line Loop (1030292) = {-30274, 30278, 30276, -30280};
> Plane Surface (30292) = {1030292};
> Line Loop (1030307) = {30294, 30299, -30295, -30298};
> Plane Surface (30307) = {1030307};
> Line Loop (1030309) = {-30296, 30300, 30297, -30301};
> Plane Surface (30309) = {1030309};
> Line Loop (1030311) = {-30294, 30302, 30296, -30303};
> Plane Surface (30311) = {1030311};
> Line Loop (1030313) = {30295, 30305, -30297, -30304};
> Plane Surface (30313) = {1030313};
> Line Loop (1030315) = {30298, 30304, -30300, -30302};
> Plane Surface (30315) = {1030315};
> Line Loop (1030317) = {-30299, 30303, 30301, -30305};
> Plane Surface (30317) = {1030317};
> Line Loop (1030332) = {30319, 30324, -30320, -30323};
> Plane Surface (30332) = {1030332};
> Line Loop (1030334) = {-30321, 30325, 30322, -30326};
> Plane Surface (30334) = {1030334};
> Line Loop (1030336) = {-30319, 30327, 30321, -30328};
> Plane Surface (30336) = {1030336};
> Line Loop (1030338) = {30320, 30330, -30322, -30329};
> Plane Surface (30338) = {1030338};
> Line Loop (1030340) = {30323, 30329, -30325, -30327};
> Plane Surface (30340) = {1030340};
> Line Loop (1030342) = {-30324, 30328, 30326, -30330};
> Plane Surface (30342) = {1030342};
> Surface Loop (1000001) = {30282, 30284, 30286, 30288, 30290, 30292, 30257, 30259, 30261, 30263, 30265, 30267};
> Volume (1) = {1000001};
> Surface Loop (1000002) = {30307, 30309, 30311, 30313, 30315, 30317, 30282, 30284, 30286, 30288, 30290, 30292};
> Volume (2) = {1000002};
> Surface Loop (1000003) = {30332, 30334, 30336, 30338, 30340, 30342, 30307, 30309, 30311, 30313, 30315, 30317};
> Volume (3) = {1000003};
> Surface Loop (1000004) = {22, 28, 26, 24, 16, 20, 18, 14, 30332, 30334, 30336, 30338, 30340, 30342};
> Volume (4) = {1000004};
> Surface Loop (1001000) = {30257, 30259, 30261, 30263, 30265, 30267, 10052, 10046, 10048, 10044, 10050, 10054, 30090, 30237, 30238, 30239, 30240, 30241, 20070, 20076, 20072, 20078, 20074, 30080, 30070, 30076, 30072, 30078, 30074};
> Volume (1000) = {1001000};
> Surface Loop (1030113) = {30106, 30107, 30108, 30109, 30110, 30111, 30087};
> Volume (30113) = {1030113};
> Surface Loop (1030139) = {30132, 30133, 30134, 30135, 30136, 30137, 30107, 30108, 30109, 30110, 30111};
> Volume (30139) = {1030139};
> Surface Loop (1030165) = {30158, 30159, 30160, 30161, 30162, 30163, 30133, 30134, 30135, 30136, 30137};
> Volume (30165) = {1030165};
> Surface Loop (1030191) = {30184, 30185, 30186, 30187, 30188, 30189, 30159, 30160, 30161, 30162, 30163};
> Volume (30191) = {1030191};
> Surface Loop (1030217) = {30210, 30211, 30212, 30213, 30214, 30215, 30185, 30186, 30187, 30188, 30189};
> Volume (30217) = {1030217};
> Surface Loop (1030243) = {30236, 30237, 30238, 30239, 30240, 30241, 30211, 30212, 30213, 30214, 30215};
> Volume (30243) = {1030243};
> Physical Surface (1) = {14, 16, 18, 20, 22, 24, 26, 28};
> Physical Surface (2) = {10044, 10046, 10048, 10050, 10052, 10054};
> Physical Surface (3) = {30087};
> Physical Surface (5) = {30070, 30072, 30074, 30076, 30078, 30080};
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Christophe Geuzaine
Applied and Computational Mathematics, Caltech
geuzaine at acm.caltech.edu - http://geuz.org