<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body>
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.WordSection1
{page:WordSection1;}
--></style>
<div class="WordSection1">
<p class="MsoNormal">According to your first reply, I extracted the points of the boundaries and I noted that, as you said, the points are ordered in a clockwise sens(Freefem only reads in a counterclockwise sens). So I created a Matlab script that reads the
points in .xyz format and creates a .geo file with the points and lines:</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">When I imported the points into a matrix, I reversed the order using the comand “flip”. Then I defined the points and lines according to the gmsh nomenclature, and when I created the mesh and exported in a .mesh format, FreeFem reads it
without problems.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks for your reply, it helped me to find another way to solve the problem.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I attach the .m script, the .xyz points and the generated .geo file.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Best Regards,</p>
<p class="MsoNormal">Jacob</p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De:</b> G. D. McBain <gdmcbain@protonmail.com><br>
<b>Enviado:</b> Monday, May 20, 2019 3:31:21 AM<br>
<b>Para:</b> Jacob Abdelfatah Ndioubnane<br>
<b>Cc:</b> gmsh<br>
<b>Asunto:</b> RE: Problem with Reading the mesh generated into freefem.</font>
<div> </div>
</div>
<div>
<div> <br>
</div>
<blockquote class="protonmail_quote" type="cite">
<div>
<p>I checked also that when the program generates the mesh it displays a warning error due to It can’t orient the normal of the Surface.<br>
</p>
</div>
</blockquote>
<div>Yes, on running 'gmsh -2 Dominio.geo', I see this too:<br>
</div>
<div><br>
</div>
<div> Warning : Could not orient normal of surface 1<br>
</div>
<div><br>
</div>
<div>I'm not immediately sure why though... I continue to suspect that this warning is associated with the error from FreeFEM about triangular elements having negative areas, but I haven't been able to work out what the problem with Dominio.brep is, i.e. why
Gmsh can't orient the Plane Surface.<br>
</div>
<div><br>
</div>
<div>I even get the same warning from Gmsh when trying to make a simple Plane Surface from either of the two Curve Loops separately. That does seem odd to me since the second, which looks like an island (or is it a lake) is simple enough to be inspected manually
in the Gmsh GUI and it looks fine.<br>
</div>
<div><br>
</div>
<div>Sorry, I'm not sure what to suggest from here...<br>
</div>
<div><br>
</div>
</div>
</body>
</html>