I understand the derivation for the boundary conditions for $\mathbf B$ and $\mathbf E$ as it was explained to me in Griffiths, but Griffiths states the following:
$$H_{\text{above}}^{\bot} - H_{\text{below}}^{\bot} = -(M_{\text{above}}^{\bot}-M_{\text{below}}^{\bot})$$
and
$$\mathbf D_{\text{above}}^{\parallel} - \mathbf D_{\text{below}}^{\parallel} = \mathbf P_{\text{above}}^{\parallel} - \mathbf P_{\text{below}}^{\parallel}$$
This feels like the right idea, due to the relations that $\nabla \times \mathbf D = \nabla \times \mathbf P$ and $\nabla \cdot \mathbf H = - (\nabla \cdot \mathbf M)$, and from what I learned of the implications by $\nabla \cdot \mathbf E = \rho / \epsilon_0 $ implies a discontinuity in the parallel component of the electric field, which is why it makes sense to me that the discontinuity in the parallel component of the $\mathbf D$ field is the way it is. However, I can't quite piece together a coherent argument justifying this boundary conditions. How is this proven?
Also, the $\mathbf H$-field perpendicular component boundary conditions are expressed in magnitude in Griffiths, yet the $\mathbf D$-field parallel component boundary conditions were expressed as vectors, as I've shown. Why is that? If I made the $\mathbf H$-field a vector in the boundary condition I wrote for it, wouldn't that still be true?