// an structured rectangle which creates an axi-symmetric
// solid cylinder of radius 0.5 and length (height) 2
SetFactory("OpenCASCADE");
Rectangle(1) = {0, 0, 0, 0.5, 2};

// entities for BCs and volume
Physical Curve("bottom", 1) = {1};
Physical Curve("top", 2) = {3};
Physical Surface("bulk", 3)  = {1};

// structured grid
Transfinite Line {1,3} = 1+0.5/0.01;
Transfinite Line {2,4} = 1+2/0.01;
Transfinite Surface "*";
Recombine Surface {1};           // quads instead of triangs
Mesh.ElementOrder = 2;           // second-order quads
Mesh.SecondOrderIncomplete = 0;  // force quads9 (if set to true, quad8 are created)