Require z to be a matrix

interpolate requires that z is an AbstractMatrix, ie rank 2, but contours and contour etc are happy to take any kind of data, as long as the first two dimensions match x and y.
I passed a rank 3 array (of eigenenergies over 2d momentum space), hoping to get contours over all the energies, but the error only appears in the interpolate function. Since z should be a matrix there, the higher-level functions should also only accept a matrix.

Weirdly, the MethodError isn't thrown if I pass a list of levels directly, or if those levels don't intersect with the data, because get_levels_cells returns an empty Dict, so the iteration in trace_contour doesn't happen and so the preallocated ContourLevel is returned without anything actually put in it:
https://github.com/JuliaGeometry/Contour.jl/blob/26994c47dcb9af88f2dc1961fb5291c8969afd1a/src/Contour.jl#L291
which relates to #80

Start with the higher-level contours and contour functions, then inspect interpolate and src/Contour.jl around line 291, including get_levels_cells and trace_contour. Reproduce the rank-3 input and the empty-level case, then add consistent validation so higher-level functions reject z values that are not matrices before contour tracing begins.