We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography, and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results ...