Point(1) = {0.0,0.0,0.0};
Point(2) = {1,0.0,0.0};
Point(3) = {1,1,0.0};
Point(4) = {0,1,0.0};

Point(5) = {0.0,0.0,1.0};
Point(6) = {1,0.0,1.0};
Point(7) = {1,1,1.0};
Point(8) = {0,1,1.0};

Line(1) = {4,3};
Line(2) = {3,2};
Line(3) = {2,1};
Line(4) = {1,4};

Line(5) = {8,7};
Line(6) = {7,6};
Line(7) = {6,5};
Line(8) = {5,8};

Line(9) = {4,8};
Line(10) = {5,1};
Line(11) = {3,7};
Line(12) = {6,2};

Line Loop(1) = {2,3,4,1};
Line Loop(2) = {5,6,7,8};
Line Loop(3) = {9,-8,10,4};
Line Loop(4) = {11,6,12,-2};
Line Loop(5) = {1,11,-5,-9};
Line Loop(6) = {-3,-12,7,10};

Plane Surface(1) = {1};
Plane Surface(2) = {2};
Plane Surface(3) = {3};
Plane Surface(4) = {4};
Plane Surface(5) = {5};
Plane Surface(6) = {6};

Physical Point("X Node") = {1,4,5,8};
Physical Point("X+ Node") = {2,3,6,7};

Surface Loop(1) = {1,2,3,4,5,6};
Physical Surface("Z 0") = {1};
Physical Surface("Z 1") = {2};
Physical Surface("X 0") = {3};
Physical Surface("X 1") = {4};
Physical Surface("Y 1") = {5};
Physical Surface("Y 0") = {6};

Volume(1) = {1};
Physical Volume("single crystal") = {1};


//+
MeshSize {6, 7, 3, 2, 8, 4, 1, 5} = 1.0;