The problem is when I do uReduce = u in line 3, the values are extrapolated automatically, which is NOT zero on ThG-Th, so that line 4 generates wrong values in uG using mpiSUM.
Is there a better method to centralise/combine results form different MPI processors?
It is explained thoroughly how to do mesh adaptation, gather a solution, so on and so forth (basically all your questions from the last couple of days…).