function show1525(directory, variable)

[t1 v1] =  readsummary(fullfile(directory, ['summary_1525_' variable '-1_mean.txt']));
v1(find(v1 > 200)) = NaN;
v1(find(v1 < -200)) = NaN;

[t2 v2] =  readsummary(fullfile(directory, ['summary_1525_' variable '-2_mean.txt']));
v2(find(v2 > 200)) = NaN;
v2(find(v2 < -200)) = NaN;

[t3 v3] =  readsummary(fullfile(directory, ['summary_1525_' variable '-3_mean.txt']));
v3(find(v3 > 200)) = NaN;
v3(find(v3 < -200)) = NaN;

[t4 v4] =  readsummary(fullfile(directory, ['summary_1525_' variable '-4_mean.txt']));
v4(find(v4 > 200)) = NaN;
v4(find(v4 < -200)) = NaN;


[r c] = size(v1);

% for i = 1:5:c
%     figure;
%     plot(t1, v1(:, i), 'b');
%     hold on
%     plot(t2, v2(:, i), 'r');
%     plot(t3, v3(:, i), 'g');
%     plot(t4, v4(:, i), 'm');
%     datetick('x');
%     title(['Cell ' num2str(i)]);
%     ylabel('velocity [mm s^-1]')
% end

MAX_ALLOWED_FIGURES = 10;
figureCount = 0;
for i = 1:floor(r/100):r

    pGood = sum(~isnan(v1(i, :))) / c * 100;
    if (pGood > 60)
        figureCount = figureCount + 1;
        if (figureCount > MAX_ALLOWED_FIGURES) 
            break
        end
        figure;
        plot(v1(i, :), 'b');
        hold on;
        plot(v2(i, :), 'r');
        plot(v3(i, :), 'g');
        plot(v4(i, :), 'm');
        set(gca, 'YLim', [-150 150], 'XLim', [0 120]);
        title(datestr(t1(i)));
    end
end