* Use boundary velocity v_j in term (v_i*v_j) both for upwinding and
  central average
* do not use upwinding for kovasznay test
* update reference solutions for kovasznay and angeli test