<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<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;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="NL-BE" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Thanks again, Jeremy.<o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-GB">You seem to have nailed the origin of the problem.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">It still feels as a bug to me that reloading a script one just created, gives a different result.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Or am I missing something on how GMSH works?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Stranger still: <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Action 1: <o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:35.4pt"><span lang="EN-GB">Load the file below (without Coherence in the script).
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:35.4pt"><span lang="EN-GB">Click: 3D then Coherence then again 3D. This gives an error.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Action 2: <o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:35.4pt"><span lang="EN-GB">Load the file below (without Coherence in the script).
<o:p></o:p></span></p>
<p class="MsoNormal" style="text-indent:35.4pt"><span lang="EN-GB">Click: 3D then Coherence then 1D then again 3D. This gives the wanted results.<o:p></o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-bottom:solid windowtext 1.0pt;padding:0cm 0cm 1.0pt 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB">This seems all unintuitive to me. Am I missing something of the GMSH logic?<o:p></o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB">e.g. doesn’t clicking on the coherence command, execute that coherence command?<o:p></o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB">Bart<o:p></o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal" style="border:none;padding:0cm"><span lang="EN-GB"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">KanR = DefineNumber[ 0.1, Name "Parameters/KanR" ];<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">WanD = DefineNumber[ 0.2, Name "Parameters/WanD" ];<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">KanMeshsize = DefineNumber[ 0.01, Name "Parameters/KanMeshsize" ];<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">SetFactory("OpenCASCADE");<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Rectangle(1) = {-2.5, -2.5, 0, 5, 5, 0};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Circle(5) = {0, 0, 0, KanR, 0, 2*Pi};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Curve Loop(2) = {5};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Plane Surface(2) = {2};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Delete; }<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Extrude {0, 0, WanD } {<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"> Surface{1}; <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">}<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Point(15) = {0, KanR, 0};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Point(16) = {0, KanR-0.05, 0};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Point(17) = {0, KanR-0.05, WanD};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Point(18) = {0, KanR, WanD};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Line(20) = {15, 16};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Line(21) = {16, 17};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Line(22) = {17, 18};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Line(23) = {18, 15};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Curve Loop(24) = {20, 21, 22, 23};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Plane Surface(8) = {24};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Extrude {{0, 0, 1}, {0, 0, 0}, Pi} {<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"> Surface{8}; Layers {10}; <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">}<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[1] = Distance;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[1].EdgesList = {5};<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[1].NNodesByEdge = 50;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2] = Threshold;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2].DistMax = 2.5;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2].DistMin = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2].IField = 1;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2].LcMax = 0.25;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Field[2].LcMin = KanMeshsize ;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Background Field = 2;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Mesh.CharacteristicLengthExtendFromBoundary = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Mesh.CharacteristicLengthFromPoints = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Mesh.CharacteristicLengthFromCurvature = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<blockquote style="border:none;border-left:solid #729FCF 1.5pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
</blockquote>
</div>
</body>
</html>