Add support for artificial viscosity - Add option to use upwind operators for discretizing the scheme - Add option to use dissipation operators constructed from undivided differences when discretizing with narrow stencil operators

function names = getVarname(varargin)
    names = cell(size(varargin));

    for i = 1:numel(varargin)
        names{i} = inputname(i);
    end
end