Parametric PDE Solving on Viscous-regularized Buckley–Leverett (ε = 10−2) (test snapshots)