r_figs_summary.m 2.01 KB
function r_figs_summary


lw = 1.5;
ms = 10;

tp = [6.5 8.0 9.5 11.0 12.5 14.0 15.5 17.0 18.5 20.0];

%% 80 floes
dy_hs2  = [-0.026 -0.007  NaN NaN    NaN    NaN    NaN   NaN   NaN    NaN];
dy_hs4  = [   NaN -0.052  NaN -0.016 -0.017 -0.009 -0.005 0.000 0.000 -0.003];
dy_hs8  = [   NaN    NaN  NaN NaN    NaN -0.016    NaN   NaN   NaN    NaN];
dy_hs10 = [   NaN    NaN  NaN NaN    NaN    NaN    NaN   NaN   NaN  0.003];

figure(1); clf
hold on
plot(tp,-100.*dy_hs2,'+k','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs4,'ok','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs8,'xk','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs10,'sqk','MarkerSize',ms,'LineWidth',lw);
plot(tp,0.*tp,'--k');
hold off
set(gca,'FontName','Times New Roman','FontSize',18);
xlabel('\itT_p \rm(s)','FontName','Times New Roman','FontSize',18);
ylabel('(\Delta\ity\rm)_{ave} (m)','FontName','Times New Roman','FontSize',18);
legend('\itH_s \rm= 2m','\itH_s \rm= 4m','\itH_s \rm= 8m','\itH_s \rm= 10m');
ylim([-1 6])

%% 40 floes
dy_hs2  = [-0.013    -0.010  NaN NaN    NaN    NaN    NaN   NaN   NaN    NaN];
dy_hs3  = [   NaN    NaN  -0.016 NaN    NaN    NaN    NaN NaN NaN NaN];
dy_hs4  = [   NaN    NaN  NaN NaN 0.000 NaN   -0.002   NaN   0.000    NaN];
dy_hs6  = [   NaN    NaN  -0.082 NaN  NaN    NaN    NaN   NaN   NaN    NaN];
dy_hs8 = [   NaN    NaN  NaN NaN    NaN   -0.008    NaN   NaN   NaN    NaN];

figure(2); clf
hold on
plot(tp,-100.*dy_hs2,'+k','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs3,'ok','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs4,'xk','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs6,'sqk','MarkerSize',ms,'LineWidth',lw);
plot(tp,-100.*dy_hs8,'dk','MarkerSize',ms,'LineWidth',lw);
plot(tp,0.*tp,'--k');
hold off
set(gca,'FontName','Times New Roman','FontSize',18);
xlabel('\itT_p \rm(s)','FontName','Times New Roman','FontSize',18);
ylabel('(\Delta\ity\rm)_{ave} (m)','FontName','Times New Roman','FontSize',18);
legend('\itH_s \rm= 2m','\itH_s \rm= 3m','\itH_s \rm= 4m','\itH_s \rm= 6m','\itH_s \rm= 8m');
ylim([-1 6])