Problem with creating vector with scalar function components

Hi @dokken , my apologies for posting on the same thread. I am sorry for that.

I had already created a new thread yesterday [(https://fenicsproject.discourse.group/t/how-to-define-bcs-on-boundaries-of-submeshes-of-a-parent-mesh/5470)], but haven’t received any response from anyone yet. So I thought of continuing this thread.

Since my problem concerns with extraction of physical tags of boundaries from gmsh to FEniCS when defining subdomains and defining BCs on them, I am not sure if a built-in mesh would be an appropriate example to address the problem. In the link mentioned above, I have tried my best to reduce the problem to a minimal working example to illustrate the issue with a mesh imported from gmsh.

Again, I am sorry for continuing this thread. My sincere apologies for any inconvenience this may have caused.