function [e,w,loc,hc,time]=xsect(e,w,loc,hc,time,varargin);


  ec=0;
if(time ~= 0.0d0)go to 11; end;
ec = 0.0d0;
e = w.*ec;
% continue;
end

%  FLOW RESISTANCE RELATIONSHIP
%(MANNING)

function [ch,deep,d]=resist(ch,deep,d,varargin);


  n=0;
n = 0.03d0;
ch =(deep.^0.166)./n;
end

%  SEDIMENT TRANSPORT RELATIONSHIP
%(ENGELUND BEDLOAD EQN.)

function [qs,u,deep,ch,rho,rhos,d,tauc]=strans(qs,u,deep,ch,rho,rhos,d,tauc,varargin);


 kn1=0;
kn2=0;
udir=0;
tau=0;
thr1=0;
thr2=0;
g=0;

g = 9.8d0;
kn1 = 5.0d0.*d.*sqrt((rhos-rho).*g.*d./rho);
kn2 = rho./((rhos-rho).*d.*ch.^2);
if(u ~= 0)udir = u./(abs(u)); end;
if(u == 0)udir = 1; end;
tau = u.^2.*kn2;
thr1 = tau - tauc;
if(thr1 <= 0)thr1=0.0d0; end;
if(thr1 > 0)thr1 =(1.0d0+(0.267d0./thr1).^4).^(-0.25); end;
thr2 = sqrt(tau)-  0.7d0.*sqrt(tauc);
if(thr2 <= 0)thr2=0.0d0; end;
qs = udir.*kn1.*thr1.*thr2;
end

