changes for HPC3 presentation
fixed riemann_solution function to handle scalar problems better, and to also print out amdq, apdq