first commit

scripts_matlab/adimentionnal_number.m

+
+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
+
+
+
+
+