%Turbular viscosity
clear
clc
clf
wsk=fopen('H:\chwila\TS\MatLab\turb_visc.dat','r');
opis1=textscan(wsk,'%s',1,'headerlines',1,'delimiter','');
linia1=textscan(wsk,'%f%f',8);
opis2=textscan(wsk,'%s',1,'delimiter','');
linia2=textscan(wsk,'%f%f',8);
opis3=textscan(wsk,'%s',1,'delimiter','');
linia3=textscan(wsk,'%f%f',7);
opis4=textscan(wsk,'%s',1,'delimiter','');
linia4=textscan(wsk,'%f%f',17);
opis5=textscan(wsk,'%s',1,'delimiter','');
linia5=textscan(wsk,'%f%f',15);
opis6=textscan(wsk,'%s',1,'delimiter','');
linia6=textscan(wsk,'%f%f',16);
hold on
dist1=linia1{1};
visc1=linia1{2};
dist2=linia2{1};
visc2=linia2{2};
dist3=linia3{1};
visc3=linia3{2};
dist4=linia4{1};
visc4=linia4{2};
dist5=linia5{1};
visc5=linia5{2};
dist6=linia6{1};
visc6=linia6{2};

plot(dist1,visc1,'.r')
plot(dist2,visc2,'.g')
plot(dist3,visc3,'.b')
plot(dist4,visc4,'.m')
plot(dist5,visc5,'.k')
plot(dist6,visc6,'.c')
hold off
xlabel('distance')
ylabel('velocity')
legend([opis1{1},opis2{1},opis3{1},opis4{1},opis5{1},opis6{1}],'location','SouthEast')
fclose(wsk);

