/[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.2 by afe, Thu Jun 10 13:45:30 2004 UTC revision 1.3 by afe, Thu Jun 10 19:25:48 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;
2    %iter=68;
3  z=15;  z=15;
4    disp=0;
5    
6    
7  %filename=sprintf('%s','../00/assimilate/pickup.in')  %filename=sprintf('%s','../00/assimilate/pickup.in')
# Line 30  ensu(:,:,:)=reshape(mean(memu,2),[nx ny Line 32  ensu(:,:,:)=reshape(mean(memu,2),[nx ny
32  ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);  ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);
33  ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);  ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);
34    
 subplot(3,2,1);imagesc(truu(:,:,z));  
 subplot(3,2,3);imagesc(truv(:,:,z));  
 subplot(3,2,5);imagesc(truw(:,:,z));  
   
 subplot(3,2,2);imagesc(ensu(:,:,z));  
 subplot(3,2,4);imagesc(ensv(:,:,z));  
 subplot(3,2,6);imagesc(ensw(:,:,z));  
 title(sprintf('%s%i','truth and ensemble mean after assim',iter));  
35    
36    rmsu=reshape(std(memu,0,2),[nx ny nz]);
37    rmsv=reshape(std(memv,0,2),[nx ny nz]);
38    rmsw=reshape(std(memw,0,2),[nx ny nz]);
39    
40    %figure(1);imagesc(rmsu(:,:,z));colorbar
41    %figure(2);imagesc(rmsv(:,:,z));colorbar
42    %figure(3);imagesc(rmsw(:,:,z));colorbar
43    
44    
45    if (disp)
46      subplot(3,2,1);imagesc(truu(:,:,z));colorbar;
47      subplot(3,2,3);imagesc(truv(:,:,z));colorbar;
48      subplot(3,2,5);imagesc(truw(:,:,z));colorbar;
49    
50      subplot(3,2,2);imagesc(ensu(:,:,z));colorbar;
51      subplot(3,2,4);imagesc(ensv(:,:,z));colorbar;
52      subplot(3,2,6);imagesc(ensw(:,:,z));colorbar;
53      title(sprintf('%s%i','truth and ensemble mean after assim',iter));
54    end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22