The work presented in this paper attempts to improve the grid (and consequently solution) for compressible-flow simulations performed with a 2D finite-volume Euler solver for unstructured grids, using mesh adaptation based on gradients of flow parameters, (pressure and pseudo-entropy, and grid geometry (cell areas, nodal distances etc.). The procedure requires solution (primitive variables) reconstruction at grid nodes, which are interpolated using linear polynomials in 2 dimensions or inverse distance based methods, and cell-averaged gradients, which are computed using the Green-Gauss method. The adaption is terminated if the global maximum displacement of any node is less than ε, where ε is a small user defined length scale. Results indicate that the method is capable of clustering the grid near an oblique shock and a contact wave, which results in sharper resolution of the discontinuities. copyright © Crown copyright (2018).All right reserved.