Multivariate Hermite interpolation with interpolation conditions induced by harmonic polynomials