Project/interpolate constants to different mesh subdomains

See How to define different materials / importet 3D geometry (gmsh)