
Merge "new_attempt.vtk"; 

Coherence Mesh;

CreateTopology;

CreateGeometry;

// Define the volume (assuming there is no hole)
s() = Surface{:};
Surface Loop(1) = s();
Volume(1) = 1;

// Since we did not create any new surface, we can easily define physical groups
Physical Surface("WHOLE", 1) = {1};
Physical Volume("Wall", 2) = {1};

// Set meshing constraints
uniform = 1;
If(uniform)
  Mesh.CharacteristicLengthMax = 400;
  Mesh.CharacteristicLengthMin = 400;
Else
  Mesh.CharacteristicLengthMax = 400;
  Mesh.CharacteristicLengthMin = 300;
  Mesh.CharacteristicLengthFromCurvature = 1;
  Mesh.MinimumElementsPerTwoPi = 20;
EndIf
Mesh.Algorithm3D = 10;
Mesh.MshFileVersion = 2.2;
