In this manuscript, singularly perturbed energies with 2, 3, or 4 preferred gradients subject to incompatible Dirichlet boundary conditions are studied. This extends results on models for martensitic microstructures in shape memory alloys (N=2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=2$$\\end{document}), a continuum approximation for the J1-J3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$J_1-J_3$$\\end{document}-model for discrete spin systems (N=4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=4$$\\end{document}), and models for crystalline surfaces with N different facets (general N). On a unit square, scaling laws are proved with respect to two parameters, one measuring the transition cost between different preferred gradients and the other measuring the incompatibility of the set of preferred gradients and the boundary conditions. By a change of coordinates, the latter can also be understood as an incompatibility of a variable domain with a fixed set of preferred gradients. Moreover, it is shown how simple building blocks and covering arguments lead to upper bounds on the energy and solutions to the differential inclusion problem on general Lipschitz domains.