% a matlab plotting program to plot float tracks 
% this version used for xyzt model, among others. 

clear all

      
load what.dat
nx = what(1)
ny = what(2)
nf = what(4)
dx = what(5);
dy = what(6);
x1 = what(7);
y1 = what(9); 


x = zeros(nx,1); y = x; 
for i=1:nx;
x(i) = (x1 + (i-1)*dx)/1000.;

end;
for i=1:ny;
y(i) = (y1 + (i-1)*dy)/1000.;
end;

       
  
 xlabel('East, km');  ylabel('North, km'); 
 
 x1 = x(1)-20
 x2 = x(nx)+20
 y1 = y(1)-20
 y2 = y(ny)+20

 load bfloat.dat 
% load floatdev.dat
% bfloat = floatdev;
 
 [m n] = size(bfloat)
 
 x = zeros(m,1); y = x;
 xc = x';
 yc = y';
 
 for i=1:nf
 
 xc(:) = bfloat(:,i);
 yc(:) = bfloat(:,i+nf);
 
 plot(xc,yc)
 plot(xc,yc,'.')
  
 axis ([x1 x2 y1 y2])
 grid on
 
 % hold on
 
                                   
 pause
 
 end