*** empty log message ***

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;