/[MITgcm]/MITgcm_contrib/afe/osse_MkII/utils/ensavg.m
ViewVC logotype

Annotation of /MITgcm_contrib/afe/osse_MkII/utils/ensavg.m

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


Revision 1.2 - (hide annotations) (download)
Fri Mar 11 19:29:00 2005 UTC (20 years, 4 months ago) by afe
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +11 -9 lines
added ability to deal with differing truth, pomo style

1 afe 1.1 nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;
2 afe 1.2 %xdata=nx;ydata=ny;zdata=nz;ndata=xdata*ydata*zdata
3     xdata=90;ydata=21;zdata=29;ndata=xdata*ydata*zdata
4 afe 1.1 %iter=68;
5     z=15;
6     disp=1;
7     disprms=0;
8     jpg=0;
9     vars=['U' 'V' 'W' 'P' 'T'];
10     U=find(vars=='U');
11     V=find(vars=='V');
12     W=find(vars=='W');
13     P=find(vars=='P');
14     T=find(vars=='T');
15     varnum= size(vars,2);
16    
17 afe 1.2 if(1)
18 afe 1.1 %filename=sprintf('%s','../da/00/assimilate/pickup.in')
19 afe 1.2 %filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out');
20     filename=sprintf('%s','../da/00/assimilate/pickup.0000016820.001.001.data')
21 afe 1.1 fid=fopen(filename,'r','ieee-be');
22 afe 1.2 foo=fread(fid,13*ndata+ydata*xdata,'float64');
23     i=0;tru(:,:,:,U)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
24     i=3;tru(:,:,:,V)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
25     i=6;tru(:,:,:,W)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
26     i=7;tru(:,:,:,T)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
27 afe 1.1 fclose(fid);
28 afe 1.2 end
29 afe 1.1
30     for k=1:nens
31     %k=0;
32     %filename=sprintf('%s%02i%s','../da/',k,'/assimilate/pickup.in')
33     filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')
34     %filename=sprintf('%s%02i%s','../da/i90/',k,'/assimilate/pickup.0000016820.001.001.data')
35     fid=fopen(filename,'r','ieee-be');
36     foo=fread(fid,13*n+ny*nx,'float64');
37     fclose(fid);
38     filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup_nh.',iter,'.out')
39     fid=fopen(filename,'r','ieee-be');
40     foonh=fread(fid,2*n,'float64');
41     fclose(fid);
42     i=0;mem(:,U,k)=foo(i*n+1:(i+1)*n);
43     i=3;mem(:,V,k)=foo(i*n+1:(i+1)*n);
44     i=6;mem(:,W,k)=foo(i*n+1:(i+1)*n);
45     i=7;mem(:,T,k)=foo(i*n+1:(i+1)*n);
46     i=0;mem(:,P,k)=foonh(i*n+1:(i+1)*n);
47     end
48    
49     ens=reshape(mean(mem,3),[nx ny nz varnum]);
50    
51     rms=reshape(std(mem,0,3),[nx ny nz varnum]);
52    
53     iterstr=sprintf('%i',iter);
54     if(disprms)
55     for i=1:varnum
56     figure(i);imagesc(rms(:,:,z,i));colorbar;
57     title(strcat(vars(i),' ensemble rms, iter ',iterstr));
58     end
59     end
60    
61    
62     if (jpg)
63     for i=1:varnum
64     set(i,'PaperPositionMode', 'manual');
65     set(i,'PaperPosition',[0.25 0.25 8 10.5]);
66     t=get(get(get(i,'CurrentAxes'),'Title'),'String');
67     outfilename=sprintf('%s%s%i','rms',t(1),iter);
68     print('-djpeg', outfilename);
69     end
70     end
71    
72     if (disp)
73     subplot(4,2,1);imagesc(tru(:,:,z,U));colorbar;
74     subplot(4,2,3);imagesc(tru(:,:,z,V));colorbar;
75     subplot(4,2,5);imagesc(tru(:,:,z,W));colorbar;
76     subplot(4,2,7);imagesc(tru(:,:,z,T));colorbar;
77    
78     subplot(4,2,2);imagesc(ens(:,:,z,U));colorbar;
79     subplot(4,2,4);imagesc(ens(:,:,z,V));colorbar;
80     subplot(4,2,6);imagesc(ens(:,:,z,W));colorbar;
81     subplot(4,2,8);imagesc(ens(:,:,z,T));colorbar;
82     title(sprintf('%s%i','truth and ensemble mean after assim ',iter));
83     end

  ViewVC Help
Powered by ViewVC 1.1.22