clear;
%close;

fig=figure('Units','centimeters','Position',[0.0 0.0 21.0 21.0]...
,'PaperUnits','centimeters','PaperPosition',[0.0 0.0 21.0 21.0]);

%len=40
%wid=40
%rt=10
len=130
wid=130
rt=30

load slip;
time_dat=slip;
L=time_dat(1,1)/1000;
W=time_dat(1,2)/1000;
NL=time_dat(2,1)
NW=time_dat(2,2)

ii=2;
for iy=1:NW
  for ix=1:NL
       ii=ii+1;
	  x(ix)=time_dat(ii,1)/1000;
	  y(iy)=time_dat(ii,2)/1000;
            t1(iy,ix)=time_dat(ii,3);
	  time1(ii)=time_dat(ii,3);
  end
end

tmax=max(time1)
mesh(x,y,t1);
xlabel('Length (km)')
ylabel('Width (km)')
zlabel('Slip (m)')
title('Slip (m)')
axis([0 len 0 wid 0 rt])
grid
view(35,50)
print slip.ps
