<div dir="auto">I agree with Nathan's advice.<div dir="auto"><br></div><div dir="auto">The source of your issue might be your geometry. If you created the surface using gmsh, you can provide the .geo file. If you imported the surface from CAD, there might be some spurious geometrical vertex that are forced to belong to the triangulation, hence giving you non-zero Z components. If you specify a Physical Surface on your .geo, it might help.</div><div dir="auto"><br></div><div dir="auto">The solution of looping on your mesh file via Python (or your script language of choice) to substitute all z by zeros is also valid.</div><div dir="auto"><br></div><div dir="auto">Take care and good luck,</div><div dir="auto"><br></div><div dir="auto">F</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 3, 2019, 4:17 PM Nathan J. Neeteson <<a href="mailto:nneeteson@rglinc.com">nneeteson@rglinc.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Ricardo,<br>
<br>
If the issue is simply that you have some vertices with non-zero z-components, and you want them to all be zero, and if the mesh file is ASCII format, could you just write a script (say, in python) to automatically loop through the $Nodes section of the file and replace the z component of every vertex position with "0.0" in a new file? I don't think any derived geometric data (such as cell areas or volumes) is in a *.mesh file, so naively I would think this should work, and the script in question would probably be no more than a few dozen lines.<br>
<br>
Thanks,<br>
<br>
Nathan Neeteson, M.Sc., E.I.T.<br>
Flow Control Research Engineer<br>
RGL Reservoir Management Inc.<br>
Corporate Head Office<br>
P 780.851.8243 | C 613.929.6283<br>
<a href="mailto:nneeteson@rglinc.com" target="_blank" rel="noreferrer">nneeteson@rglinc.com</a> | <a href="http://rglinc.com" rel="noreferrer noreferrer" target="_blank">rglinc.com</a><br>
API Q1™ and ISO 9001:2015 certified facilities.<br>
<br>
-----Original Message-----<br>
From: Ricardo Ruiz Baier [mailto:<a href="mailto:ricardo.ruiz.baier@gmail.com" target="_blank" rel="noreferrer">ricardo.ruiz.baier@gmail.com</a>]<br>
Sent: June 3, 2019 11:32 AM<br>
To: Nathan J. Neeteson <<a href="mailto:nneeteson@rglinc.com" target="_blank" rel="noreferrer">nneeteson@rglinc.com</a>>; Ricardo Ruiz Baier <<a href="mailto:ruizbaier@maths.ox.ac.uk" target="_blank" rel="noreferrer">ruizbaier@maths.ox.ac.uk</a>><br>
Cc: <a href="mailto:gmsh@onelab.info" target="_blank" rel="noreferrer">gmsh@onelab.info</a><br>
Subject: Re: [Gmsh] dilate to flatten a mesh<br>
<br>
Dear Nathan,<br>
<br>
thank you very much. This indeed resolves the syntax problem. However in my case it does not do anything since transformations can only be applied to the geometry, and in this case I want to apply the dilation to the mesh. I'll keep looking for other solutions.<br>
<br>
All the best,<br>
<br>
Ricardo<br>
<br>
> Hi Ricardo,<br>
><br>
> The syntax for anisotropic dilation is as follows:<br>
><br>
> Dilate { { expression-list }, { expression, expression, expression } }<br>
> { transform-list }<br>
><br>
> The first expression-list is the homethetic center of the dilation.<br>
><br>
> The set of three expressions are the X,Y,Z scaling factors of the dilation.<br>
><br>
> The transform-list is the entities to be acted on.<br>
><br>
> I think your syntax should be:<br>
><br>
> Dilate{ {0,0,0}, {1,1,0} }{ Point{all_points[]}; }<br>
><br>
> Assuming you want all points on the x-y plane at z=0.<br>
><br>
> Nathan Neeteson, M.Sc., E.I.T.<br>
> Flow Control Research Engineer<br>
> RGL Reservoir Management Inc.<br>
> Corporate Head Office<br>
> P 780.851.8243 | C 613.929.6283<br>
> <a href="mailto:nneeteson@rglinc.com" target="_blank" rel="noreferrer">nneeteson@rglinc.com</a> | <a href="http://rglinc.com" rel="noreferrer noreferrer" target="_blank">rglinc.com</a><br>
> API Q1(tm) and ISO 9001:2015 certified facilities.<br>
><br>
> -----Original Message-----<br>
> From: gmsh [mailto:<a href="mailto:gmsh-bounces@ace20.montefiore.ulg.ac.be" target="_blank" rel="noreferrer">gmsh-bounces@ace20.montefiore.ulg.ac.be</a>] On Behalf<br>
> Of Ricardo Ruiz Baier<br>
> Sent: June 3, 2019 5:51 AM<br>
> To: <a href="mailto:gmsh@onelab.info" target="_blank" rel="noreferrer">gmsh@onelab.info</a><br>
> Subject: [Gmsh] dilate to flatten a mesh<br>
><br>
> Dear all,<br>
><br>
> I have a .mesh triangulation which is supposed to be 2D. However it has some vertices with (spurious) crazy large z-components. How can I flatten these from e.g. a geo file?<br>
><br>
> I've tried with "Dilate" in the manner mentioned below, but it does not work (the compilation complains about the syntax and I have not found any examples).<br>
><br>
> Merge "testMesh.mesh";<br>
><br>
> all_points[] = Point '*';<br>
><br>
> Dilate { { all_points[] }, { 1, 1, 0} } {Surface(1) };<br>
><br>
> Thank you so much!<br>
><br>
> Kind regards,<br>
><br>
><br>
><br>
> --<br>
> Ricardo Ruiz Baier<br>
> Mathematical Institute<br>
> University of Oxford<br>
> <a href="http://people.maths.ox.ac.uk/ruizbaier" rel="noreferrer noreferrer" target="_blank">people.maths.ox.ac.uk/ruizbaier</a><br>
><br>
> _______________________________________________<br>
> gmsh mailing list<br>
> <a href="mailto:gmsh@onelab.info" target="_blank" rel="noreferrer">gmsh@onelab.info</a><br>
> <a href="https://can01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fonela" rel="noreferrer noreferrer" target="_blank">https://can01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fonela</a><br>
> <a href="http://b.info" rel="noreferrer noreferrer" target="_blank">b.info</a>%2Fmailman%2Flistinfo%2Fgmsh&amp;data=02%7C01%7Cnneeteson%40rgli<br>
> <a href="http://nc.com" rel="noreferrer noreferrer" target="_blank">nc.com</a>%7C47bb41a5f6d1431e0c9d08d6e84967cf%7C1f2819e2b6f842cd9d720c61c1<br>
> a98107%7C0%7C0%7C636951799311766160&amp;sdata=UWGz6UOQI3HwjrCcF032KsJo<br>
> TPbQSNtv5R7Naoe6aLY%3D&amp;reserved=0<br>
> Email disclaimer located at<br>
> <a href="https://can01.safelinks.protection.outlook.com/?url=http%3A%2F%2Frglin" rel="noreferrer noreferrer" target="_blank">https://can01.safelinks.protection.outlook.com/?url=http%3A%2F%2Frglin</a><br>
> <a href="http://c.com" rel="noreferrer noreferrer" target="_blank">c.com</a>%2Fdisclaimer&amp;data=02%7C01%7Cnneeteson%<a href="http://40rglinc.com" rel="noreferrer noreferrer" target="_blank">40rglinc.com</a>%7C47bb41a<br>
> 5f6d1431e0c9d08d6e84967cf%7C1f2819e2b6f842cd9d720c61c1a98107%7C0%7C1%7<br>
> C636951799311766160&amp;sdata=0UeeOq8ZUzxetIh20vxRSt4AxnqJLVbenaD9tcmO<br>
> sD4%3D&amp;reserved=0<br>
<br>
<br>
--<br>
Ricardo Ruiz Baier<br>
Mathematical Institute<br>
University of Oxford<br>
<a href="http://people.maths.ox.ac.uk/ruizbaier" rel="noreferrer noreferrer" target="_blank">people.maths.ox.ac.uk/ruizbaier</a><br>
Email disclaimer located at <a href="http://rglinc.com/disclaimer" rel="noreferrer noreferrer" target="_blank">http://rglinc.com/disclaimer</a><br>
_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info" target="_blank" rel="noreferrer">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer noreferrer" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a><br>
</blockquote></div>