1 |
prefix='dynStD'; |
prefix='dynStD'; |
2 |
|
%prefix='iceStD'; |
3 |
namA={'g02','b02'}; |
namA={'g02','b02'}; |
4 |
%namA={'r17a'}; |
%namA={'r17a'}; |
5 |
Nexp=size(namA,2); |
Nexp=size(namA,2); |
33 |
[ntA(n),rList,tim,vv1,listV] = ... |
[ntA(n),rList,tim,vv1,listV] = ... |
34 |
read_StD(prefix,namA(n),listV); |
read_StD(prefix,namA(n),listV); |
35 |
nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3); |
nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3); |
|
vv1(find(vv1==undef))=NaN; |
|
36 |
%- set global dims: & load vvA --> vvB |
%- set global dims: & load vvA --> vvB |
37 |
nbV=size(listV,2); |
nbV=size(listV,2); |
38 |
nrec=nIt; n3d=nk; nReg=nRg; |
nrec=nIt; n3d=nk; nReg=nRg; |
43 |
[ntA(n),rList,tim,vv1,listV] = ... |
[ntA(n),rList,tim,vv1,listV] = ... |
44 |
read_StD(prefix,namA(n),listV); |
read_StD(prefix,namA(n),listV); |
45 |
nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3); |
nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3); |
|
vv1(find(vv1==undef))=NaN; |
|
46 |
if (nrec < nIt), |
if (nrec < nIt), |
47 |
fprintf('\n'); |
fprintf('\n'); |
48 |
error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]); |
error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]); |
129 |
figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);clf; |
figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);clf; |
130 |
var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp); |
var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp); |
131 |
for n=1:Nexp, |
for n=1:Nexp, |
132 |
dd(:,n)=max(var(isA(n):ieA(n),:,n))-min(var(isA(n):ieA(n),:,n)); |
tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0; |
133 |
av(:,n)=mean(var(isA(n):ieA(n),:,n)); |
av(:,n)=sum(tmp); |
134 |
|
if length(I) < size(tmp,1), av(:,n)=av(:,n)/(size(tmp,1)-length(I)); end |
135 |
|
tmp(I,:)=NaN; dd(:,n)=max(tmp)-min(tmp); |
136 |
end |
end |
137 |
|
var(find(var==undef))=NaN; |
138 |
for nv=1:4, |
for nv=1:4, |
139 |
subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:'; |
subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:'; |
140 |
for n=1:Nexp, |
for n=1:Nexp, |