%
%  
%
g=9.81;
load EXP15_WDL25.txt
t=EXP15_WDL25(:,1);
N=length(t);
tN=t(N);
i=1:N;
x=EXP15_WDL25(:,2);
y=EXP15_WDL25(:,3);
z=EXP15_WDL25(:,4);
cal=sqrt(x.*x+y.*y+z.*z)/g;
%
%   plot each channel
%
figure(1)
subplot(4,1,1);plot(t,x);
axis([680,730,-20.,20.])
axis([0,tN,-12.,12.])
ylabel('acceleration m/s**2')
title('Nov. 18, EXP15_WDL25')
subplot(4,1,2);plot(t,y);
axis([680,730,-20.,20.])
axis([0,tN,-12.,12.])
ylabel('acceleration m/s**2')
title('CH2')
subplot(4,1,3);plot(t,z);
axis([680,730,-20.,20.])
axis([0,tN,-2.,22.])
ylabel('acceleration m/s**2')
title('CH3')
subplot(4,1,4);plot(t,cal);
axis([680,730,0.0,3.0])
axis([0,tN,0.5,1.5])
ylabel('|g|/g0')
xlabel('time (s) ')
%
