osse/EnKF/field.m

nx=10;ny=5;nz=5;
f=load('fields.dat');
size(f)
% in to out, bottom to top
for ijim=1:length(f(:,1))
 it=0;
 for i=1:nz
  for j=1:ny
   for k=1:nx
    it=it+1;   
    state(k,j,i)=f(ijim,it);
   end
  end
 end 
 toplot(:,:)=state(1,:,:);
 figure(1);contour(toplot');drawnow
end
1	afe	1.1	nx=10;ny=5;nz=5;
2			f=load('fields.dat');
3			size(f)
4			% in to out, bottom to top
5			for ijim=1:length(f(:,1))
6			it=0;
7			for i=1:nz
8			for j=1:ny
9			for k=1:nx
10			it=it+1;
11			state(k,j,i)=f(ijim,it);
12			end
13			end
14			end
15			toplot(:,:)=state(1,:,:);
16			figure(1);contour(toplot');drawnow
17			end