Use maximum of left and right speeds to compute dissipation, rather than evaluating on some average state.
Fixed mistake in mesh size. n1,n2,n3,n4 are number of points, not number of intervals.