Nicer implementation of level-slicing functionality.

NOTE: "level" parameter should now be a string, e.g. "1000 hPa".

- unit conversion handled automatically
- Subset is done before data actually gets read, so things happen
  much faster when loading a single level
- Data probe issue should be fixed

[1589]