data=load('BenchmarkOrthogonal_Cum25_2D-DIC_TS_MP(1x32x32_50ov)_Vy_y60-y-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,1)
plot(data(:,1),data(:,k+1))

title('Across-strike extension (y = 60 mm)')
xlabel('position x (mm)')
ylabel('cumulative y-displacement (mm)')
xlim([-50 850])
ylim([-5 15])

hold on

end


data=load('BenchmarkOrthogonal_Cum25_2D-DIC_TS_MP(1x32x32_50ov)_Vy_x360-y-t.prm');

[m n]=size(data)

for k=1:2:n-1
    
subplot(2,3,4)
plot(data(:,1),data(:,k+1))

title('Across-strike extension (x = 360 mm)')
xlabel('position y (mm)')
ylabel('cumulative y-displacement (mm)')
xlim([-200 200])
ylim([-20 20])

hold on

end