After this, the initialized global domain (with a droplet in the middle) needs to be divided between the processes (the local u is set to a subset of the global u 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: