SC, and from $r_{dS}=2.25\,\rm{mm}$ to $r_{dS}=2.5\,\rm{mm}$ in the last $10$ columns. In this way, the sound waves are more confined into the waveguides and, therefore, it is easier to correctly excite the S and A supermodes with the two sound sources and, in turn, allows for a better identification of the final intensity at each output. The fact that the supermodes are more confined when they are surrounded by rows of radius $r_{dS}=2.5\,\rm{mm}$ has been checked by comparing the band diagrams and the S and the A supermodes obtained with these extra thicker rows surrounding the linear defects plotted in [Fig. \[sc\_5defectes\]]{} with the diagrams and supermodes shown in Figs. \[sc\_fig2\](a) and (b). These new band diagrams have been determined by considering the surrounding rows as new linear defects of radius $r_{dS}$ and recalculating [Eq. [(\[sc\_coeff2\])]{}]{}: $$\begin{aligned}
\rho_{\mathbf{\bar{G}}}^{-1}=&\left(\frac{\rho_h}{\rho_c}-1\right)\frac{2\pi}{9a^2}\notag
\\
&\left[r_{dU}\frac{J_1(|\mathbf{\bar{G}}|r_{dU})}{|\mathbf{\bar{G}}|}e^{i\bar{G}_ya}+r_{dL}\frac{J_1(|\mathbf{\bar{G}}|r_{dL})}{|\mathbf{\bar{G}}|}e^{-i\bar{G}_ya}+\right.\notag
\\
&\left.2r_{dS}\frac{J_1(|\mathbf{\bar{G}}|r_{dS})}{|\mathbf{\bar{G}}|}\left(1/2+\cos(2\bar{G}_ya)\right)+\right.\notag
\\
&\left. 2r_0\frac{J_1(|\mathbf{\bar{G}}|r_0)}{|\mathbf{\bar{G}}|}\sum_{j=3}^{4}{\cos(j\bar{G}_ya)}\right], \,\rm{\ for\ } \mathbf{\bar{G}}\neq0.
\label{sc_coeff2_c}\end{aligned}$$ The expressions for $b_{\mathbf{\bar{G}}}^{-1}$ can be analogously calculated.

[Fig. \[sc\_figfase\]]{} shows the results for the sound propagation obtained by numerical integrating [Eq. [(\[sc\_linearized\])]{}]{} using the Finite Element method and considering that the sound sources have a Gaussian transverse pressure profile with a full-width at half maximum of $5\,\rm{mm}$, and are 