for iYear = 1:size(MCresult,1)
%     ReturnPeriods   = [1e3, 1e4, 1e5];
%     ReturnFreqs     = [1e-3, 1e-4, 1e-5];
%     ReturnPeriods   = [4e3, 4e4, 1e5];
%     ReturnFreqs     = [2.5e-4, 2.5e-5, 1e-5];
    ReturnPeriods   = [2e3, 2e4, 1e5];
    ReturnFreqs     = [5e-4, 5e-5, 1e-5];
    
    semilogy(MCresult{iYear,2}.settings.x2zVariables{12}, MCresult{iYear,2}.Output.P_f,'b.-')
    hold on
    for iFORM = 1:numel(FORMresult(iYear,:,2))
        scatter(FORMresult{iYear,iFORM,3}.settings.x2zVariables{12},FORMresult{iYear,iFORM,3}.Output.P_f,'r')
        line([min(MCresult{iYear,2}.settings.x2zVariables{12}) max(MCresult{iYear,2}.settings.x2zVariables{12})],...
            [ReturnFreqs(iFORM) ReturnFreqs(iFORM)], 'Color', 'k', 'LineStyle', '--')
        RetreatDistances(iFORM) = FORMresult{iYear,iFORM,3}.settings.x2zVariables{12};
        Pfs(iFORM)              = FORMresult{iYear,iFORM,3}.Output.P_f;
    end
    
    n = numel(RetreatDistances);
    if size(FORMresult,3) > 3
        for iIteration = 4:(size(FORMresult,3))
%             PolyFit         = polyfitn(log10(Pfs),RetreatDistances,1);
%             fitRDs          = polyvaln(PolyFit,log10([1e-2 5e-2 1e-3 5e-3 1e-4 5e-4 1e-5]));
%             scatter(fitRDs, [1e-2 5e-2 1e-3 5e-3 1e-4 5e-4 1e-5],'g*')
            for iFORM = 1:size(FORMresult,2)
                if size(FORMresult,3) > 3 && ~isempty(FORMresult{iYear,iFORM,iIteration})
                    n = n+1;
                    scatter(FORMresult{iYear,iFORM,iIteration}.settings.x2zVariables{12},FORMresult{iYear,iFORM,iIteration}.Output.P_f,'m*')
                    RetreatDistances(n)  = FORMresult{iYear,iFORM,iIteration}.settings.x2zVariables{12};
                    Pfs(n) = FORMresult{iYear,iFORM,iIteration}.Output.P_f;
                end
            end
        end
    end   
    hold off
    title(num2str(MCresult{iYear,1}))
    keyboard
end