investigate flux artifacts after adaptive grid refinement
This issue is already reported in the PDELab repository.
After a refinement step, the flux shows artifacts whose appearance changes depending on the FEM order. After another time step, they disappear. It has to be determined whether this is due to the diffusive dynamics of the Richards equation or due to a visualization error.