<div dir="ltr"><div>Hi Bart,<br><br></div>Could you please check what happens if you remove the CreateTopology command ?<br><br>V.<br></div><div class="gmail_extra"><br><div class="gmail_quote">2014-09-03 5:17 GMT-05:00 Bart Vermeulen <span dir="ltr"><<a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all,<br>
<br>
I am trying to create a mesh starting with a STL model. This model consists of 6 different surfaces. I use the file included below (quite similar to one of the tutorials). When I load this into gmsh I can see the geometry with the different surfaces. At the corners of each of the surfaces there is a point (4 points for each surface). The Compound lines connect these points directly instead of following the geometry (see attached screenshot). Obviously the meshing fails.<br>
<br>
Any idea how to properly define the compound lines? (The surfaces ar OK)<br>
<br>
Thanks in advance,<br>
Bart Vermeulen<br>
<br>
<br>
<br>
Merge "mmbendo.stl";<br>
CreateTopology;<br>
<br>
// We can now define a compound line (resp. surface) for each discrete<br>
// line (resp. surface) in the model<br>
ll[] = Line "*";<br>
For j In {0 : #ll[]-1}<br>
  Compound Line(newl) = ll[j];<br>
EndFor<br>
ss[] = Surface "*";<br>
s = news;<br>
For i In {0 : #ss[]-1}<br>
  Compound Surface(s+i) = ss[i];<br>
EndFor<br>
<br>
// And we can create the volume based on the new compound entities<br>
Surface Loop(1) = {s : s + #ss[]-1};<br>
Volume(1) = {1};<br>
<br>
Physical Surface(1) = {s : s + #ss[]-1};<br>
Physical Volume(1) = 1;<br>
<br>
uniform = 1;<br>
If(uniform)<br>
  // uniform mesh size...<br>
  Mesh.CharacteristicLengthMin = 15;<br>
  Mesh.CharacteristicLengthMax = 15;<br>
EndIf<br>
If(!uniform)<br>
  // ... or apply a funny mesh size field, just because we can :-)<br>
  Field[1] = MathEval;<br>
  Field[1].F = "2*Sin((x+y)/5) + 3";<br>
  Background Field = 1;<br>
EndIf<br>
<br>
Mesh.RemeshAlgorithm = 0; // (0) no split (1) automatic (2) automatic only with metis<br>
Mesh.RemeshParametrization = 7; // (0) harmonic (1) conformal spectral (7) conformal finite element<br>
Geometry.HideCompounds = 0; // don't hide the compound entities<br>
Mesh.Algorithm = 6; // Frontal<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Bart Vermeulen<br>
Wageningen University, Department of Environmental Sciences<br>
Hydrology and Quantitative Water Management Group<br>
Postal address: Postbus 47, 6700AA Wageningen,The Netherlands<br>
Visiting address: Room 0.C027, Droevendaalsesteeg 3,6708PB Wageningen, The Netherlands<br>
Tel. <a href="tel:%2B31-317-482765" value="+31317482765" target="_blank">+31-317-482765</a> Fax <a href="tel:%2B31-317-419000" value="+31317419000" target="_blank">+31-317-419000</a><br>
email: <a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a><br>
<br>
</font></span><br>_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br>
<a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
<br></blockquote></div><br></div>