I'm using CFD (Computational fluid dynamics) to solve numerically the shallow water equations.
I'm writing a code based on roe solver, I'm supossed to get two vectors of solutions: one for the water depth ($h$) and the other one for the discharge ($hu$)...But when I compare the numerical solution with the analytical one, it comes that the water depth doesn't match, while the discharge it does match. Is there any way that I could have an error with the implementation if I'm able to get one of the two solutions?