%  overplot_ship and ao 
%% load data
load ship.mat
load ../Argos.mat
% load ../sonic1m.mat

%% truncate for better view
ix=find (ship.mday > sys1.mday(1) & ship.mday <= sys1.mday(end));
ship=struct_trunc(ship,ix);
ix=find (sys1.mday > ship.mday(1) & sys1.mday <= ship.mday(end));
sys1=struct_trunc(sys1,ix);
ix=find (sys2.mday > ship.mday(1) & sys2.mday <= ship.mday(end));
sys2=struct_trunc(sys2,ix);


%% Now plot 

% 1.atmp hrh  2.bpr  3. wind
% ship plot 'bpr_c'    'hrh'    'atmp'    'wdir'    'wspd'  
figure(1); clf; 
subplot(2,1,1);    %atmp
plot(ship.mday,ship.atmp,sys1.mday,sys1.atmp,sys2.mday,sys2.atmp )
tlabel
ylabel('atmp');

subplot(2,1,2);  %hrh
plot(ship.mday,ship.hrh,sys1.mday,sys1.hrh,sys2.mday,sys2.hrh )
tlabel
ylabel('hrh');

% 2. bpr - and bpr_c
figure(2); clf; 
subplot(2,1,1);    %bpr
plot(ship.mday,ship.bpr,sys1.mday,sys1.bpr,sys2.mday,sys2.bpr )
tlabel
ylabel('bpr');

subplot(2,1,2);  %bpr_c
plot(ship.mday,ship.bpr_c,sys1.mday,sys1.bpr,sys2.mday,sys2.bpr )
tlabel
ylabel('bpr(c)');

% 3. wind: speed, dir???  all uncorrected???
figure(3);clf;  % wind spd & dir

subplot(2,1,1);   % wind speed wnds
% ship 'wdir,wspd'  lgrs wndd wnds
plot(ship.mday,ship.wspd,sys1.mday,sys1.wnds,sys2.mday,sys2.wnds)
ylabel('relwspd');
tlabel
 
subplot(2,1,2);  %wdir
plot(ship.mday,ship.wdir,sys1.mday,sys1.wndd,sys2.mday,sys2.wndd)
ylabel('relwdir');
tlabel

