def select_domain(self, x):
inq = self.in_domain(x)<EOL>if   isinstance(x, Mesh): return x.submesh(inq)<EOL>elif pimms.is_vector(x):  return (None if inq[<NUM_LIT:0>] else x)<EOL>x = np.asarray(x)<EOL>return x[:,inq] if x.shape[<NUM_LIT:0>] == <NUM_LIT:3> else x[inq]<EOL>
proj.select_domain(x) yields a subset of the coordinates in x that lie in the domain of
  the given map projection proj, assuming x is a 3D coordinate matrix.
proj.select_domain(x0) yields either the point x0 or None.
proj.select_domain(mesh) yields a copy of mesh that has been sub-sampled using
  mesh.submesh().