clear 
close all

% In this script we look at records for the Sumatra earthquake from
% different sensor, look at the spectrum and try to estimate the
% peak-to-peak displacement recorded in Germany

% Load data
disp(' Loading data ... please wait! ')
load sumatra

% sampling interval
dt=1/20;

Zdata=Zdata/1e9;   % conversion from nanometers/s to meters/s
Rdata=Rdata/1e9;
Tdata=Tdata/1e9;

Rtime=Rtime-54*60+59.988998;  % conversion of the rotation time


% Plot all observations
subplot(6,1,1),plot(Ttime,Zdata),ylabel(' Z '),axis([800 7200 -Inf Inf])
title(' The Sumatra M9.0 Main Shock, 26 Dec 2004 ')
subplot(6,1,2),plot(Ttime,Rdata),ylabel(' R '),axis([800 7200 -Inf Inf])
subplot(6,1,3),plot(Ttime,Tdata),ylabel(' T '),axis([800 7200 -Inf Inf])
subplot(6,1,4),plot(Rtime,RZdata),ylabel(' Rot_Z '),axis([800 7200 -Inf Inf])
subplot(6,1,5),plot(Rtime,TNSdata),ylabel(' Tilt N-S '),axis([800 7200 -Inf Inf])
subplot(6,1,6),plot(Rtime,TEWdata),ylabel(' Tilt E-W '),axis([800 7200 -Inf Inf])
xlabel(' Time (s) ')

drawnow