/[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.4 by afe, Sat Jun 12 22:15:40 2004 UTC revision 1.5 by afe, Wed Jun 16 15:29:50 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=1;  disp=0;
5  disprms=0;  disprms=0;
6    
7    
8  %filename=sprintf('%s','../00/assimilate/pickup.in')  %filename=sprintf('%s','../da/00/assimilate/pickup.in')
9  filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out')  filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out')
10  fid=fopen(filename,'r','ieee-be');  fid=fopen(filename,'r','ieee-be');
11  foo=fread(fid,13*n+ny*nx,'float64');  foo=fread(fid,13*n+ny*nx,'float64');
# Line 19  i=6;truw=reshape(foo(i*n+1:(i+1)*n),[nx Line 19  i=6;truw=reshape(foo(i*n+1:(i+1)*n),[nx
19  for k=1:nens  for k=1:nens
20  %k=0;  %k=0;
21  filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')  filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')
22  %filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')  %filename=sprintf('%s%02i%s','../da/',k,'/assimilate/pickup.in')
23  fid=fopen(filename,'r','ieee-be');  fid=fopen(filename,'r','ieee-be');
24  foo=fread(fid,13*n+ny*nx,'float64');  foo=fread(fid,13*n+ny*nx,'float64');
25  fclose(fid);  fclose(fid);
26  filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup_nh.',iter,'.out')  %filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup_nh.',iter,'.out')
27  %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));
28  fid=fopen(filename,'r','ieee-be');  %fid=fopen(filename,'r','ieee-be');
29  foonh=fread(fid,2*n,'float64');  %foonh=fread(fid,2*n,'float64');
30  fclose(fid);  %fclose(fid);
31  i=0;memu(:,k)=foo(i*n+1:(i+1)*n);  i=0;memu(:,k)=foo(i*n+1:(i+1)*n);
32  i=3;memv(:,k)=foo(i*n+1:(i+1)*n);  i=3;memv(:,k)=foo(i*n+1:(i+1)*n);
33  i=6;memw(:,k)=foo(i*n+1:(i+1)*n);  i=6;memw(:,k)=foo(i*n+1:(i+1)*n);
# Line 46  rmsw=reshape(std(memw,0,2),[nx ny nz]); Line 46  rmsw=reshape(std(memw,0,2),[nx ny nz]);
46  rmsp=reshape(std(memp,0,2),[nx ny nz]);  rmsp=reshape(std(memp,0,2),[nx ny nz]);
47    
48  if(disprms)  if(disprms)
49  figure(1);imagesc(rmsu(:,:,z));colorbar;title('U ensemble rms');  iterstr=sprintf('%i',iter);
50  figure(2);imagesc(rmsv(:,:,z));colorbar;title('V ensemble rms');  figure(1);imagesc(rmsu(:,:,z));colorbar;
51  figure(3);imagesc(rmsw(:,:,z));colorbar;title('W ensemble rms');  title(strcat('U ensemble rms, iter ',iterstr));
52  figure(4);imagesc(rmsp(:,:,z));colorbar;title('Phi ensemble rms');  figure(2);imagesc(rmsv(:,:,z));colorbar;
53    title(strcat('V ensemble rms, iter ',iterstr));
54    figure(3);imagesc(rmsw(:,:,z));colorbar;
55    title(strcat('W ensemble rms, iter ',iterstr));
56    figure(4);imagesc(rmsp(:,:,z));colorbar;
57    title(strcat('Phi ensemble rms, iter ',iterstr));
58    end
59    
60    if (jgp)
61      for i=1:3
62        set(i,'PaperPositionMode', 'manual');
63        set(i,'PaperPosition',[0.25 0.25 8 10.5]);
64        t=get(get(get(i,'CurrentAxes'),'Title'),'String');
65        outfilename=sprintf('%s%s%i','rms',t(1),iter);
66        print('-djpeg', outfilename);
67      end
68  end  end
69    
70  if (disp)  if (disp)
# Line 60  if (disp) Line 75  if (disp)
75    subplot(3,2,2);imagesc(ensu(:,:,z));colorbar;    subplot(3,2,2);imagesc(ensu(:,:,z));colorbar;
76    subplot(3,2,4);imagesc(ensv(:,:,z));colorbar;    subplot(3,2,4);imagesc(ensv(:,:,z));colorbar;
77    subplot(3,2,6);imagesc(ensw(:,:,z));colorbar;    subplot(3,2,6);imagesc(ensw(:,:,z));colorbar;
78    title(sprintf('%s%i','truth and ensemble mean after assim',iter));    title(sprintf('%s%i','truth and ensemble mean after assim ',iter));
79  end  end

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22