--- MITgcm_contrib/osse/utils/ensavg.m	2004/06/16 21:03:04	1.6
+++ MITgcm_contrib/osse/utils/ensavg.m	2004/07/01 13:45:09	1.7
@@ -1,14 +1,15 @@
 nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;
 %iter=68;
 z=15;
-disp=0;
-disprms=1;
+disp=1;
+disprms=0;
 jpg=0;
-vars=['U' 'V' 'W' 'P'];
+vars=['U' 'V' 'W' 'P' 'T'];
 U=find(vars=='U');
 V=find(vars=='V');
 W=find(vars=='W');
 P=find(vars=='P');
+T=find(vars=='T');
 varnum= size(vars,2);
 
 
@@ -20,6 +21,7 @@
 i=0;tru(:,:,:,U)=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]);
 i=6;tru(:,:,:,W)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
+i=7;tru(:,:,:,T)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);
 
 
 
@@ -38,6 +40,7 @@
 i=0;mem(:,U,k)=foo(i*n+1:(i+1)*n);
 i=3;mem(:,V,k)=foo(i*n+1:(i+1)*n);
 i=6;mem(:,W,k)=foo(i*n+1:(i+1)*n);
+i=7;mem(:,T,k)=foo(i*n+1:(i+1)*n);
 i=0;mem(:,P,k)=foonh(i*n+1:(i+1)*n);
 end
 
@@ -65,12 +68,14 @@
 end
 
 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;
+  subplot(4,2,1);imagesc(tru(:,:,z,U));colorbar;
+  subplot(4,2,3);imagesc(tru(:,:,z,V));colorbar;
+  subplot(4,2,5);imagesc(tru(:,:,z,W));colorbar;
+  subplot(4,2,7);imagesc(tru(:,:,z,T));colorbar;
+
+  subplot(4,2,2);imagesc(ens(:,:,z,U));colorbar;
+  subplot(4,2,4);imagesc(ens(:,:,z,V));colorbar;
+  subplot(4,2,6);imagesc(ens(:,:,z,W));colorbar;
+  subplot(4,2,8);imagesc(ens(:,:,z,T));colorbar;
   title(sprintf('%s%i','truth and ensemble mean after assim ',iter));
 end