


FD_coefs = ReadData('tcyl3.out');

F = FD_coefs.F;
w = FD_coefs.w;

figure
subplot(2,2,1),plot(w,real(F{1,1}),w,real(F{2,2}),w,real(F{3,3}))
xlabel('w');
ylabel('Added Mass');
legend('Surge','Sway','Heave');


subplot(2,2,2),plot(w,imag(F{1,1}),w,imag(F{2,2}),w,imag(F{3,3}))
xlabel('w');
ylabel('Damping');
legend('Surge','Sway','Heave');


subplot(2,2,3),plot(w,real(F{7,1}),w,real(F{7,2}),w,real(F{7,3}));
xlabel('w');
ylabel('Re[Exciting Force]');
legend('Surge','Sway','Heave');


subplot(2,2,4),plot(w,imag(F{7,1}),w,imag(F{7,2}),w,imag(F{7,3}));
xlabel('w');
ylabel('Im[Exciting Force]');
legend('Surge','Sway','Heave');