MagparWiki: attachment:cylinder_fixed.geo of BugTracker/GmshPhysicalVolumeNumbering

   1 lc = 0.25;
   2 Point(1) = {0.0, 0.0, -1.0, lc};
   3 Point(2) = {0.2, 0.0, -1.0, lc};
   4 Point(3) = {0.0, 0.2, -1.0, lc};
   5 Point(4) = {-0.2, 0.0, -1.0, lc};
   6 Point(5) = {0.0, -0.2, -1.0, lc};
   7 
   8 Point(6) = {0.0, 0.0, 1.0, lc};
   9 Point(7) = {0.2, 0.0, 1.0, lc};
  10 Point(8) = {0.0, 0.2, 1.0, lc};
  11 Point(9) = {-0.2, 0.0, 1.0, lc};
  12 Point(10) = {0.0, -0.2, 1.0, lc};
  13 
  14 Circle(13) = {2,1,3};
  15 Circle(14) = {3,1,4};
  16 Circle(15) = {4,1,5};
  17 Circle(16) = {5,1,2};
  18 
  19 Circle(17) = {8,6,7};
  20 Circle(18) = {9,6,8};
  21 Circle(19) = {10,6,9};
  22 Circle(20) = {7,6,10};
  23 
  24 Line(30) = {2, 7};
  25 Line(31) = {3, 8};
  26 Line(32) = {4, 9};
  27 Line(33) = {5, 10};
  28 
  29 Line(34) = {7, 2};
  30 Line(35) = {8, 3};
  31 Line(36) = {9, 4};
  32 Line(37) = {10, 5};
  33 
  34 Line Loop(38) = {34,13,-35,17};
  35 Ruled Surface(39) = {38};
  36 Line Loop(40) = {-36,18,35,14};
  37 Ruled Surface(41) = {40};
  38 Line Loop(42) = {36,15,-37,19};
  39 Ruled Surface(43) = {42};
  40 Line Loop(44) = {37,16,-34,20};
  41 Ruled Surface(45) = {44};
  42 Line Loop(46) = {18,17,20,19};
  43 Plane Surface(47) = {46};
  44 Line Loop(48) = {15,16,13,14};
  45 Plane Surface(49) = {48};
  46 Surface Loop(50) = {47,41,43,49,45,39};
  47 Volume(51) = {50};
  48 Physical Volume (1) = {51};