This paper presents a new parallel mesh generation method leading to subdomains of shape well-suited to Schur based domain decomposition methods such as the FETI and BDD solvers. Starting from a coarse mesh, subdomains meshes are created in parallel through hierarchical mesh refinement and morphing techniques. The proposed methodology aims at limiting the occurrence of known pathological situations (jagged interfaces, misplaced heterogeneity with respect to the interfaces, . . . ) that penalize the convergence of the solver. Furthermore, it enables to distribute and parallelize the mesh generation step in the early phases of the whole analysis. Besides its good behavior towards convergence, the mesh generation is thus distributed. The method is assessed, on several academical and industrial test cases, for both its parallel efficiency when creating the mesh and its capability to generate decomposition resulting in less FETI iterations.