The grid is composed by 120 elements along each direction and the time step is such that the maximum Courant number \(C \approx 118\)  and the maximum value of advective Courant number \(C_u\)  is around \(0.03\) . Figures REF ,REF  and REF  show the results at \(t \in \left\lbrace 10, 15, 20\right\rbrace \)  both in terms of contours and plots along the same specific cuts along \(x-\) axis chosen in [1]}. All the results are in good agreement with the reference ones and, also in this case, we are able to recover the Kelvin-Helmholtz instability.
<FIGURE><FIGURE><FIGURE>