 SetOrder order;
 Mesh.MshFileVersion = 2.2;
 
-Save Sprintf("periodic-annulus-sector-o%01g.msh",order);
+Save Sprintf("periodic-annulus-sector-t%01g-o%01g.msh",type,order);

SetFactory("OpenCASCADE");

// Set the geometry order (1, 2, or 3)
order = 1;

R1 = 1.0;
R2 = 2.0;


Curve Loop(5) = {1, 4, -2, -3};
Plane Surface(1) = {5};

Transfinite Curve{1} = nrad+1;
Transfinite Curve{2} = nrad+1;
Transfinite Curve{3} = nazm1+1;
Transfinite Curve{4} = nazm2+1;

If (nazm1 == nazm2)
   Transfinite Surface{1};
EndIf
If (type == 4)
   Recombine Surface {1};
EndIf


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