// box2.geo - Gmsh test file
// Unit cube with "Physical Surfaces"

// Corners
Point(1) = {0, 0, 0, 1};
Point(2) = {1, 0, 0, 1};
Point(3) = {0, 1, 0, 1};
Point(4) = {1, 1, 0, 1};
Point(5) = {0, 0, 1, 1};
Point(6) = {1, 0, 1, 1};
Point(7) = {0, 1, 1, 1};
Point(8) = {1, 1, 1, 1};

// Bottom face
Line(1) = {1, 2};
Line(2) = {2, 4};
Line(3) = {4, 3};
Line(4) = {3, 1};
Line Loop(5) = {1, 2, 3, 4};
Ruled Surface(5) = {5};

// Top face
Line(6) = {5, 6};
Line(7) = {6, 8};
Line(8) = {8, 7};
Line(9) = {7, 5};
Line Loop(10) = {6, 7, 8, 9};
Ruled Surface(10) = {10};

// Left face
Line(11) = {1, 3};
Line(12) = {3, 7};
Line(13) = {7, 5};
Line(14) = {5, 1};
Line Loop(15) = {11, 12, 13, 14};
Ruled Surface(15) = {15};

// Right face
Line(16) = {2, 4};
Line(17) = {4, 8};
Line(18) = {8, 6};
Line(19) = {6, 2};
Line Loop(20) = {16, 17, 18, 19};
Ruled Surface(20) = {20};

// Front face
Line(21) = {1, 2};
Line(22) = {2, 6};
Line(23) = {6, 5};
Line(24) = {5, 1};
Line Loop(25) = {21, 22, 23, 24};
Ruled Surface(25) = {25};

// Back face
Line(26) = {3, 4};
Line(27) = {4, 8};
Line(28) = {8, 7};
Line(29) = {7, 3};
Line Loop(30) = {26, 27, 28, 29};
Ruled Surface(30) = {30};

// Volume
Surface Loop(1) = {5, 10, 15, 20, 25, 30};
Volume(1) = {1};

// Physical Surfaces and Volumes
Physical Surface(1) = {5};
Physical Surface(2) = {10};
Physical Surface(3) = {15};
Physical Surface(4) = {20};
Physical Surface(5) = {25};
Physical Surface(6) = {30};
//Physical Volume(7) = {1};
