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 save_frame = setup_fig_mov(figure_handle,dirname)

    fresh_dir(dirname)

    figure_handle.Units = 'centimeters';
    figure_handle.PaperUnits = 'centimeters';
    figure_handle.PaperPosition(3:4) = figure_handle.Position(3:4);


    frame_nr = 0;
    function save_frame_F()
        saveas(figure_handle,sprintf('%s/%s%04d.png',dirname,'fig',frame_nr));
        frame_nr = frame_nr + 1;
    end
    save_frame = @save_frame_F;
end