examples/scuff-em-benchmarks/shapes/sphere.geo.template

   1 //
   2 // gmsh geometry specification for a sphere of radius R = RADIUS nm
   3 //
   4
   5 //************************************************************
   6 //* input parameters
   7 //************************************************************
   8 R = RADIUS/1000.;    // radius
   9
  10 //************************************************************
  11 //* meshing finenesses ***************************************
  12 //************************************************************
  13 l3 = ELEMSIZ/1000.;  // fineness at north pole
  14 l2 = ELEMSIZ/1000.;  // fineness at equator
  15 l1 = ELEMSIZ/1000.;  // fineness at south pole
  16
  17 //************************************************************
  18 //* upper sphere *********************************************
  19 //************************************************************
  20 Point(1) = {  0 ,    0,  0.0,  l2};
  21 Point(2) = {  R,    0,  0.0,  l2};
  22 Point(3) = {  0 ,   R,  0.0,  l2};
  23 Circle(1) = {2,1,3};
  24 Point(4) = { -R,    0,  0.0,  l2};
  25 Point(5) = {   0,  -R,  0.0,  l2};
  26 Circle(2) = {3,1,4};
  27 Circle(3) = {4,1,5};
  28 Circle(4) = {5,1,2};
  29 Point(6) = {   0,    0,  0.0+R, l3};
  30 Point(7) = {   0,    0,  0.0-R, l1};
  31 Circle(5) = {3,1,6};
  32 Circle(6) = {6,1,5};
  33 Circle(7) = {5,1,7};
  34 Circle(8) = {7,1,3};
  35 Circle(9) = {2,1,7};
  36 Circle(10) = {7,1,4};
  37 Circle(11) = {4,1,6};
  38 Circle(12) = {6,1,2};
  39 Line Loop(13) = {2,8,-10};
  40 Ruled Surface(14) = {13};
  41 Line Loop(15) = {10,3,7};
  42 Ruled Surface(16) = {15};
  43 Line Loop(17) = {-8,-9,1};
  44 Ruled Surface(18) = {17};
  45 Line Loop(19) = {-11,-2,5};
  46 Ruled Surface(20) = {19};
  47 Line Loop(21) = {-5,-12,-1};
  48 Ruled Surface(22) = {21};
  49 Line Loop(23) = {-3,11,6};
  50 Ruled Surface(24) = {23};
  51 Line Loop(25) = {-7,4,9};
  52 Ruled Surface(26) = {25};
  53 Line Loop(27) = {-4,12,-6};
  54 Ruled Surface(28) = {27};
  55 Physical Surface(1) = {28,26,16,14,20,24,22,18};
  56
  57 //************************************************************
  58 //* reference point to get outward-pointing surface normals right
  59 //************************************************************
  60 Physical Point(1) = {1};