function delft3dfm_save_boundary_ext_file(extfile,boundary)
%ddb_DFlowFM_saveExtForcing

%% Copyright notice
%   --------------------------------------------------------------------
%   Copyright (C) 2013 Deltares
%       Maarten van Ormondt
%
%       Maarten.vanOrmondt@deltares.nl
%
%       P.O. Box 177
%       2600 MH Delft
%       The Netherlands
%
%   This library is free software: you can redistribute it and/or modify
%   it under the terms of the GNU General Public License as published by
%   the Free Software Foundation, either version 3 of the License, or
%   (at your option) any later version.
%
%   This library is distributed in the hope that it will be useful,
%   but WITHOUT ANY WARRANTY; without even the implied warranty of
%   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%   GNU General Public License for more details.
%
%   You should have received a copy of the GNU General Public License
%   along with this library.  If not, see <http://www.gnu.org/licenses/>.
%   --------------------------------------------------------------------

% This tool is part of <a href="http://www.OpenEarth.eu">OpenEarthTools</a>.
% OpenEarthTools is an online collaboration to share and manage data and
% programming tools in an open source, version controlled environment.
% Sign up to recieve regular updates of this function, and to contribute
% your own tools.

%% Version <http://svnbook.red-bean.com/en/1.5/svn.advanced.props.special.keywords.html>
% Created: 29 Nov 2011
% Created with Matlab version: 7.11.0.584 (R2010b)

% $Id: ddb_DFlowFM_writeComponentsFile.m 9233 2013-09-19 09:19:19Z ormondt $
% $Date: 2013-09-19 11:19:19 +0200 (Thu, 19 Sep 2013) $
% $Author: ormondt $
% $Revision: 9233 $
% $HeadURL: https://svn.oss.deltares.nl/repos/openearthtools/trunk/matlab/applications/DelftDashBoard/models/DFlowFM/fileio/ddb_DFlowFM_writeComponentsFile.m $
% $Keywords: $
    
% Saves Delft3D-FM boundary ext file

nb=0;
for ib=1:length(boundary)
    
    bnd=boundary(ib).boundary;
    
    switch bnd.type
        case{'water_level'}
            iz=1;
            in=0;
            it=0;
        case{'water_level_plus_normal_velocity'}
            iz=1;
            in=1;
            it=0;
        case{'water_level_plus_normal_velocity_plus_tangential_velocity'}
            iz=1;
            in=1;
            it=1;
    end
    
    if iz==1
        
        % waterlevelbnd
        nfrc=0;
        frcfil={''};
        
        if bnd.water_level.time_series.active
            fil=bnd.water_level.time_series.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.water_level.astronomic_components.active
            fil=bnd.water_level.astronomic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.water_level.harmonic_components.active
            fil=bnd.water_level.harmonic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        for jj=1:length(frcfil)
            nb=nb+1;
            b(nb).quantity='waterlevelbnd';
            b(nb).locationfile=bnd.location_file;
            b(nb).forcingfile=frcfil{jj};
        end
        
    end
    
    if in==1
        
        % waterlevelbnd
        nfrc=0;
        frcfil={''};
        
        if bnd.normal_velocity.time_series.active
            fil=bnd.normal_velocity.time_series.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.normal_velocity.astronomic_components.active
            fil=bnd.normal_velocity.astronomic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.normal_velocity.harmonic_components.active
            fil=bnd.normal_velocity.harmonic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        for jj=1:length(frcfil)
            nb=nb+1;
            b(nb).quantity='normalvelocitybnd';
            b(nb).locationfile=bnd.location_file;
            b(nb).forcingfile=frcfil{jj};
        end
        
    end
    
    if it==1
        
        % waterlevelbnd
        nfrc=0;
        frcfil={''};
        
        if bnd.tangential_velocity.time_series.active
            fil=bnd.tangential_velocity.time_series.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.tangential_velocity.astronomic_components.active
            fil=bnd.tangential_velocity.astronomic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        if bnd.tangential_velocity.harmonic_components.active
            fil=bnd.tangential_velocity.harmonic_components.forcing_file;
            ifil=strmatch(lower(fil),lower(frcfil));
            if isempty(ifil)
                nfrc=nfrc+1;
                frcfil{nfrc}=fil;
            end
        end
        
        for jj=1:length(frcfil)
            nb=nb+1;
            b(nb).quantity='tangentialvelocitybnd';
            b(nb).locationfile=bnd.location_file;
            b(nb).forcingfile=frcfil{jj};
        end
        
    end
    
end

% Boundaries in new ext
fid=fopen(extfile,'wt');
for ip=1:length(b)
    fprintf(fid,'%s\n','[boundary]');
    fprintf(fid,'%s\n',['quantity       = ' b(ip).quantity]);
    fprintf(fid,'%s\n',['locationfile   = ' b(ip).locationfile]);
    fprintf(fid,'%s\n',['forcingfile    = ' b(ip).forcingfile]);
    fprintf(fid,'%s\n','');
end
fclose(fid);
