 // Select periodic mesh by setting this to either 0 - standard, 1 - periodic
 periodic = 0;
 
-// Set the geometry order (1, 2, or 3)
+// Set the geometry order (1, 2, ..., 9)
 order = 3;
 
 // Set the element type (3 - triangles, 4 - quadrilaterals)


// Set a rotation periodicity constraint:
If (periodic)
   Periodic Line{1} = {2} Rotate{{0,0,1}, {0,0,0}, -Pi/3};
EndIf

// Tag surfaces and volumes with positive integers

R1 = 1.0;
R2 = 2.0;

// Angular size of the sector
Phi = Pi/3.0;

Point(1) = {0.0, 0, 0, 1.0};
Point(2) = {R1, 0, 0, 1.0};
Point(3) = {R2, 0, 0, 1.0};
Point(4) = {R1*Cos(Phi), R1*Sin(Phi), 0, 1.0};
Point(5) = {R2*Cos(Phi), R2*Sin(Phi), 0, 1.0};
Line(1) = {2, 3};
Line(2) = {4, 5};
Circle(3) = {2, 1, 4};
