Dear all, I am computing the one dimensional Shu-Osher problem using DG of the CompressibleFlowSolver. I found one tricky problem. The initial condition contains a shock. When I set the initial conditions like the following <FUNCTION NAME="InitialConditions"> <E VAR="rho" VALUE="3.857143*(x<=-4.0) + (1.0+0.2*sin(5*x))*(x>-4.0)"/> <E VAR="rhou" VALUE="3.857143*2.629369*(x<=-4.0) + 0.0*(x>-4.0)" /> <E VAR="rhov" VALUE="0.0" /> <E VAR="E" VALUE="(10.33333/(Gamma-1) + 0.5*3.857143*2.629369*2.629369)*(x<=-4.0) + (1.0/(Gamma-1))*(x>-4.0)"/> </FUNCTION> spurious oscillations would occur in the initial condition around x=-4.0 where the shock is located. Although I can arrange the mesh to make the shock lie just at the interface of two elements, only one value can be specified at x=-4.0, causing the spurious oscillations. However, for the DG method, two different values should be accepted at one element interface. Can anybody help me with this ? P.S.: Where can I find the version history of Nektar++, illustrating what is new in the newer version software compared with the older one? Best, Jian Yu