Thanks for responding quickly! Looking at that example more closely, it seems they first sample draws from the GP hyperprior at the X coordinates in the dataset. That works for me, though the result is just linearly-interpolated retrodictions:
The conditional method is where things go haywire.
(As a side note, the second application in that example samples from conditional, but then just plots the retrodictions. Maybe that should be more consistent?)