Note that \(h_{\mathbf {a}}\)  is the element length in the direction of the advection flow; see [1]} and [2]}. As a reference solution we are considering the one obtained by classical tensor-product B-splines of the same degrees on a very fine mesh of \(1500 \times 1500\)  elements, as shown in figadvsqrbspline.
<TABLE>