_h = plt.contourf(xp, xp, fs2)