After this, the initialized global domain (with a droplet in the middle) needs to be divided between the processes (the local is set to a subset of the global where the position depends on the processor rank).
Finally for the main function, the local water elevations (after removing the ghostlines) need to be gathered on the main process: