High-resolution extensions to six Riemann solvers and three flux vector splitting schemes are developed within the framework of a reconstruction-evolution approach. Third-order spatial accuracy is achieved using two different piecewise parabolic reconstructions and a weighted essentially nonoscillatory scheme. A three-stage TVD Runge–Kutta time stepping is employed for temporal integration. The modular development of solvers provides an ease in selecting a reconstruction scheme and/or a Riemann solver/flux vector splitting scheme. The performances of these high-resolution solvers are compared for several one- and two-dimensional test cases. Based on a comprehensive assessment of the solutions obtained with all solvers, it is found that the use of the weighted essentially nonoscillatory reconstruction with the van Leer flux vector splitting scheme provides solutions for a variety of problems with acceptable accuracy.