%partikkelmodel
clear
%Indata

dt=30;        %tidsskridt
Tmax=6690;  %varighed af simulering i sek.
Q=0.3423;        %vandføring (m3/s)
Areal=2.06;    %m2
D=0.7;          %m2/s
Nmax=20000;   %antal partikler
Masse=0.00225;    %kg masse af stof tilsat
Masse_p=Masse/Nmax;    %massen af en partikkel
k1=Q*dt/Areal;
k2=2*sqrt(6*D*dt);

%Strækning hvor koncentrationen ønskes beregnet
st_ovre=450;
st_nedre=490;

%startbetingelser
step=1;
x=zeros(Nmax,1);

%tidsløkke
while step<Tmax/dt
    for n=1:Nmax
        x(n)=x(n)+k1+k2*(rand-0.5);
    end
    antal=find(x>st_ovre & x<st_nedre);
    c(step)=size(antal,1)*Masse_p/Areal/(st_nedre-st_ovre);
step=step+1;
end

 graf=xlsread('rodamin.xls');
% plot(graf(:,1),graf(:,2),'r')
% legend('Målte koncentrationer')

hold on
plot(graf(:,1),graf(:,2),'r',graf(:,1),c)
legend('Målte koncentrationer','Partikelmodel')
Afig1=sum(graf(:,2))*dt    %målte
Afig2=sum(c)*dt   %beregnede