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
| load('MAT_BARILOCHE_STAGE.mat')
SAOD_BAR=MAT_JOURS(:,7)'; % VISIBLE (550 nm)
Vect=numel(SAOD_BAR(1,:));
MAT_YEAR=MAT_JOURS(:,1);
MAT_MOIS=MAT_JOURS(:,2);
u=0;
for j=2012:2017
u=u+1;
BLOC_1=find(MAT_YEAR(:,1)==j);
BLOC_ANNEE_T=SAOD_BAR(:,BLOC_1);
k=0;
for t=1:12
k=k+1;
x1=find(MAT_MOIS(BLOC_1,1)==t);
if numel(x1)==0;
MOY_J=NaN;
DEV_J=NaN;
elseif numel(x1)==1;
MOY_J=NaN;
DEV_J=NaN;
else
MOY_J1=nanmean(BLOC_ANNEE_T(:,x1));
DEV_J1=nanstd(BLOC_ANNEE_T(:,x1));
MOY_J=MOY_J1;
DEV_J=DEV_J1;
end
MOY_T(k,1)=MOY_J;
DEV_T(k,1)=DEV_J;
end
Matrice(:,u)=MOY_T;
Matrice_dev(:,u)=DEV_T;
end
r1=numel(Matrice(1,:)); % Stockage des valeurs de "Valeur" dans une matrice colonne.
o1=0;
for m1=1:r1
for m11=1:numel(Matrice(:,1));
o1=o1+1;
values11=cat(1,Matrice(m11,m1));
SAOD_MENS(o1,:)=values11;
end
end
r11=numel(Matrice_dev(1,:)); % Stockage des valeurs de "Valeur" dans une matrice colonne.
b1=0;
for m11=1:r11
for m111=1:numel(Matrice_dev(:,1));
b1=b1+1;
values111=cat(1,Matrice_dev(m111,m11));
SAOD_DISP(b1,:)=values111;
end
end |
Partager