Issues with marking subdomains for 3D geometry

Perhaps this relates to your issue: A problem about the function near() - #2 by dokken

also Explicitly defining the subdomains and boundaries return warning Found no facets matching domain for boundary condition