clear
close all
R=input('What is the name of the data file (eg, mbari_design)?  ','s');
eval(['load ', R])

a=find(t==30):length(t);
%a=1:length(t);
dt=t(2)-t(1);
t=t(1:length(a));
Ten=T_t(a,1:3);
mean_ten=round(mean(Ten))
max_ten=mean(Ten)+4*std(Ten)

ay=ceil(max(max(Ten))/1000)*1000;
figure(1)
subplot(311)
plot(t,Ten(:,2))
axis([0 t(end) 0 ay])
title(['Tension below the surface buoy (lbf), mean=',num2str(mean_ten(2)), ...
      ', and max=', num2str(max_ten(2))])
zoom on
subplot(312)
plot(t,Ten(:,3))
axis([0 t(end) 0 ay])
title(['Tension at the junction (lbf), mean=',num2str(mean_ten(3)), ...
      ', and max=', num2str(max_ten(3))])
zoom on
subplot(313)
plot(t,Ten(:,1))
axis([0 t(end) 0 ay])
title(['Tension at the anchor (lbf), mean=',num2str(mean_ten(1)), ...
      ', and max=', num2str(max_ten(1))])
xlabel('Time (s)')
zoom on

figure(2)
plot(x,z)
grid on
title('Static configuration')
ylabel('Depth (m)')
xlabel('Horizontal Coordinate (m)')
zoom on

figure(3)
subplot(211)
plot(s, T)
grid on
title('Static tension (lbf)')
ylabel('Depth (m)')
xlabel('Horizontal Coordinate (m)')
zoom on
subplot(212)
plot(s, phi)
grid on
title('Static angle (rad)')
ylabel('Lagrangian Coordinate (m)')
xlabel('angle (rad)')
zoom on



