subdomain_data should be a MeshFunction and subdomain_id=3
subdomain_data
MeshFunction
subdomain_id=3