Calculate the buoyancy over the entire surface slice including ghost
cells
Communicate the boundary value for scalars in case of fixed fluxes for a
proper boundary buoyancy calculation
Fixed the ustar bug and with that the parallel restarts