function [GHW,GLW] = bepaalGHWGLW(maindir,raai)

%% Script voor bepalen GHW- en GLW-stand voor een bepaalde raai
%
% Programmer: Vincent Vuik
% Company:    HKV lijn in water
% Project:    PR2063.20 - Kustlijnzorg
% Date:       19 mei 2011
%
% Dit script wordt aangestuurd vanuit bepaalStrandbreedte.m
% Raai: PC-ring codering, bijv. 7002300 (kustvak 7, km 23.00)
%                         bijv. 13001250 (kustvak 13, km 12.50)

%% Locatie benodigde file

GHW_GLW_file  = [maindir,'Strandbreedte\hw-lw-lijn\gemiddeld hoog en laag water 2008 edit.txt'];

%% Openen file

GHW_GLW_fid   = fopen(GHW_GLW_file,'r');

%% Bestand met GLW- en GHW-niveau langs de kust inlezen

% kolommen in file: 
% - kustvak
% - DONAR-codering kustvak
% - beschrijving kustvak
% - gehele kilometer
% - GHW [cm]
% - GLW [cm]

InputText     = textscan(GHW_GLW_fid,'%n %s %s %n %n %n','delimiter',' ','multipleDelimsAsOne',1,'headerlines',4);
kustvak       = InputText{1,1};
metingnaam    = InputText{1,3};
metinglocatie = InputText{1,4};
GHW_niveau    = InputText{1,5};
GLW_niveau    = InputText{1,6};

%% Bepalen GHW en GLW bij opgevraagde raainummer

kustvak_raai     = floor(raai/1e6);                 
km_raai          = (raai-kustvak_raai*1e6)/100;     

elements_kustvak = find(kustvak==kustvak_raai);
km_kustvak       = metinglocatie(elements_kustvak);
GHW_kustvak      = GHW_niveau(elements_kustvak);
GLW_kustvak      = GLW_niveau(elements_kustvak);


GHW              = interp1(km_kustvak,GHW_kustvak,km_raai,'linear','extrap')/100;   % cm --> m+NAP
GLW              = interp1(km_kustvak,GLW_kustvak,km_raai,'linear','extrap')/100;   % cm --> m+NAP

%% Sluiten file

close_GHW_GLW_fid = fclose(GHW_GLW_fid);