Fix compressible solver with NUMMODES=1
This MR includes a few bug fixes required to run simulations with the compressible solver using just one mode per element.
I will not mark this as a patch because I think some of the changes are not compatible with the release branch.