load static_minus
depth=4000;
Z_minus=Z_minus-depth;
Z_minus2=Z_minus2-depth;
Z_minus3=Z_minus3-depth;
Z_minus4=Z_minus4-depth;
Z_design=Z_design-depth;
ns=nodes(end);
figure(1)
clf
%plot(X_design,Z_design,'linewidth',2.5)
hold
y=Z_design;
y(ns+1:end)=y(ns)+(y(ns+1:end)-y(ns))*y(ns)/(y(ns)-y(end));
%plot(X_design(end),y(end),'s','MarkerSize',10,'LineWidth',2.5);
%plot(X_design(ns),y(ns),'o','MarkerSize',8,'LineWidth',2.5);

plot(X_minus,Z_minus,'linewidth',2.5)
y=Z_minus;
y(ns+1:end)=y(ns)+(y(ns+1:end)-y(ns))*y(ns)/(y(ns)-y(end));
plot(X_minus(end),y(end),'s','MarkerSize',10,'LineWidth',2.5);
plot(X_minus(ns),y(ns),'o','MarkerSize',8,'LineWidth',2.5);

plot(X_minus2,Z_minus2,'m','linewidth',2.5)
y=Z_minus2;
y(ns+1:end)=y(ns)+(y(ns+1:end)-y(ns))*y(ns)/(y(ns)-y(end));
plot(X_minus2(end),y(end),'s','MarkerSize',10,'LineWidth',2.5);
plot(X_minus2(ns),y(ns),'o','MarkerSize',8,'LineWidth',2.5);

plot(X_minus3,Z_minus3,'g','linewidth',2.5)
y=Z_minus3;
y(ns+1:end)=y(ns)+(y(ns+1:end)-y(ns))*y(ns)/(y(ns)-y(end));
plot(X_minus3(end),y(end),'s','MarkerSize',10,'LineWidth',2.5);
plot(X_minus3(ns),y(ns),'o','MarkerSize',8,'LineWidth',2.5);

plot(X_minus4,Z_minus4,'r','linewidth',2.5)
y=Z_minus4;
y(ns+1:end)=y(ns)+(y(ns+1:end)-y(ns))*y(ns)/(y(ns)-y(end));
plot(X_minus4(end),y(end),'s','MarkerSize',10,'LineWidth',2.5);
plot(X_minus4(ns),y(ns),'o','MarkerSize',8,'LineWidth',2.5);

plot([0 1700],[0 0],'--')
hold

xlabel('Horizontal Offset (m)','fontsize',12,'fontangle','italic');
ylabel('Depth (m)','fontsize',12,'fontangle','italic');
axis([400 1600 -600 30])
title('4000-Meter Benthic Mooring (w/ 25-Meter Snubber & 55 Meter Pull-Up) ')
grid

figure(2)
clf

Tb=[T_minus(end) T_minus2(end) T_minus3(end) T_minus4(end)];
Xb=[X_minus(end) X_minus2(end) X_minus3(end) X_minus4(end)];

Ta=[T_minus(1) T_minus2(1) T_minus3(1) T_minus4(1)];
Xa=Xb;

plot(Xb,Tb,'k','LineWidth',2.5);
hold
plot(Xb,Tb,'d','MarkerSize',11,'LineWidth',2.5);
plot(Xa,Ta,'k','LineWidth',2.5);
plot(Xa,Ta,'s','MarkerSize',10,'LineWidth',2.5);
hold
xlabel('Horizontal Offset (m)','fontsize',12,'fontangle','italic');
ylabel('Tension (lbs)','fontsize',12,'fontangle','italic');
text(Xb(4)+20,Tb(4),'Surface Buoy','fontsize',10)
text(Xb(4)+20,Ta(4),'Anchor','fontsize',10)
axis([1000 1600 0 5000])




