SUMMARYThe numerical simulation of physical phenomena represented by non-linear hyperbolic systems of conservation laws presents speciÿc di culties mainly due to the presence of discontinuities in the solution. State of the art methods for the solution of such equations involve high resolution shock capturing schemes, which are able to produce sharp proÿles at the discontinuities and high accuracy in smooth regions, together with some kind of grid adaption, which reduces the computational cost by using ÿner grids near the discontinuities and coarser grids in smooth regions. The combination of both techniques presents intrinsic numerical and programming di culties. In this work we present a method obtained by the combination of a high-order shock capturing scheme, built from Shu-Osher's conservative formulation (J. Comput. Phys. 1988; 77:439-471; 83:32-78), a ÿfth-order weighted essentially non-oscillatory (WENO) interpolatory technique (J. Comput. Phys. 1996; 126:202-228) and Donat-Marquina's ux-splitting method (J. Comput. Phys. 1996; 125:42-58), with the adaptive mesh reÿnement (AMR) technique of Berger and collaborators (Adaptive mesh reÿnement for hyperbolic partial di erential equations.