%% Script to create plots similar to figures 10 and 11 of the report

%clear all;
%close all;
%clear

% choose different transects, years

test_number = 240; %numbers from 1 to 334
year_test   = 17;

visual
tr          = id_DF(test_number); 
yy          = year_test;
[x_df, z_df, x, z, z_mean_high] = dunefoot_plot(tr,yy);  
cell        = find(idID_STb==test_number)
yy          = year_test;

%% Figure    
figure()
set(gcf,'Position',[5 35 1040 760]); 
set(gca,'FontSize',14); grid off
plot(x, squeeze(z(year_test,:,:)),'b'); 
hold on
plot([x(1) x(end)],[z_mean_high z_mean_high],'--c');

df = plot(x_df(year_test),z_df(year_test),'bs','MarkerFaceColor','r','MarkerSize',6);

xlabel('x [m]','FontSize',14); ylabel('z[m]','FontSize',14)
xlim  ([-1000 1000]);
ylim  ([-2 20]);
xlabel('crossshore [m]')
ylabel('z [m]')
hold on    
df_vis=plot([DF_cross(cell,year_test) DF_cross(cell,year_test)],[-3 14],'--r');  
legend('dune profile', 'low tide', '2nd derivative method','visual observation','Location', 'NorthEast','FontSize',14)
