dg1d: Fix convention for Bernstein polynomials (!211)

We now establish the convention that all basis functions are indexed by one.

Also had to update the `modal_bernstein_advection_sine` reftest data, because !201 skipped running CI on main and me cancelling some other CI runs.