Positive and S-convex C1-interpolation of gridded three dimensional data