function [Ero Sed t] = decc_ero_sed_time(varargin)
%DECC_ERO_SED_TIME  One line description goes here.
%
%   More detailed description goes here.
%
%   Syntax:
%   varargout = decc_ero_sed_time(varargin)
%
%   Input:
%   varargin  =
%
%   Output:
%   varargout =
%
%   Example
%   decc_ero_sed_time
%
%   See also

%% Copyright notice
%   --------------------------------------------------------------------
%   Copyright (C) 2012 Delft University of Technology
%       Kees den Heijer
%
%       C.denHeijer@TUDelft.nl
%
%       Faculty of Civil Engineering and Geosciences
%       P.O. Box 5048
%       2600 GA 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: 07 Feb 2012
% Created with Matlab version: 7.12.0.635 (R2011a)

% $Id$
% $Date$
% $Author$
% $Revision$
% $HeadURL$
% $Keywords: $

%%
OPT = struct(...
...    'ncfile', 'F:\Sources\OpenEarthModels\tudelft\dune_erosion_coastline_curvature\simulations\run24\radius01900\rev2585M\xboutput.nc',...
    'ncfile', 'd:\heijer\sources\dune_erosion_coastline_curvature\simulations\00000000\xbeach\run26\radius01900\xboutput.nc',...
    'angle', 0);

OPT = setproperty(OPT, varargin);

%%
ncfile = OPT.ncfile;
t = nc_varget(ncfile, 'globaltime');
y = nc_varget(ncfile, 'globaly')';
x = nc_varget(ncfile, 'globalx')';

g = xb_stagger(x, y);

phi = -rad2deg(mean(g.alfaz,1));
id1 = find(phi < OPT.angle, 1, 'last');
id = round(interp1(phi(id1+[0 1]), id1+[0 1], OPT.angle));

dsz = g.dsz;
dsz(1,:) = 0;
s = cumsum(dsz, 1);
s = s(:,id);

zb = nc_varget(ncfile, 'zb', [0 id-1 0], [-1 1 -1]);
zb0 = repmat(zb(1,:), size(zb,1), 1);
svp = ones(size(zb)) * 5;
zb1 = cat(3, zb, svp);

dsdnz = repmat(g.dsdnz(:,id)', size(zb,1), 1);
dnz = repmat(g.dnz(:,id)', size(zb,1), 1);

Ero = sum((max(zb1, [], 3) - max(zb0, 5)) .* dsdnz ./ dnz, 2);
Sed = sum((min(zb1, [], 3) - min(zb0, 5)) .* dsdnz ./ dnz, 2);