Improve inaccuracies and performance issues in cube_window attributes

[jcrivenaes]

The method in RegularSurface().slice_cube_window() that extracts various cube attributes between surfaces has a long story and several algorithms. From external input and my own benchmarking, it seems that algorithm 2 is not sufficient reliable and and there precision issues with algorithm 3 also. It is also evident that there are too many options, and some options are confusing.

I think it time to deprecate the current function and replace with a new, from which comparisons with similar attribute surfaces from e.g. OneSeismic shown only minor differences. Computation speed is also important.