osse/utils/ensavg.m

nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;iter=68;


filename=sprintf('%s','../00/assimilate/pickup.in')
%filename=sprintf('%s%03i%s','../00/assimilate/pickup.',iter,'.out')
fid=fopen(filename,'r','ieee-be');
foo=fread(fid,13*n+ny*nx,'float64');
%s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
i=0;truu=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
i=3;truv=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
i=6;truw=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);


for k=1:nens
%k=0;
%filename=sprintf('%s%02i%s%03i%s','../',k,'/assimilate/pickup.',iter,'.out')
filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
fid=fopen(filename,'r','ieee-be');
foo=fread(fid,13*n+ny*nx,'float64');
%s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
i=0;memu(:,k)=foo(i*n+1:(i+1)*n);
i=3;memv(:,k)=foo(i*n+1:(i+1)*n);
i=6;memw(:,k)=foo(i*n+1:(i+1)*n);
fclose(fid);
end
ensu(:,:,:)=reshape(mean(memu,2),[nx ny nz]);
ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);
ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);
1	afe	1.1	nx=120;ny=31;nz=29;n=nxnynz;nens=30;iter=68;
2
3
4			filename=sprintf('%s','../00/assimilate/pickup.in')
5			%filename=sprintf('%s%03i%s','../00/assimilate/pickup.',iter,'.out')
6			fid=fopen(filename,'r','ieee-be');
7			foo=fread(fid,13n+nynx,'float64');
8			%s=0; mn(k+1)=mean(foo(sn+1:sn+n));
9			i=0;truu=reshape(foo(in+1:(i+1)n),[nx ny nz]);
10			i=3;truv=reshape(foo(in+1:(i+1)n),[nx ny nz]);
11			i=6;truw=reshape(foo(in+1:(i+1)n),[nx ny nz]);
12
13
14
15			for k=1:nens
16			%k=0;
17			%filename=sprintf('%s%02i%s%03i%s','../',k,'/assimilate/pickup.',iter,'.out')
18			filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
19			fid=fopen(filename,'r','ieee-be');
20			foo=fread(fid,13n+nynx,'float64');
21			%s=0; mn(k+1)=mean(foo(sn+1:sn+n));
22			i=0;memu(:,k)=foo(in+1:(i+1)n);
23			i=3;memv(:,k)=foo(in+1:(i+1)n);
24			i=6;memw(:,k)=foo(in+1:(i+1)n);
25			fclose(fid);
26			end
27			ensu(:,:,:)=reshape(mean(memu,2),[nx ny nz]);
28			ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);
29			ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);