//characteristic length
lc = 8e-2;

//domain = [0,1] x [0,1] x [0,1]

//domain corners
Point(101) = {0, 0, 0, lc};
Point(102) = {1, 0, 0, lc};
Point(103) = {1, 1, 0, lc};
Point(104) = {0, 1, 0, lc};
Point(105) = {0, 0, 1, lc};
Point(106) = {1, 0, 1, lc};
Point(107) = {1, 1, 1, lc};
Point(108) = {0, 1, 1, lc};

//edges
Line(111) = {101, 102};
Line(112) = {102, 103};
Line(113) = {103, 104};
Line(114) = {104, 101};
Line(115) = {105, 106};
Line(116) = {106, 107};
Line(117) = {107, 108};
Line(118) = {108, 105};
Line(119) = {101, 105};
Line(120) = {102, 106};
Line(121) = {103, 107};
Line(122) = {104, 108};
Line(123) = {102, 104};
Line(124) = {104, 105};
Line(125) = {105, 102};

//domain surfaces
Line Loop(131) = {111, 123, 114};
Line Loop(132) = {112, 113, -123};
Line Loop(133) = {-114, 124, -119};
Line Loop(134) = {-124, 122, 118};
Line Loop(135) = {115:118};
Line Loop(136) = {-120, 112, 121, -116};
Line Loop(137) = {-117, -121, 113, 122};
Line Loop(138) = {115, -120, -125};
Line Loop(139) = {-111, 119, 125};
Plane Surface(151) = {131};
Plane Surface(152) = {132};
Plane Surface(153) = {133};
Plane Surface(154) = {134};
Plane Surface(155) = {135};
Plane Surface(156) = {136};
Plane Surface(157) = {137};
Plane Surface(158) = {138};
Plane Surface(159) = {139};

//interface surface
Line Loop(140) = {123:125};
Plane Surface(160) = {140};


//domain
Surface Loop(171) = {151:159};
Volume(172) = {171};
Physical Volume(1) = {172};

//interface
Surface{160} In Volume{172};
Physical Surface(2) = {160};


