function runStrandbreedte(kustvak,jaren,strandtype)

%% Script voor de communicatie met bepaling_strandbreedte.m
%
% Auteur:   Vincent Vuik
% Bedrijf:  HKV lijn in water
% Project:  PR2063.20 - Kustlijnzorg
% Datum:    19 mei 2011
%
% Dit script stuurt de juiste raaien en jaren naar het script
% bepaling_strandbreedte.m. Daarnaast worden alle strandbreedtes
% gecombineerd weggeschreven in een mat-file.

close all;
clc;

%% Locaties en instellingen

kustvaknr  = str2num(kustvak(1:2));
opslaan    = 'ja';

maindir    = 'g:\Projecten\pr2063.30 - Kustlijnzorg\Werkmap\';
datadir    = [maindir,'Profielen\naStap01\kustvak',n2svol(kustvaknr,2),'\'];
inputdir   = [maindir,'Doorgerekende vakken\'];
outputdir  = [maindir,'Strandbreedte\output\'];
xlsfile    = [outputdir,'strandbreedtes_',strandtype,'.xls'];

% Lijst met te analyseren raaien
lijst_file = [inputdir,'vakken_',kustvak,'.txt'];

%% Bepalen en wegschrijven strandbreedtes

raaien      = load(lijst_file);

%waitbar(0,strrep(kustvak,'_',''));
for r = 1:length(raaien)

    raai = raaien(r);

    % inlezen mat-files met strandbreedtes indien al aanwezig
    if strcmp(strandtype,'droog')==1 && r==1 && exist([outputdir,'strandbreedtes_droog_strand.mat'],'file')
        load([outputdir,'strandbreedtes_droog_strand.mat']);
    elseif strcmp(strandtype,'totaal')==1 && r==1 && exist([outputdir,'strandbreedtes_totaal_strand.mat'],'file')
        load([outputdir,'strandbreedtes_totaal_strand.mat']);
    end

    for j = 1:length(jaren)
        jaar = jaren(j);
        % bepaal of er voor het betreffende jaar data bestaat
        if exist([datadir,num2str(jaar),'.txt'],'file')==2
            % stuur script bepaalStrandbreedte.m aan
            [strandbreedte,duingrens,zeegrens] = ...
                bepaalStrandbreedte(maindir,datadir,raai,jaar,strandtype);

            strandbreedtes(r,j) = strandbreedte;
            duingrenzen(r,j)    = duingrens;
            zeegrenzen(r,j)     = zeegrens;
        else % geen file met data voor het betreffende jaar
            strandbreedtes(r,j) = NaN;
            duingrenzen(r,j)    = NaN;
            zeegrenzen(r,j)     = NaN;
        end
    end

    % opslaan mat-files met matrices, dimensies [raaien,jaren]
    if strcmp(opslaan,'ja')==1     && strcmp(strandtype,'droog')==1  && r==length(raaien)
        save([outputdir,'strandbreedtes_droog_',kustvak,'.mat'], 'strandbreedtes','duingrenzen','zeegrenzen','jaren','raaien');
        xlswrite(xlsfile,strandbreedtes,kustvak,'B2');
        xlswrite(xlsfile,jaren,kustvak,'B1');
        xlswrite(xlsfile,raaien,kustvak,'A2');
    elseif strcmp(opslaan,'ja')==1 && strcmp(strandtype,'totaal')==1 && r==length(raaien)
        save([outputdir,'strandbreedtes_totaal_',kustvak,'.mat'],'strandbreedtes','duingrenzen','zeegrenzen','jaren','raaien');
        xlswrite(xlsfile,strandbreedtes,kustvak,'B2');
        xlswrite(xlsfile,jaren,kustvak,'B1');
        xlswrite(xlsfile,raaien,kustvak,'A2');
    end
    %waitbar(r/length(raaien));
end

close force all;