<div dir="ltr"><div>Thanks for that Christophe. Yes I have seen that already, but it does not tell me which lines are labelled 2, 3 etc. My current solution is to take the Ingram gmsh, where the extrusion of the airfoil  has been identified, presumably using the GUI, and replace the airfoil coordinates by my new values. I just hope that changing the number of nodes in the airfoil curve will not change the index. Short of reading the source code this is not documented?</div><div><br></div><div>Tom<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 8 Feb 2019 at 17:59, Christophe Geuzaine <<a href="mailto:cgeuzaine@uliege.be">cgeuzaine@uliege.be</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
> On 1 Dec 2018, at 12:46, Tom Peachey <<a href="mailto:tcp.free@gmail.com" target="_blank">tcp.free@gmail.com</a>> wrote:<br>
> <br>
> <br>
> I am planning to run OpenFoam batch jobs, airfoil simulations, using gmsh for the meshing. For this I need to automate the creation of .geo files. It’s straightforward, following the Ingram tutorial<br>
> <a href="https://community.dur.ac.uk/g.l.ingram/download/gmshtutorial2012.pdf" rel="noreferrer" target="_blank">https://community.dur.ac.uk/g.l.ingram/download/gmshtutorial2012.pdf</a><br>
>  until the final stage where one defines the “Physical Surfaces”<br>
> Physical Surface("back") = {1027};<br>
> presumably  a surface created by the previous “Extrude” command. As Ingram says: “Finding the number associated with a particular surface is a bit of a pain.” It seems he uses the gmsh GUI which is not available in my automated scenario.<br>
> <br>
> The Octave code by César A. Vecchio Toloy gives some more detail:<br>
> fprintf (fid, 'j5[] = Extrude {0,0,%.10g} {Surface{5};Layers{1};Recombine;};\n', cuerda/10);<br>
> fprintf (fid, 'j6[] = Extrude {0,0,%.10g} {Surface{6};Layers{1};Recombine;};\n', cuerda/10);<br>
> followed by<br>
> fprintf (fid, 'Physical Surface("inlet") = {j5[5],j6[5]};\n');<br>
> But again it is not clear how to choose the correct components of j5 and j6.<br>
<br>
Extrude returns the list of generated entities: the “top” of the extruded entity at index 0 and the extruded entity at index 1, followed by the “sides” of the extruded entity at indices 2, 3, etc. <br>
<br>
<a href="http://gmsh.info/dev/doc/texinfo/gmsh.html#Extrusions" rel="noreferrer" target="_blank">http://gmsh.info/dev/doc/texinfo/gmsh.html#Extrusions</a><br>
<br>
Christophe<br>
<br>
<br>
> <br>
> Any help greatly appreciated.<br>
> _______________________________________________<br>
> gmsh mailing list<br>
> <a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
> <a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/listinfo/gmsh</a><br>
<br>
— <br>
Prof. Christophe Geuzaine<br>
University of Liege, Electrical Engineering and Computer Science <br>
<a href="http://www.montefiore.ulg.ac.be/~geuzaine" rel="noreferrer" target="_blank">http://www.montefiore.ulg.ac.be/~geuzaine</a><br>
<br>
<br>
<br>
</blockquote></div>