diff --git a/tutorial/t6.geo b/tutorial/t6.geo index f0e5e05e693366fa797e0d4bd60f3658ba360217..0df4a2d27dfd720ed2b2a12fb6749da21e5a0028 100644 --- a/tutorial/t6.geo +++ b/tutorial/t6.geo @@ -10,26 +10,42 @@ // one Include "t1.geo"; -// Put 20 equidistant points on curve 4 -Transfinite Line{4} = 20 ; +// Delete the left line and create replace it with 3 new ones +Delete{ Surface{6}; Line{4}; } + +p1 = newp; Point(p1) = {-0.05, 0.05, 0, lc}; +p2 = newp; Point(p2) = {-0.05, 0.1, 0, lc}; + +l1 = newl; Line(l1) = {1, p1}; +l2 = newl; Line(l2) = {p1, p2}; +l3 = newl; Line(l3) = {p2, 4}; + +// Create surface +Line Loop(1) = {2, -1, l1, l2, l3, -3}; +Plane Surface(1) = {1}; // Put 20 points with a refinement toward the extremities on curve 2 Transfinite Line{2} = 20 Using Bump 0.05; +// Put 20 points total on combination of curves l1, l2 and l3 +Transfinite Line{l1} = 6; +Transfinite Line{l2} = 6; +Transfinite Line{l3} = 10; + // Put 30 points following a geometric progression on curve 1 // (reversed) and on curve 3 Transfinite Line{-1,3} = 30 Using Progression 1.2; // Define the Surface as transfinite, by specifying the four corners // of the transfinite interpolation -Transfinite Surface{6} = {1,2,3,4}; +Transfinite Surface{1} = {1,2,3,4}; // (Note that the list on the right hand side refers to points, not // curves. The way triangles are generated can be controlled by // appending "Left", "Right" or "Alternate" after the list.) // Recombine the triangles into quads -Recombine Surface{6}; +Recombine Surface{1}; // Apply an elliptic smoother to the grid Mesh.Smoothing = 100;