--- MITgcm_contrib/osse/utils/ensavg.m	2004/06/09 14:20:09	1.1
+++ MITgcm_contrib/osse/utils/ensavg.m	2004/06/10 19:25:48	1.3
@@ -1,8 +1,11 @@
-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;
+z=15;
+disp=0;
 
 
-filename=sprintf('%s','../00/assimilate/pickup.in')
-%filename=sprintf('%s%03i%s','../00/assimilate/pickup.',iter,'.out')
+%filename=sprintf('%s','../00/assimilate/pickup.in')
+filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out')
 fid=fopen(filename,'r','ieee-be');
 foo=fread(fid,13*n+ny*nx,'float64');
 %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
@@ -14,8 +17,8 @@
 
 for k=1:nens
 %k=0;
-%filename=sprintf('%s%02i%s%03i%s','../',k,'/assimilate/pickup.',iter,'.out')
-filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
+filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')
+%filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
 fid=fopen(filename,'r','ieee-be');
 foo=fread(fid,13*n+ny*nx,'float64');
 %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
@@ -24,6 +27,28 @@
 i=6;memw(:,k)=foo(i*n+1:(i+1)*n);
 fclose(fid);
 end
+
 ensu(:,:,:)=reshape(mean(memu,2),[nx ny nz]);
 ensv(:,:,:)=reshape(mean(memv,2),[nx ny nz]);
 ensw(:,:,:)=reshape(mean(memw,2),[nx ny nz]);
+
+
+rmsu=reshape(std(memu,0,2),[nx ny nz]);
+rmsv=reshape(std(memv,0,2),[nx ny nz]);
+rmsw=reshape(std(memw,0,2),[nx ny nz]);
+
+%figure(1);imagesc(rmsu(:,:,z));colorbar
+%figure(2);imagesc(rmsv(:,:,z));colorbar
+%figure(3);imagesc(rmsw(:,:,z));colorbar
+
+
+if (disp)
+  subplot(3,2,1);imagesc(truu(:,:,z));colorbar;
+  subplot(3,2,3);imagesc(truv(:,:,z));colorbar;
+  subplot(3,2,5);imagesc(truw(:,:,z));colorbar;
+
+  subplot(3,2,2);imagesc(ensu(:,:,z));colorbar;
+  subplot(3,2,4);imagesc(ensv(:,:,z));colorbar;
+  subplot(3,2,6);imagesc(ensw(:,:,z));colorbar;
+  title(sprintf('%s%i','truth and ensemble mean after assim',iter));
+end