Derivative-constrained gaussian processes?

Is there a way that I can apply constraints on GP at certain x points?
I found this notebook (Jupyter Notebook Viewer) in which a derivative constraint at a point is given with potential. But the notebook is from 2012, and I have no idea how to reproduce it with the recent pymc3 version.