Extending Functional kriging to a multivariate context