xSliceCoords = [];
ySliceCoords = [];
zSliceCoords = [];
and isosurface values
IsosurfValues = [0.1 0.5];
IsosurfAlphas = [0.1 1.0];
IsosurfColors = ['b' 'r'];
This shows the isosurface at value 0.5 as a solid red surface, and the isosurface at value 0.1 as a wispy blue surface.
Note that the surfaces remain close together at the leading edge (the shock) and grow farther apart at the trailing edge (the rarefaction wave).
Data at t = 0:
Solution at t = 0.4:
Solution at t = 0.8: