1
// gmsh -2 my_cube_PH_gmsh.mshcad -o my_cube_PH-10.msh
2
// ../pbin/msh2geo < my_cube_PH-10.msh | ../pbin/geo -upgrade -geo - > my_cube_PH-10.geo
3
// Omega = ]a,b[x]c,d[x]f,g[
8
Point(1) = {a, c, f, hloc};
9
Point(2) = {b, c, f, hloc};
10
Point(6) = {b, c, g, hloc};
11
Point(5) = {a, c, g, hloc};
16
Line Loop(300) = {1,6,-9,-5};
18
Plane Surface(3) = {300} ;
20
Mesh.RecombinationAlgorithm = 0;
21
Recombine Surface {3} = angle; // mixt qua & tri
22
extr[] = Extrude{0,d-c,0}{ Surface{3}; Layers{n}; Recombine;};
23
Physical Surface("boundary")= {-extr[2], -extr[4], 3, -extr[3], -extr[0], -extr[5]};
24
Physical Surface("bottom") = {-extr[2]};
25
Physical Surface("top") = {-extr[4]};
26
Physical Surface("left") = {3};
27
Physical Surface("front") = {-extr[3]};
28
Physical Surface("right") = {-extr[0]};
29
Physical Surface("back") = {-extr[5]};
30
Physical Volume("interior") = {extr[1]};