//gmsh-4.4.1-Windows64
SetFactory("Built-in");

//Geometry parameters
R = DefineNumber[ 5, Name "Parameters/R" ];

//Mesh parameter
h = DefineNumber[ 0.25, Name "Parameters/h" ];

//Points of the negative tape \cap sigma of the cylinder
Point(1) = {-0.5, 0., 0., h};
Point(2) = {-0.5, 0., R, h};
Point(3) = {-0.5, R, 0., h};
Point(4) = {-0.5, 0., -R, h};
Point(5) = {-0.5, -R, 0., h};
//Points of the positive tape \cap sigma of the cylinder
Point(6) = {0.5, 0., 0., h};
Point(7) = {0.5, 0., R, h};
Point(8) = {0.5, -R, 0., h};
Point(9) = {0.5, 0., -R, h};
Point(10) = {0.5, R, 0., h};

//Lines of the negative tape \cap sigma of the cylinder
Circle(1) = {2,1,3};
Circle(2) = {3,1,4};
Circle(3) = {4,1,5};
Circle(4) = {5,1,2};
//Lines of the positive tape \cap sigma of the cylinder
Circle(5) = {7,6,8};
Circle(6) = {8,6,9};
Circle(7) = {9,6,10};
Circle(8) = {10,6,7};
//Lines for the wall of the cylinder
Line(9) = {2,7};
Line(10) = {3,10};
Line(11) = {4,9};
Line(12) = {5,8};

//Curve-loops of the negative and positive tape \cap sigma of the cylinder
Curve Loop(1) = {1,2,3,4};
Curve Loop(2) = {5,6,7,8};
//Curve-loops for the wall of the cylinder
Curve Loop(3) = {-1,9,-8,-10};
Curve Loop(4) = {-2,10,-7,-11};
Curve Loop(5) = {-3,11,-6,-12};
Curve Loop(6) = {-4,12,-5,-9};

//Surfaces of the tapes, negative and positive
Plane Surface(1) = {1};
Plane Surface(2) = {2};
//Surfaces for the wall of the cylinder
Surface(3)={3};
Surface(4)={4};
Surface(5)={5};
Surface(6)={6};

//Surface Loop for the Exterior Volume
Surface Loop(1)={1,3,4,5,6,2};

//Exterior Volume
Volume(1) = {1};

//Physical groups
Physical Surface("Sigm",101) = {3,4,5,6};
Physical Surface("TapaNeg",102) = {1};
Physical Surface("TapaPos",103) = {2};
Physical Volume("Ye",201) = {1};