L  = 10              ;
H  = 5               ;
W  = 5               ;
d  = 0.5             ;

// ==================================
//         G E O M E T R Y 
// ==================================
 
Point(1)    = {0        , 0     , 0 , d};
Point(2)    = {L        , 0     , 0 , d};
Point(3)    = {L        , H     , 0 , d};
Point(4)    = {0        , H     , 0 , d};
Point(5)    = {0        , 0     , W , d};
Point(6)    = {L        , 0     , W , d};
Point(7)    = {L        , H     , W , d};
Point(8)    = {0        , H     , W , d};

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

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

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


Line Loop(1)  = {-4,-3,-2,-1};
Line Loop(2)  = {5,6,7,8};
Line Loop(3)  = {-11,3,12,-7};
Line Loop(4)  = {-12,4,9,-8};
Line Loop(5)  = {11,-6,-10,2};
Line Loop(6)  = {-5,-9,1,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};

Surface Loop(1) = {1,3,2,4,6,5};
Volume(1) = {1};

// ================================
//           G R O U P S  
// ================================

// DIRICHLET BC 
Physical Point("Fix"  ,        1)   = {1,4,5,8};
Physical Curve("Fix"  ,        1)   = {4,8,9,12};
Physical Surface("Fix",        1)   = {4};

// NEUMANN BC 
Physical Point("ForceP",        2)   = {2,6};
//Physical Curve("ForceC",        3)   = {10};
//Physical Surface("ForceS",      4)   = {5};

// MATERIAL GROUP
Physical Volume("Solid",        5)   = {1}; 

// ================================
//         M E S H I N G  
// ================================

Mesh 3;

Save "mesh_tetra.msh";