--- MITgcm_contrib/osse/utils/ensavg.m	2004/06/12 22:15:40	1.4
+++ MITgcm_contrib/osse/utils/ensavg.m	2004/06/16 15:29:50	1.5
@@ -1,11 +1,11 @@
 nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;
 %iter=68;
 z=15;
-disp=1;
+disp=0;
 disprms=0;
 
 
-%filename=sprintf('%s','../00/assimilate/pickup.in')
+%filename=sprintf('%s','../da/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');
@@ -19,15 +19,15 @@
 for k=1:nens
 %k=0;
 filename=sprintf('%s%02i%s%03i%s','../da/',k,'/assimilate/pickup.',iter,'.out')
-%filename=sprintf('%s%02i%s','../',k,'/assimilate/pickup.in')
+%filename=sprintf('%s%02i%s','../da/',k,'/assimilate/pickup.in')
 fid=fopen(filename,'r','ieee-be');
 foo=fread(fid,13*n+ny*nx,'float64');
 fclose(fid);
-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')
 %s=0; mn(k+1)=mean(foo(s*n+1:s*n+n));
-fid=fopen(filename,'r','ieee-be');
-foonh=fread(fid,2*n,'float64');
-fclose(fid);
+%fid=fopen(filename,'r','ieee-be');
+%foonh=fread(fid,2*n,'float64');
+%fclose(fid);
 i=0;memu(:,k)=foo(i*n+1:(i+1)*n);
 i=3;memv(:,k)=foo(i*n+1:(i+1)*n);
 i=6;memw(:,k)=foo(i*n+1:(i+1)*n);
@@ -46,10 +46,25 @@
 rmsp=reshape(std(memp,0,2),[nx ny nz]);
 
 if(disprms)
-figure(1);imagesc(rmsu(:,:,z));colorbar;title('U ensemble rms');
-figure(2);imagesc(rmsv(:,:,z));colorbar;title('V ensemble rms');
-figure(3);imagesc(rmsw(:,:,z));colorbar;title('W ensemble rms');
-figure(4);imagesc(rmsp(:,:,z));colorbar;title('Phi ensemble rms');
+iterstr=sprintf('%i',iter);
+figure(1);imagesc(rmsu(:,:,z));colorbar;
+title(strcat('U ensemble rms, iter ',iterstr));
+figure(2);imagesc(rmsv(:,:,z));colorbar;
+title(strcat('V ensemble rms, iter ',iterstr));
+figure(3);imagesc(rmsw(:,:,z));colorbar;
+title(strcat('W ensemble rms, iter ',iterstr));
+figure(4);imagesc(rmsp(:,:,z));colorbar;
+title(strcat('Phi ensemble rms, iter ',iterstr));
+end
+
+if (jgp)
+  for i=1:3
+    set(i,'PaperPositionMode', 'manual');
+    set(i,'PaperPosition',[0.25 0.25 8 10.5]);
+    t=get(get(get(i,'CurrentAxes'),'Title'),'String');
+    outfilename=sprintf('%s%s%i','rms',t(1),iter);
+    print('-djpeg', outfilename);
+  end
 end
 
 if (disp)
@@ -60,5 +75,5 @@
   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));
+  title(sprintf('%s%i','truth and ensemble mean after assim ',iter));
 end