Blame view

scripts_matlab/adimentionnal_number.m 1.13 KB
68586e03   Jérémy Baudry   new release
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64

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