Add the "flux" diagnostic (uflux, vflux).

These are components of the vertically-averaged horizontal flux of ice.
Improve the "flux" diagnostic.

- Now we use "flux" to compute "velbar", not the other way around.
- Don't ignore a little bit of ice above grid.zlevels[ks].