Polish resample script.

- improve documentation
- allow different left and right shift
- preserve original datasets in geometry
- prefix interpolated arrays with "interpol_"
- compute differences with prefix "diff_"