/[MITgcm]/MITgcm_contrib/osse/utils/ensavg.m
ViewVC logotype

Contents of /MITgcm_contrib/osse/utils/ensavg.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Wed Jun 9 14:20:09 2004 UTC (21 years, 1 month ago) by afe
Branch: MAIN
scripts for ensemble diagnostics

1 nx=120;ny=31;nz=29;n=nx*ny*nz;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,13*n+ny*nx,'float64');
8 %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
9 i=0;truu=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
10 i=3;truv=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
11 i=6;truw=reshape(foo(i*n+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,13*n+ny*nx,'float64');
21 %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
22 i=0;memu(:,k)=foo(i*n+1:(i+1)*n);
23 i=3;memv(:,k)=foo(i*n+1:(i+1)*n);
24 i=6;memw(:,k)=foo(i*n+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]);

  ViewVC Help
Powered by ViewVC 1.1.22