  We adapt Glimm's approximation method to the framework of convex integration
to show density of wild data for the (complete) Euler system of gas dynamics.
The desired infinite family of entropy admissible solutions emanating from the
same initial data is obtained via convex integration of suitable Riemann
problems pasted with local smooth solutions. In addition, the wild data belong
to BV class.
