How to manually generate a rectangular mesh with subdomains

Ah-ha, that's the cause....I was very puzzled by this FaceDescriptor...
So, domin, domout, and surfnr flags are never used in 2D mesh, right?

domin and domout are used if you have a 2D surface mesh in 3D. In your case they are not used.