% COMPUTE_LIQUID.M %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Compute liquid flows for a section

function L = compute_liquid(secnr,L_in,x,LO,xO)
global coldata hydraulic
i=1:coldata.Nt(secnr);L(i)  = LO + (x(i)-xO(i))*hydraulic.k1 ;
L=L';
