Create residual viscosity package +rv and generalize the ResidualViscosity class - Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update. - Create separate RungekuttaRV specifically using interior RV updates - Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class

function [srfHandle, cbHandle] = logsurf(X,Y,Z, lim)
    absLogZ = log10(abs(Z));
    srfHandle = surf(X,Y,absLogZ);

    cbHandle = colorbar();
    colormap(hot(256));
    ah = gca();
    ah.CLim = lim;

    oldTickLabels = cbHandle.TickLabels;

    newTickLabels = {};

    for i = 1:length(oldTickLabels)
        newTickLabels{i} = sprintf('10^{%s}',oldTickLabels{i});
    end

    cbHandle.TickLabels = newTickLabels;
end