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	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