% plot_ts.m - show t/s plots for CTD data vs sled data

cmax = 7; % last CTD cast number
smax = 6; % last sample tow number
% note: cast 8 was outside of the eddy
figure(cmax+1)
clf

for ii = [1:3 5:cmax]
  
  if ii <10
    fn = ['oc41540',int2str(ii)];
  elseif cast <100
    fn = ['oc4154',int2str(ii)];
  end

  eval(['load ',mat_path,fn])

  h1=plot(s0,t0,'.');
  if ~ishold, hold on; end
  
end

for sam=1:smax
   if sam <10
    fn = ['SAM_00',int2str(sam),'_fl'];
  elseif cast <100
    fn = ['SAM_0',int2str(sam),'_fl'];
  end
 
  eval(['load /mnt/oc415-4a/oc415_4/sled/',fn])
  
  h2=plot(sal0,temp0,'r.');
keyboard

end

axis([36.68 36.715 19.15 19.5])
xlabel('Salinity (PSU)')
ylabel('Temperature (^oC)')
legend([h1 h2],'CTD','Sled')
title('OC415-4 T/S for CTD1:7 and Sled tow 1:6')
