function sFlow(ini,filename,clr);
Y = statEvalFlow(ini,100,0);
S = size(Y);

l = Y(:,1); Kc = Y(:,2);
if length(ini) >5
    y = Y(:,6); Db = Y(:,7); Df = Y(:,8);
else
    y = Y(:,4); Db = Y(:,5); Df = Y(:,6);
end

if strcmp(filename,'nf') == 0
	dlmwrite(filename,Y,' ');
end

% subplot(2,2,1)
hold on;
plot(Kc,Db,clr);
title('K_{c}-D_{b} phase space');
xlabel('K_{c}');
ylabel('D_{b}');
grid on;
