function plot_strandbreedte_tijd(raai)

%% Script voor plotten strandbreedte in tijd
%
% Auteur:   Vincent Vuik
% Bedrijf:  HKV lijn in water
% Project:  PR2063.20 - Kustlijnzorg
% Datum:    19 mei 2011
%
% Raai: PC-ring codering, bijv. 7002300 (kustvak 7, km 23.00)
%                         bijv. 13001250 (kustvak 13, km 12.50)

savefigures    = 'ja';    % ja/nee
closefigures   = 'ja';    % ja/nee
figure_visible = 'nee';   % ja/nee (nee gaat sneller, dan wel savefigures='ja', anders plotfigures zinloos)

maindir    = 'g:\Projecten\pr2063.20 - Kustlijnzorg\Strandbreedte\';
outputdir  = [maindir,'output\'];
figdir     = [maindir,'figuren\strandbreedte_tijd\beide\'];

cd(outputdir)
load('strandbreedtes_droog_strand.mat');

element                                                 = find(raaien==raai);
strandbreedte_plot_droog                                = strandbreedtes(element,:);
strandbreedte_plot_droog(strandbreedte_plot_droog==0)   = NaN;
duingrens_plot                                          = -duingrenzen  (element,:);
duingrens_plot          (strandbreedte_plot_droog==0)   = NaN;
zeegrens_plot_droog                                     = -zeegrenzen   (element,:);
zeegrens_plot_droog     (strandbreedte_plot_droog==0)   = NaN;

load('strandbreedtes_totaal_strand.mat');

element                                                 = find(raaien==raai);
strandbreedte_plot_totaal                               = strandbreedtes(element,:);
strandbreedte_plot_totaal(strandbreedte_plot_totaal==0) = NaN;
zeegrens_plot_totaal                                    = -zeegrenzen   (element,:);
zeegrens_plot_totaal     (strandbreedte_plot_totaal==0) = NaN;

ylims(1) = 0;
ylims(2) = 10*round((max(strandbreedte_plot_totaal)+15)/10);
xlims(1) = 1960;
xlims(2) = 2015;

if strcmp(figure_visible,'nee')==1
    figure('visible','off');
else
    figure(1);
end

plot(jaren,strandbreedte_plot_droog,'LineStyle','-','Color','r','Marker','s','MarkerEdgeColor','r','MarkerFaceColor','y','MarkerSize',4);
hold on;
plot(jaren,strandbreedte_plot_totaal,'LineStyle','-','Color','b','Marker','s','MarkerEdgeColor','b','MarkerFaceColor','y','MarkerSize',4);
hold off;
grid on;
xlabel('tijd [jaren]');
ylabel('strandbreedte [m]');
title(['Verloop van de strandbreedte in de tijd voor raai ',num2str(raai)]);
ylim(ylims);
xlim(xlims);
set(gca,'box','on','XTick',[xlims(1):5:xlims(2)]);
legend('Strand bij GHW','Strand bij GLW');

if strcmp(savefigures,'ja')==1
    figname = [figdir,'strandbreedte_tijd_raai_',num2str(raai),'.png'];
    print(figname,'-dpng');
end

ylims(1) = 10*round((nanmin(zeegrens_plot_totaal)-5 )/10);
ylims(2) = 10*round((nanmax(duingrens_plot      )+35)/10);

if strcmp(figure_visible,'nee')==1
    figure('visible','off');
else
    figure(2);
end

plot(jaren,duingrens_plot,'LineStyle','-','Color','r','Marker','s','MarkerEdgeColor','r','MarkerFaceColor','y','MarkerSize',4);
hold on;
plot(jaren,zeegrens_plot_droog ,'LineStyle','-','Color','b','Marker','s','MarkerEdgeColor','b','MarkerFaceColor','y','MarkerSize',4);
plot(jaren,zeegrens_plot_totaal ,'LineStyle','-','Color','g','Marker','s','MarkerEdgeColor','g','MarkerFaceColor','y','MarkerSize',4);
hold off;
grid on;
xlabel('tijd [jaren]');
ylabel('Oriëntatie t.o.v. RSP-lijn [m]');
title(['Verloop van de strandbegrenzingen in de tijd voor raai ',num2str(raai)]);
ylim(ylims);
xlim(xlims);
set(gca,'box','on','XTick',[xlims(1):5:xlims(2)]);
legend('Duingrens','Zeegrens bij GHW','Zeegrens bij GLW');

if strcmp(savefigures,'ja')==1
    figname = [figdir,'strandbegrenzingen_tijd_raai_',num2str(raai),'.png'];
    print(figname,'-dpng');
end

if strcmp(closefigures,'ja')==1
    close all;
end

cd(maindir);
