Facet markers created in gmesh not correct in 2016.1.0 and 2016.2.0

Hi all!

I have a problem with the facet regions generated in gmesh.
If I create the meshfunction in 1.6.0 using the following line there is no problem.

boundaries = MeshFunction('size_t',mesh, 'my_facet_region.xml')

However, in 2016.1.0 and up the same line creates incorrect boundaries.
Any idea why this would happen?
If needed I can upload the mesh, it is pretty straightforward though.

Thanks!

closed with the note: Answer was found

asked Feb 1, 2017 by meron FEniCS User (2,340 points)
closed Feb 3, 2017 by meron

Hi, please post the mesh.

commented Feb 1, 2017 by MiroK FEniCS Expert (80,920 points)

Hi, thank you for responding.

The mesh can be downloaded here:
https://drive.google.com/open?id=0B2qbsJ0qxL6tZm9fYkRaOVZ6c2M
The facets can be found here:
https://drive.google.com/open?id=0B2qbsJ0qxL6tbVdmLURZUS1STEk

commented Feb 2, 2017 by meron FEniCS User (2,340 points)

Hi, sorry about the confusion - I meant the geo file and also your gmsh version.

commented Feb 3, 2017 by MiroK FEniCS Expert (80,920 points)

Hi MiroK, thank you for being patient with me.
This actually tipped me of on where things went wrong. I used the xml created by doflin-convert in version 1.6.0. This apparently needs to be done again with the newer version of dolfin-convert, else the facets are messed up.
Shall I close this question or could it be of interest to others?

commented Feb 3, 2017 by meron FEniCS User (2,340 points)

I think you can close it - you found the answer :)

commented Feb 3, 2017 by MiroK FEniCS Expert (80,920 points)

...