Is it possible to train a gaussian process regression (GPR) using input and output multidimensional arrays?

I want to use a GPR to reproduce the perfomance of a mechanistic model. I have already conduct simulations and gather input and output data. I have 100 time series each with 10k data and 20 characteristics, resulting in a multidimensional array of 10k x 20 x 100.

Thanks for your help, Rafael