I would put a print statement in your point method (Riemann solver) and print all the incoming values, face center, normal, L/R states, etc. Then 'sort' the output, because the ordering will be different, and 'diff' the good and bad runs. See what changes.
>
