Gaussian Process with 2 dimension input

There is an earlier discussion here: Cox Process in PyMC3