/[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.6 by afe, Wed Jun 16 21:03:04 2004 UTC revision 1.8 by afe, Thu Jul 1 13:48:05 2004 UTC
# Line 1  Line 1 
1  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;
2  %iter=68;  %iter=68;
3  z=15;  z=15;
4  disp=0;  disp=1;
5  disprms=1;  disprms=0;
6  jpg=0;  jpg=0;
7  vars=['U' 'V' 'W' 'P'];  vars=['U' 'V' 'W' 'P' 'T'];
8  U=find(vars=='U');  U=find(vars=='U');
9  V=find(vars=='V');  V=find(vars=='V');
10  W=find(vars=='W');  W=find(vars=='W');
11  P=find(vars=='P');  P=find(vars=='P');
12    T=find(vars=='T');
13  varnum= size(vars,2);  varnum= size(vars,2);
14    
15    
# Line 20  foo=fread(fid,13*n+ny*nx,'float64'); Line 21  foo=fread(fid,13*n+ny*nx,'float64');
21  i=0;tru(:,:,:,U)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=0;tru(:,:,:,U)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
22  i=3;tru(:,:,:,V)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=3;tru(:,:,:,V)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
23  i=6;tru(:,:,:,W)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=6;tru(:,:,:,W)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
24    i=7;tru(:,:,:,T)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
25    fclose(fid);
26    
27    
28  for k=1:nens  for k=1:nens
# Line 38  fclose(fid); Line 40  fclose(fid);
40  i=0;mem(:,U,k)=foo(i*n+1:(i+1)*n);  i=0;mem(:,U,k)=foo(i*n+1:(i+1)*n);
41  i=3;mem(:,V,k)=foo(i*n+1:(i+1)*n);  i=3;mem(:,V,k)=foo(i*n+1:(i+1)*n);
42  i=6;mem(:,W,k)=foo(i*n+1:(i+1)*n);  i=6;mem(:,W,k)=foo(i*n+1:(i+1)*n);
43    i=7;mem(:,T,k)=foo(i*n+1:(i+1)*n);
44  i=0;mem(:,P,k)=foonh(i*n+1:(i+1)*n);  i=0;mem(:,P,k)=foonh(i*n+1:(i+1)*n);
45  end  end
46    
# Line 65  if (jpg) Line 68  if (jpg)
68  end  end
69    
70  if (disp)  if (disp)
71    subplot(3,2,1);imagesc(truu(:,:,z));colorbar;    subplot(4,2,1);imagesc(tru(:,:,z,U));colorbar;
72    subplot(3,2,3);imagesc(truv(:,:,z));colorbar;    subplot(4,2,3);imagesc(tru(:,:,z,V));colorbar;
73    subplot(3,2,5);imagesc(truw(:,:,z));colorbar;    subplot(4,2,5);imagesc(tru(:,:,z,W));colorbar;
74      subplot(4,2,7);imagesc(tru(:,:,z,T));colorbar;
75    subplot(3,2,2);imagesc(ensu(:,:,z));colorbar;  
76    subplot(3,2,4);imagesc(ensv(:,:,z));colorbar;    subplot(4,2,2);imagesc(ens(:,:,z,U));colorbar;
77    subplot(3,2,6);imagesc(ensw(:,:,z));colorbar;    subplot(4,2,4);imagesc(ens(:,:,z,V));colorbar;
78      subplot(4,2,6);imagesc(ens(:,:,z,W));colorbar;
79      subplot(4,2,8);imagesc(ens(:,:,z,T));colorbar;
80    title(sprintf('%s%i','truth and ensemble mean after assim ',iter));    title(sprintf('%s%i','truth and ensemble mean after assim ',iter));
81  end  end
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22