
  %I used the import wizard and got an excellent table of the data from the 
  % october_svd.reloc file!
  save('october_svd.mat', 'october_svd_data');
  % to get it back!
  load('october_svd.mat');
  % each event # is a row
  eventno = squeeze(october_svd_data(:,1));
  deptharr = squeeze(october_svd_data(:,4));
  montharr = squeeze(october_svd_data(:,12));
  dayarr = squeeze(october_svd_data(:,13));
  hourarr = squeeze(october_svd_data(:,14));
  minarr = squeeze(october_svd_data(:,15));
  secarr = squeeze(october_svd_data(:,16));
  deptherrarr = squeeze(october_svd_data(:,10));
  latarr= squeeze(october_svd_data(:,2));
  longarr= squeeze(october_svd_data(:,3));
  
  ysecarr=getysec(october_svd_data);  % times in secs from j
  booe_stn = [-40.52308 12.67527];  % station position
  distarr=getdist(october_svd_data,booe_stn); % distances in km from stn
  edroparr = ysecarr*0.0;
  cdroparr = ysecarr*0.0;
  e0 = ysecarr*0.0;
  c0 = ysecarr*0.0;
  
sacfilename='data/2005.295.00.36.25.0000.YZ.BOOE..BHZ.D.SAC';
[tracebooe,fs,sy,jsecbooe]=loadsac(sacfilename);
framehz=0.5;
[ebooe]=dobank(tracebooe,fs,framehz); 

fitcolorloop;

sacfilename='data/2005.296.00.36.25.0000.YZ.BOOE..BHZ.D.SAC';
[tracebooe,fs,sy,jsecbooe]=loadsac(sacfilename);
[ebooe]=dobank(tracebooe,fs,framehz); 

fitcolorloop;

sacfilename='data/2005.297.00.36.25.0000.YZ.BOOE..BHZ.D.SAC';
[tracebooe,fs,sy,jsecbooe]=loadsac(sacfilename);
[ebooe]=dobank(tracebooe,fs,framehz); 

fitcolorloop;




