paths = ["standard_0spill", "standard_0.5spill", "standard_1spill", "standard_1.5spill"];
marks = ["o" "x" "<" "diamond"];

hold on;
for i=1:4
    bl_plot(paths(i), true, marks(i));
end
hold off
legend("0.0", "0.5", "1.0", "1.5")
figure()

hold on;
for i=1:4
    wall_units(paths(i), true, "/dev/null", marks(i));
end
p1 = plot(linspace(0, 11), linspace(0, 11));
p2 = plot( linspace(11, 1000), (1/0.4)*log(linspace(11, 1000))+5);
xlabel("$y^+$", "Interpreter", "latex")
ylabel("$U^+$", "Interpreter", "latex")
hold off
legend("0.0", "0.5", "1.0", "1.5")