close all
clear all
% /chalmers/users/lada/python-DES-code/channel-5200-IDD-PANS-inlet-synt-MTF271/u-time-history-i70.dat

n1=7500;
n2=15000;
datau=load('u-time-history-i70.dat');
t=datau(n1:n2,1);
u1=datau(n1:n2,2);   %v_1 at point 1
u2=datau(n1:n2,3);   %v_1 at point 2
u3=datau(n1:n2,4);   %v_1 at point 3
u4=datau(n1:n2,5);   %v_1 at point 4
u5=datau(n1:n2,6);   %v_1 at point 5

dataw=load('w-time-history-i70.dat');
w1=dataw(n1:n2,2);   %v_3 at point 1
w2=dataw(n1:n2,3);   %v_3 at point 2
w3=dataw(n1:n2,4);   %v_3 at point 3
w4=dataw(n1:n2,5);   %v_3 at point 4
w5=dataw(n1:n2,6);   %v_3 at point 5

dt=t(2)-t(1)


%%%%%%%%%%%%%%%% plotting section %%%%%%%%%%%%%%%%%%%%%%%%%%
% plot u
figure(1)
plot(t,u1,'b-')
hold
% plot every 30th time step
plot(t(1:30:end),u1(1:30:end),'bo')
plot(t,u3,'r--')
axis([10  11  9 25])
xlabel('t')
ylabel('u')
handle=gca;
set(handle,'fontsize',[20])
print u_time.ps  -depsc
