Loop for plotting severa variables

I am trying to plot Depth against Salinity, Temperature and Oxygen collected by a CTD profiler at several stations. Lenght of data is variable depending of water depth at eact station.

My current code is as follow:

subplot(2,3,1)
hold on;
plot(T1,D1, 'b -');
plot(T2,D2, 'r --');
title('Perfil Temperatura (°C)')
yticks(0:5:100);
set(gca,'xaxisLocation','top');
set(gca,'YDir','reverse');
ylabel('m');
legend('Est.1','Est.2','Location','south' )
 
subplot(2,3,2)
hold on
plot(S1,D1, 'b -');
plot(S2,D2, 'r --');
title ('Perfil Salinidad (PSU)')
yticks(0:5:100);
set(gca,'xaxisLocation','top');
set(gca,'YDir','reverse');
ylabel('m');
legend('Est.1','Est.2', 'Location','south' )
 
subplot(2,3,3)
hold on
plot(DO1,D1, 'b -');
plot(DO2,D2, 'r --');
title ('Perfil Oxígeno (mg/L)')
yticks(0:5:100);
set(gca,'xaxisLocation','top');
set(gca,'YDir','reverse');
ylabel('m');
legend('Est.1','Est.2', 'Location','south' )

How can I use a loop to come up with a shorter code??

Answer 1

[Walter Roberson]

X1 = {T1, S1, D01};
X2 = {T2, S2, D02};
Y1 = D1;
Y2 = D2;
titles = {'Perfil Temperatura (°C)', 'Perfil Salinidad (PSU)', 'Perfil Oxígeno (mg/L)'};

for K = 1 : 3
    ax = subplot(2, 3, K);
    hold(ax, 'on');
    plot(ax, X1{K}, Y1, 'b -');
    plot(ax, X2{K}, Y2, 'r --');
    title(ax, titles{K});
    yticks(0:5:100:
    set(ax, 'xaxisLocation', 'top');
    set(ax, 'YDir', 'reverse');
    ylabel(ax, 'm');
    legend(ax, 'Est.1','Est.2', 'Location','south' );
end

[Felipe Rios]

Works like a charm! Many thanks Walter.