avebox: do not base detection of depth dimensions of e3t etc.
Because if e.g. z_l
is not in any of the variables dimensions, but in the coordinates, it will not be averaged over.
Another way of solving this could be the inclusion of the W-grid into the grid definitions...