% plot_ctd.m 
% creado por: John Leonard, el 8 de agosto de 1998 para el curso 13.010, otoño 1998
% archivo de comandos de matlab para trazar la salinidad, la densidad y la velocidad del sonido 
% como función de la profundidad a partir de los datos de conductividad, temperatura, 
% y presión.

% este archivo supone que tenemos un array llamado datos,

load('sample.dat');      % carga los datos y los ubica en la variable llamada "sample" (muestra) 
D = sample(:,1);  % la primera columna de datos corresponde a la profundidad
T = sample(:,2);  % la segunda columna de datos corresponde a la temperatura
S = sample(:,3);  % la tercera columna de datos corresponde a la salinidad

PR = 10;          % utilice 10 decibares para la presión de referencia (1 atm)
LAT = 0;          % utilice la latitud cero en todos sus cálculos (el error es mínimo)
P = sw_pres(D,LAT); % calcule la presión

THETA = sw_ptmp(S,T,P,PR);
SIGMA = sw_dens(S,T,P) - 1000,0;         % calcule sigma



subplot(131)
plot(T,D)
hold on
plot(THETA, D, 'r-.')
title('temperatura frente a profundidad')hold on(gca, 'Ydir', 'reverso')
grid on
   (añadir líneas divisorias o rejillas  a la gráfica)
zoom on
labelx('temperatura (deg C)')
labely('profundidad (m)')

subplot(132)
plot(S,D)
title('salinidad frente a profundidad')
grid on
zoom on
labelx('salinidad(PSU)')
labely('profundidad (m)')
hold on(gca, 'Ydir', 'reverse')


subplot(133)plot(SIGMA, D)
hold on(gca, 'Ydir', 'reverso')grid on
zoom on
labelx('\sigma (PSU)')
labely('profundidad (m)')
title('densidad potencial frente a profundidad')