Non-flux condition doesn't work

[RemDelaporteMathurin]

The solution yields nan when non flux condition is put on one of the sides

[RemDelaporteMathurin]

Found the problem. It comes from the fact that we split u into 4 sub-solutions but don't solve the fourth one.