Use a softened version of the interpolant as a likelihood.

def lnlike(self, X):
        """
        Use a softened version of the interpolant as a likelihood.
        """
        return -3.5*np.log(self._interpolant(X[0], X[1], grid=False))