Meshio not managing cells having several physical tags

[bragostin]

Dear Meshio users,

Using:

# Name                    Version                   Build  Channel
meshio                    5.3.4              pyhd8ed1ab_0    conda-forge

I though this issue had been adressed in #756
However, when I create a mesh with gmsh_4.10.5 in 4.1 format:

$MeshFormat
4.1 0 8
$EndMeshFormat
$PhysicalNames
8
2 1 "ev"
2 2 "cd"
2 3 "vc_in"
2 4 "bs"
2 5 "sy"
2 6 "sz"
2 7 "sc"
3 8 "ch"
$EndPhysicalNames

where basically vc_in = ev + cd in terms of cells.
Then through msh2xdmf, I get

+-----------------------------------------+
|     GMSH label     | MeshFunction value |
+--------------------+--------------------+
|     ev             |         1          |
|     cd             |         1          |
|     vc_in          |         1          |
|     bs             |         4          |
|     sy             |         5          |
|     sz             |         6          |
|     sc             |         6          |
|     ch             |         8          |
+-----------------------------------------+

where ev, cd and vc_in have the same values instead of distinct ones.

Any idea?