function RD_durosp_xbeach(info, dirs)

d = dir(dirs.output);

RDduros = [];
RDdpp = [];
RDxbeach = [];

for i = 1:numel(d)
    if numel(d(i).name) > 3 && strcmpi(d(i).name(end-3:end),'.mat')
        load(d(i).name);
        if isfield(RD, 'duros') && isfield(RD, 'xbeach') && isfield(RD, 'dpp')
            RDduros = [RDduros RD.duros];
            RDdpp = [RDdpp RD.dpp];
            RDxbeach = [RDxbeach RD.xbeach];
        end
    end
end

RDall = [RDxbeach(:); RDduros(:); RDdpp(:)];
limits = [min(RDall) max(RDall)];
plot(RDxbeach, RDduros, 'o', 'DisplayName', 'XBeach vs. DUROS+')%, ...
hold on
plot(RDxbeach, RDdpp, '+', 'DisplayName', 'XBeach vs. D++')%, ...
    plot(limits, limits, '-', 'DisplayName', '1:1');
legend show

xlabel('XBeach retreat distance [m]')
ylabel('DUROS+ resp. D++ retreat distance [m]')
xlim(limits)
ylim(limits)

hold off
xb_write_plot(gcf, dirs.output, mfilename);