

figure('Position',[227    412    951   1027],'paperorientation','portrait','papersize',[8.5 11]);
ah1 = subplot(3,1,1)
plot(WindAndWaterResults.t,WindAndWater.seacurrentX(:,2))
hold all
plot(WindAndWaterResults.t,WindAndWater.seacurrentX(:,3))
plot(WindAndWaterResults.t,WindAndWater.seacurrentX(:,4))
legend('z = 0m','z = -5m ','z = -10m')
xlabel('time (s)');
ylabel('X Sea Velocity (m/s)');


ah2 = subplot(3,1,2)
plot(WindAndWaterResults.t,WindAndWater.seacurrentY(:,2))
hold all
plot(WindAndWaterResults.t,WindAndWater.seacurrentY(:,3))
plot(WindAndWaterResults.t,WindAndWater.seacurrentY(:,4))
legend('z = 0m','z = -5m ','z = -10m')
xlabel('time (s)');
ylabel('Y Sea Velocity (m/s)');


ah3 = subplot(3,1,3)
plot(WindAndWaterResults.t,WindAndWater.seacurrentZ(:,2))
hold all
plot(WindAndWaterResults.t,WindAndWater.seacurrentZ(:,3))
plot(WindAndWaterResults.t,WindAndWater.seacurrentZ(:,4))
legend('z = 0m','z = -5m ','z = -10m')
xlabel('time (s)');
ylabel('Z Sea Velocity (m/s)');


print('-dpdf','SeaVelocity.pdf');


figure
plot(MooringRiserResults.Y',-MooringRiserResults.Z');
hold all
set(gca,'ColorOrderIndex',1);
plot(HoseTetherResults.Y',-HoseTetherResults.Z');

lh = plot(MooringRiserResults.MeanY',-MooringRiserResults.MeanZ','-k');
set(lh,'linewidth',3);
lh = plot(HoseTetherResults.MeanY',-HoseTetherResults.MeanZ','-k');
set(lh,'linewidth',3);
