#28 fix grid fix to values
the check 'grids do not fit to valueOriginal' [quantity.Discrete, line 102] only checks for absolute number of grid points and not if each grid has the right number!
The following should not work!
z = linspace(0,1,31).';
zeta = linspace(0,1,51);
[zGrid, zetaGrid] = ndgrid(z, zeta);
Z = quantity.Domain("z", z);
Zeta = quantity.Domain("zeta", zeta);
c = quantity.Discrete(zGrid, [Z, Zeta], 'name', 'a');
d = quantity.Discrete(zetaGrid', [Z, Zeta], 'name', 'b');