Vertical gradients using libmpdata++

[pdziekan]

slvr_dim.hpp has functions for calculating vertical gradients that also do some halo filling.
Replace that with proper halo filling (e.g. xchng_scalar) and nabla::grad calls.