% mk_mat_pod.m

% for deployment 1 - hardwired start time 7/26(yd=206) 2100Z
st_time=206+21/24;
dt=.5/60/60/24;

figure(1)
% RBR DR-1050  3.50 009791
% get start time and interval

% this produces a vector of data
!tail +15 s1_009791.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(1)=plot(ttt,xxx,'b');
hold on

% RBR DR-1050  3.50 009792 
% get start time and interval

% this produces a vector of data
!tail +15 s1_009792.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(2)=plot(ttt,xxx,'g');

legend(h,'9791','9792')

% Temp.
figure(2)
clf

% RBR TR-1050  3.50 010983 (Windows: 5.14)
% get start time and interval

% this produces a vector of data
!tail +15 s1_010983.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(1)=plot(ttt,xxx,'b');
hold on

% RBR TR-1050  3.50 010984
% get start time and interval

% this produces a vector of data
!tail +15 s1_010984.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(2)=plot(ttt,xxx,'g');


% RBR TR-1050  3.50 010986
% get start time and interval

% this produces a vector of data
!tail +15 s1_010986.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(3)=plot(ttt,xxx,'r');


% RBR TR-1050  3.50 011177
% get start time and interval

% this produces a vector of data
!tail +15 s1_011177.dat >xxx
load xxx
ttt=[0:length(xxx)-1]*dt+st_time;
h(4)=plot(ttt,xxx,'c');


legend(h,'10983','10984','10986','11177')

xlabel('Yearday 2005')
ylabel('Temperature (^oC)')
title(['Tpod data from Sample cast: ',int2str(cast)])