Code Matlab

 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%       Calcul des flux dans le bassin de lagunage      %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
close all
clear all
 
load Data_Kaedi.mat
 
% Données
A=1.7E-3;  % surface de la lagune en km2
uz=1; % vitesse du vent en m/s
In=40; % Flux entrant en m3/j
Inf=1; % Flux d'infiltration mm/j
Volinit=In*60;
Soleil=[8,9,9,10,10,9,9,9,9,9,8,8];
 
for i=1:1:12
    Tmax(i)=Data_Kaedi(i,1);   % Température maximale en °C
    Tmin(i)=Data_Kaedi(i,2);   % Température minimale en °C
    P(i)=Data_Kaedi(i,3);      % Précipitation en mm
    RH(i)=Data_Kaedi(i,4);     % Humidité relative en %
end
 
% Calcul de la pression de vapeur saturante de l'air en z=O e0
 
for i=1:1:12
    e0max(i)=exp(13.7)*exp(-5120/(Tmax(i)+273.3))*1013*100;
    e0min(i)=exp(13.7)*exp(-5120/(Tmin(i)+273.3))*1013*100;
end
 
% Calcul de la pression de vapeur saturante de l'eau es
 
for i=1:1:12
    esmax(i)=610.8*exp(17.17*Tmax(i)/(273.3+Tmax(i)));
    esmin(i)=610.8*exp(17.17*Tmin(i)/(273.3+Tmin(i)));
end
 
% Calcul (EVP constitue l'évaporation)
 
for i=1:1:12
    EVPmax(i)=1.69E-11*A^(-0.05)*(esmax(i)-RH(i)*e0max(i)/100)*86400000; % en mm/J (2592)
    EVPmin(i)=1.69E-11*A^(-0.05)*(esmin(i)-RH(i)*e0min(i)/100)*86400000; % en mm/J (2592)
    %multiplié par le nombre de secondes ensoleillées dans une journée !!
  
end
 
figure
plot((1:1:12),EVPmax(:),'r')
hold on
plot((1:1:12),EVPmin(:),'b')
title('Evaporation mensuelle à Kaedi')
xlabel('Mois')
ylabel('Evaporation en mm')
legend('Evaporation maximale','Evaporation minimale')
%plot((1:1:12),P(:),'c')
 
%Calcul de flux pour un jour d'un mois spécifique
 
for i=1:12
    delta(i)=In+P(i)*A*1000/30-EVPmax(i)*A*1000*Soleil(i)/24-Inf*A*1000;
end
 
figure
plot((1:1:12),delta(:),'g')
title('Différence de flux')
xlabel('Mois')
ylabel('Delta (m3)')
 
Eff=sum(delta(:))/12 %Valeur du débit constant d'effluent journalier nécessaire pour retrouver un même volume en 1 an 
 
%Calcul du volume final sur 1 an
 
Vol1=Volinit;
for i=1:12
Vol(i)=Vol1+30*(delta(i)-Eff);
Vol1=Vol(i);
end
 
figure
plot((1:1:12),Vol(:),'m')
title('Evoluvion du volume')
xlabel('Mois')
ylabel('Volume (m3)')
set(gca,'YLim',[0 3000],'Fontsize',12);
 
figure
plot((1:1:12),Vol(:)/(A*10^4),'b')
title('Evoluvion de la hauteur d eau dans la lagune')
xlabel('Mois')
ylabel('Hauteur d eau (cm)')
set(gca,'YLim',[0 150],'Fontsize',12);
 
%Scénario de remplissage commencant début mars
 
Vol2=31*delta(3)+30*delta(4); %mois de mars et avril
j=0;
while Vol2 < Volinit
     Vol2=Vol2+delta(5);
     j=j+1; %nombre de jours en mai
end
 
if j > 31
    Vol2=31*delta(3)+30*delta(4)+31*delta(5); %mois de mars, avril et mai
    k=0;
    while Vol2 < Volinit
     Vol2=Vol2+delta(6);
     k=k+1; %nombre de jours en juin
    end
    tr=31+30+31+k;
else
    tr=31+30+j;
end