s = readtable("2505020301ODV.txt");
%%
n = max(s.Station);
cast = {'U','D'};
for i = 1:20
    tl = tiledlayout(1,2);
    title(tl, ['Station ', num2str(i)])
    for j = 1:2
        idx = s.Station == i & ismember(s.Cast_Direction, char(cast(j)));
        try
            nexttile
            scatter(s.pHin(idx),s.pres(idx),'.');
            set(gca,'YDir','reverse')
            ylim([0 max(s.pres(idx))])
            title(char(cast(j)))
        catch
            continue;
        end
    end 
    pause(1)
end