close all
clear all

basedir = '../HamiltonRuns/';

configdir = 'ASIS_Mooring_80_d110/';

Depth = 110;
Wind = 0;
WaveAmp = 0;

currents = [.50 .60 .70 .8 .9 1.0];


filename = 'ResultsSummary.mat';

for i = 1:length(currents)
dir = [basedir configdir 'Results_c' num2str(100*currents(i)) 'u_d' num2str(Depth) '_w' num2str(Wind) '_H' num2str(WaveAmp) '/'];
load([dir filename]);
 
AnchorHorizForce(i) =  MooringRiserResults.NodeNHorizForce(end)/4.45;
AnchorVertForce(i) = -MooringRiserResults.NodeNVertForce(end)/4.45;

end

figure
plot(currents,AnchorHorizForce,'-x');
hold all
plot(currents,AnchorVertForce,'-x');
xlabel('Current Speed (m/s)');
ylabel('Anchor Force (lbs)');
title([configdir '  Anchor Loads:  Depth = ' num2str(Depth) 'm   Wind = ' num2str(Wind) 'kts'],'interpreter','none');
grid on;
legend('Horizontal Anchor Force','Vertical Anchor Force','Location','NorthWest');
print([basedir configdir 'StaticAnchorLoads.pdf'],'-landscape','-bestfit');
