s1 = 64850;
s2 = 65050;
ii = s1:s2;
filler = ones(size(c_time));
filler = filler *0;
ticker = c_time(s1:24:s2);


figure;
subplot(4,1, 1);
plot(c_time(ii), c_speed(ii));
ylabel('Current (m/s)');
set(gca, 'Xlim', [c_time(s1) c_time(s2)]);
set(gca,'Xtick',ticker)
set(gca,'XTickLabel',datestr(ticker,2))

subplot(4,1, 2);
plot(wind_time(ii), wind_speed(ii));
ylabel('Wind (m/s)');
set(gca, 'Xlim', [c_time(s1) c_time(s2)]);
set(gca,'Xtick',ticker)
set(gca,'XTickLabel',datestr(ticker,2))

subplot(4,1, 3);
plot(wind_time(ii), power(ii));
hold on;
plot(wind_time(ii), wind_power(ii)*effg, 'r');
hold off
ylabel('Power (W)');
set(gca, 'Xlim', [c_time(s1) c_time(s2)]);
legend('Power Needed','Power Generated');
set(gca,'Xtick',ticker)
set(gca,'XTickLabel',datestr(ticker,2))
    
subplot(4,1, 4);
quiver(c_time(ii), filler(ii), fu(ii), fv(ii));
hold on;
quiver(c_time(ii), filler(ii), u_20m(ii), v_20m(ii), 'r');
quiver(c_time(ii), filler(ii), wind_u(ii), wind_v(ii), 'g');
axis equal;
set(gca, 'Xlim', [c_time(s1) c_time(s2)]);
set(gca,'Xtick',ticker)
set(gca,'XTickLabel',datestr(ticker,2))
legend('Total force on buoy','Current Direction', 'Wind Direction');


