NEW FEATURE: Add a formatoption to calculate data-driven formatoptions based on the entire data

Summary

It should be possible to calculate bounds, ylim, etc. based on the entire data

Reason

Particularly if one would want to update a dimension (e.g. choose the next time-step), it would be helpful to see the evolution without changing x-limits, y-limits, etc.

Detailed explanation

The following formatoptions need to be considered:

xlim, ylim
bounds, vbounds
xticks, yticks, cticks, vcticks

Don't know, what the best name would be, maybe use_full_data? @MuellerSeb, do you have a good idea? It should also support the specification of certain dimensions, e.g.

sp.update(use_full_data=True)  # to consider all dimensions
sp.update(use_full_data='time') # to only consider one dimension
sp.update(use_full_data='t')       # to only consider the time dimension, whatever it is named
sp.update(use_full_data=['time', 'lev'])  # to consider multiple dimensions

The default should be False