function ddb_openDelft3DFLOW(opt)
%DDB_OPENDELFT3DFLOW  One line description goes here.
%
%   More detailed description goes here.
%
%   Syntax:
%   ddb_openDelft3DFLOW(opt)
%
%   Input:
%   opt =
%
%
%
%
%   Example
%   ddb_openDelft3DFLOW
%
%   See also

%% Copyright notice
%   --------------------------------------------------------------------
%   Copyright (C) 2011 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: $
% $Date: $
% $Author: $
% $Revision: $
% $HeadURL: $
% $Keywords: $

%%
handles=getHandles;

switch opt
    case {'opendomains'}
        % DD
        [filename, pathname, filterindex] = uigetfile('*.ddb', 'Select ddbound file');
        if pathname~=0
            pathname=pathname(1:end-1); % Get rid of last file seperator
            if ~strcmpi(pathname,handles.workingDirectory)
                cd(pathname);
                handles.workingDirectory=pathname;
            end
            ddb_plotDelft3DFLOW('delete');
            handles.Model(md).Input=[];
            handles=ddb_readDDBoundFile(handles,filename);
            for i=1:handles.Model(md).nrDomains
                handles.activeDomain=i;
                runid=handles.Model(md).Input(i).runid;
                handles=ddb_initializeFlowDomain(handles,'all',i,runid);
                filename=[runid '.mdf'];
                handles=ddb_readMDF(handles,filename,i);
                handles=ddb_readAttributeFiles(handles,i);
            end
            % Get coordinates of DD boundaries
            handles=ddb_getDDBoundCoordinates(handles);
            handles.activeDomain=1;
            setHandles(handles);
            ddb_plotDelft3DFLOW('plot','active',0,'visible',1,'domain',0);
            gui_updateActiveTab;
        end
    case {'openpresent'}
        % One Domain
        [filename, pathname, filterindex] = uigetfile('*.mdf', 'Select MDF file');
        if pathname~=0
            pathname=pathname(1:end-1); % Get rid of last file seperator
            if ~strcmpi(pathname,handles.workingDirectory)
                cd(pathname);
                handles.workingDirectory=pathname;
            end
            ddb_plotDelft3DFLOW('delete');
            id=handles.activeDomain;
            handles.Model(md).Input=clearStructure(handles.Model(md).Input,id);
            runid=filename(1:end-4);
            handles.Model(md).domains{id}=runid;
            handles.Model(md).DDBoundaries=[];
            handles=ddb_initializeFlowDomain(handles,'all',id,runid);
            filename=[runid '.mdf'];
            handles=ddb_readMDF(handles,filename,id);
            handles=ddb_readAttributeFiles(handles,id);

            xl(1)=min(min(handles.Model(md).Input(id).gridX));
            xl(2)=max(max(handles.Model(md).Input(id).gridX));
            yl(1)=min(min(handles.Model(md).Input(id).gridY));
            yl(2)=max(max(handles.Model(md).Input(id).gridY));
            handles=ddb_zoomTo(handles,xl,yl,0.1);
            
            setHandles(handles);
            ddb_plotDelft3DFLOW('plot','active',0,'visible',1,'domain',ad);
            ddb_updateDataInScreen;
            gui_updateActiveTab;
        end
    case {'opennew'}
        % One Domain
        [filename, pathname, filterindex] = uigetfile('*.mdf', 'Select MDF file');
        if pathname~=0
            pathname=pathname(1:end-1); % Get rid of last file seperator
            if ~strcmpi(pathname,handles.workingDirectory)
                cd(pathname);
                handles.workingDirectory=pathname;
            end
            ddb_plotDelft3DFLOW('delete');
            handles.Model(md).nrDomains=handles.Model(md).nrDomains+1;
            handles.activeDomain=handles.Model(md).nrDomains;
            id=handles.activeDomain;
            handles.Model(md).Input=appendStructure(handles.Model(md).Input);
            runid=filename(1:end-4);
            handles.Model(md).domains{id}=runid;
            handles.Model(md).DDBoundaries=[];
            handles=ddb_initializeFlowDomain(handles,'all',id,runid);
            filename=[runid '.mdf'];
            handles=ddb_readMDF(handles,filename,id);
            handles=ddb_readAttributeFiles(handles,id);

            xl(1)=min(min(handles.Model(md).Input(id).gridX));
            xl(2)=max(max(handles.Model(md).Input(id).gridX));
            yl(1)=min(min(handles.Model(md).Input(id).gridY));
            yl(2)=max(max(handles.Model(md).Input(id).gridY));
            handles=ddb_zoomTo(handles,xl,yl,0.1);
            
            setHandles(handles);
            ddb_plotDelft3DFLOW('plot','active',0,'visible',1,'domain',0);
            ddb_updateDataInScreen;
            gui_updateActiveTab;
        end
end

ddb_refreshDomainMenu;
