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

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

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

revision 1.1 by afe, Wed Jun 9 14:20:09 2004 UTC revision 1.2 by afe, Thu Jun 10 13:45:30 2004 UTC
# Line 1  Line 1 
1  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;iter=68;  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;iter=68;
2    z=15;
3    
4    
5  filename=sprintf('%s','../00/assimilate/pickup.in')  %filename=sprintf('%s','../00/assimilate/pickup.in')
6  %filename=sprintf('%s%03i%s','../00/assimilate/pickup.',iter,'.out')  filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out')
7  fid=fopen(filename,'r','ieee-be');  fid=fopen(filename,'r','ieee-be');
8  foo=fread(fid,13*n+ny*nx,'float64');  foo=fread(fid,13*n+ny*nx,'float64');
9  %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));  %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
# Line 14  i=6;truw=reshape(foo(i*n+1:(i+1)*n),[nx Line 15  i=6;truw=reshape(foo(i*n+1:(i+1)*n),[nx
15    
16  for k=1:nens  for k=1:nens
17  %k=0;  %k=0;
18  %filename=sprintf('%s%02i%s%03i%s','../',k,'/assimilate/pickup.',iter,'.out')  filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')
19  filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')  %filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
20  fid=fopen(filename,'r','ieee-be');  fid=fopen(filename,'r','ieee-be');
21  foo=fread(fid,13*n+ny*nx,'float64');  foo=fread(fid,13*n+ny*nx,'float64');
22  %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));  %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
# Line 24  i=3;memv(:,k)=foo(i*n+1:(i+1)*n); Line 25  i=3;memv(:,k)=foo(i*n+1:(i+1)*n);
25  i=6;memw(:,k)=foo(i*n+1:(i+1)*n);  i=6;memw(:,k)=foo(i*n+1:(i+1)*n);
26  fclose(fid);  fclose(fid);
27  end  end
28    
29  ensu(:,:,:)=reshape(mean(memu,2),[nx ny nz]);  ensu(:,:,:)=reshape(mean(memu,2),[nx ny nz]);
30  ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);  ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);
31  ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);  ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);
32    
33    subplot(3,2,1);imagesc(truu(:,:,z));
34    subplot(3,2,3);imagesc(truv(:,:,z));
35    subplot(3,2,5);imagesc(truw(:,:,z));
36    
37    subplot(3,2,2);imagesc(ensu(:,:,z));
38    subplot(3,2,4);imagesc(ensv(:,:,z));
39    subplot(3,2,6);imagesc(ensw(:,:,z));
40    title(sprintf('%s%i','truth and ensemble mean after assim',iter));
41    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22