<div dir="ltr"><div dir="ltr">Hello all,<div><br></div><div>I am new to gmsh and I am struggling to find what points make up a line from an imported .iges file. The goal is to import an .iges file then use a script to find which points are on the top layer and then find which lines make that surface so I can refine the mesh along that surfaces boundary. So far I can extract the point's x-y-z coordinates but I cannot find which points make up a line. I have tried this:</div><div><div>***********</div><div>//+<br></div><div>SetFactory("OpenCASCADE");</div><div>ShapeFromFile("box.iges")</div><div>//+</div><div>coords[] = Point{7};</div><div><b>point_nums[] = Line{1} // this does not work</b></div><div>c = #a[];</div><div><br></div><div>Printf("x coord = %g", coords[0]);</div><div>Printf("y coord = %g", coords[1]);</div><div>Printf("z coord = %g", coords[2]);</div></div><div>**************</div><div>The point_nums comment does not work.</div><div>I was wondering if it was possible to output the label of the two points that define a line.  I plan to find out 1) which points define the 1-D boundary I care about 2) which lines are made with those points so I can refine along the line. Or maybe there is a better way to do this all together? I essentially need the mesh to be fine along the defined line and coarse everywhere else. </div><div><br></div><div>The reason I am not doing this manually is because I have 1000+ .iges files to mesh due to the scope of my project. </div><div><br></div><div>any help would be greatly appreciated!</div><div>Thanks,</div><div>Terrence </div><div><br></div></div></div>