Consider Coupling the equations defined on two different subdomains.
You may want to look How to interpolate between non-matching meshes in parallel? and Interpolation on non-matching meshes.
See Mesh generation and conversion with GMSH and PYGMSH tutorial.
I also don’t know how to do that, you might have to open a new post to ask that or can continue on How to rotate a mesh in FEniCs?
Consider How to match the nodes at the boundaries of two different domains, Multiple domains with GMSH: Any tutorial available? and Cantilever with two “domains”.