<div dir="ltr">Dear Duners,<div><br></div><div style>  I use GmshReader to load a mesh with "Physical Surfaces" into an ALUGrid, then use the GridFactory's methods wasInserted(), and insertionIndex() to retrieve them. This I do into the functors I'll use for my boundary conditions.<br>

</div><div style><br></div><div style>  Apparently, so far, so good. I guess this should be pretty standard. But in the tests for the functors, I plot the support and export it to vtk and see that it's not what it should be. In this test:</div>

<div style><br></div><div style>   <a href="https://dl.dropboxusercontent.com/u/15831115/wrongsupport.png">https://dl.dropboxusercontent.com/u/15831115/wrongsupport.png</a><br></div><div style><br></div><div style>the support of the functor is set to be two of the lateral sides of the prism but it's obviously wrong (the front side is *not* among the physical surfaces in the support).</div>

<div style><br></div><div style>I've uploaded a minimal dune project here:</div><div style><br></div><div style>   <a href="https://dl.dropboxusercontent.com/u/15831115/gmshreadertest.tar.bz2">https://dl.dropboxusercontent.com/u/15831115/gmshreadertest.tar.bz2</a><br>

</div><div style><br></div><div style>Am I doing anything wrong? Is this a problem in GmshReader? I'd be most grateful for any comments or ideas.</div><div style><br></div><div style>Best,</div><div><div><div dir="ltr">

--<br>Miguel de  Benito.</div></div>
</div></div>