adimentionnal_number.m 1.13 KB

clear all;
close all;


fig=figure;

% fig.Units='inches';
% fig.Position=[6 6 9 3.2];
% fig.OuterPosition=[6 6 9 3.2];
% 
% fig.PaperUnits = 'inches';
% fig.PaperPosition = [0 0 9 3.2];
% fig.PaperPositionMode = 'manual';
% fig.PaperOrientation='landscape';
% fig.PaperSize=[9 3.2];
% 


for n=3
titre1=['10^{-6}';'10^{-4}';'10^{-3}';'10^{-2}'];
titre2=cellstr(titre1);

dos1=['1e6      ';'1e4      ';'test     ';'1e2      '];
dos2=cellstr(dos1);


dossier=[num2str(n),'mu1';num2str(n),'mu2';num2str(n),'mu3';dos2(n);num2str(n),'mu4'];
dossier=cellstr(dossier);


diffus=['1e-2_uniforme'];

mix=[1e-6 1e-4 1e-3 1e-2];
mu=[0.5 1 1.5 2 2.5]/86400;

col=['mo';'ro';'bo';'ko'];
col=cellstr(col);

for i=1:5
dos=char(dossier(i));

number(:,n)=sqrt(mu/mix(n));


lagbiomass=load(['~/projets/correction_memo/output/',dos,'/weight_',diffus,'.dat']);
eulbiomass=load(['~/projets/correction_memo/output/',dos,'/eulerian_',diffus,'.dat']);

biolag_tot=sum(lagbiomass,2)*0.5;
bioeul_tot=sum(eulbiomass,2)*0.5;

std=sqrt((biolag_tot-bioeul_tot).^2./bioeul_tot.^2)*100;
std(1:3)=0;

maxstd(i,n)=max(std);
maxbio(i,n)=max(biolag_tot);


plot(std)
hold on

end

end