variable rho implemented
rho*u fluxes implemented