If you want an example where $|mD|$ has no fixed codim-1 components, you can do the following. Let $X \dashrightarrow X^+$ be an Atiyah flop of a curve $C$, and let $H \subset X$ be a general member of a very ample linear system with $H \cdot C \geq 2$. The strict transform $\tilde{H} \subset X^+$ is still big. Let $Y$ be the resolution of the flop with exceptional divisor $E$, and $\bar{H}$ the strict transform of $H$ on $Y$.
On $Y$, we have $K_Y = g^* K_{X^+} +E$, and the strict transform is $\bar{H} = f^*H$ (since $H$ is general, it doesn't contain $C$). It's easy to check $f^*H+aE = g^*\tilde{H}$, where $a = H \cdot C$. Then $K_Y + \bar{H} = (g^*K_{X^+} + E) + f^\ast H = (g^* K_{X^+} + E) + (g^*\tilde{H} - aE)$ i.e. $K_Y + \bar{H} = g^*(K_{X^+} + \tilde{H}) + (1-a)E$. Since $a \geq 2$, the pair $(X^+,\tilde{H})$ is not numerically klt (or even numerically lc, if you take $a>2$). The problem is that $|\tilde{H}|$ is too singular along the flopped curve $C^+$; $\tilde{H}$ is already general in its linear series, so we can't get the multiplicity along $C^+$ to be any smaller, and passing to a multiple doesn't help. (Hopefully I got all my signs right!)