a = 1.0; coarse = 0.1; fine = 0.01; Point (1) = {0, 0, 0, coarse}; Point (2) = {0, 0, a, coarse}; Point (3) = {0, a, 0, coarse}; Point (4) = {0, a, a, coarse}; Point (5) = {a, 0, 0, coarse}; Point (6) = {a, 0, a, coarse}; Point (7) = {a, a, 0, coarse}; Point (8) = {a, a, a, coarse}; Point (10)= {0.5*a, 0.5*a, a, fine}; Line (1) = {1, 2}; Line (2) = {1, 3}; Line (3) = {3, 4}; Line (4) = {4, 8}; Line (5) = {4, 2}; Line (6) = {2, 6}; Line (7) = {6, 8}; Line (8) = {8, 7}; Line (9) = {7, 3}; Line (10) = {7, 5}; Line (11) = {5, 1}; Line (12) = {5, 6}; Line (13)= {10, 2}; Line (14)= {10, 4}; Line (15)= {10, 6}; Line (16)= {10, 8}; Line Loop (14) = {10, 12, 7, 8}; Plane Surface (14) = {14}; Line Loop (16) = {11, 1, 6, -12}; Plane Surface (16) = {16}; Line Loop (18) = {5, -1, 2, 3}; Plane Surface (18) = {18}; Line Loop (20) = {9, 3, 4, 8}; Plane Surface (20) = {20}; Line Loop (22) = {10, 11, 2, -9}; Plane Surface (22) = {22}; Line Loop (31)= {14,4,-16}; Line Loop (32)= {5,-13,14}; Line Loop (33)= {6,-15,13}; Line Loop (34)= {7,-16,15}; Plane Surface (31) = {31}; Plane Surface (32) = {32}; Plane Surface (33) = {33}; Plane Surface (34) = {34}; Surface Loop (26) = {22, 14, 16, 18, 31, 32, 33, 34, 20}; Volume (26) = {26}; Mesh 2; Mesh 3; OptimizeMesh "Gmsh"; RefineMesh; OptimizeMesh "Gmsh"; Save "krychle_variable.msh2";