Adding Constrained Discontinuities to Gaussian Process Models of Wind Fields