PSN-L Email List Message

Subject: IRIS Wilber II
From: Bobhelenmcclure@.......
Date: Sat, 10 Feb 2007 20:11:15 EST


A memo to Larry Conklin and to others who might be interested.
 
Hi Larry,
 
  You mentioned that you compare your data with PAL. I do also since  they 
are only 21km from me. Have you tried looking at BINY? They are closer to  you.
 
  WinQuake and Amaseis will open SAC binary files. You can download  such 
files from major events from almost any network station in the world easily  from 
IRIS at
 
  _http://www.iris.washington.edu/_ (http://www.iris.washington.edu/) 
 
  Click on "Earthquakes" and "Wilber II" to get to the download  page. BINY 
belongs to the US network. Do not try to select an event using  the picture 
display. Go to the listings of events instead. After you have  selected your 
event and stations (you have to unselect almost everything), be  sure you select 
"SAC BINARY individual files" as the file format. Select BHZ,  BHN, and BHE as 
the channels. Give a user name and email address. Elect to have  an email sent 
when processing is complete. If you do not ask for many files, you  will have 
a high priority and will usually get your files prepared in a few  minutes 
after the request is entered.
 
  I occasionally download events from PAL. My correspondence with  their data 
is remarkably close, once I have processed my data using my WQFilter  
utility. Something is not right with the IIR filters in WinQuake. They have more  
ringing and sharper cutoffs than they should, so I avoid using them. The FFT  
filters are OK. Cochrane did not write the filter code himself. He lifted the  
code from a book  "Analog and Digital Filter Design using C" by Les  Thede, 
which is no longer in print.
 
  The filters in Amaseis are OK. Both Dr. Jones and I worked out  the 
algorithms for Butterworth filters and ended up with identical results. I  gave him 
my algorithm for the period extending filter, and his filter  performance 
agrees with mine. I do have an optional back-and-forth filter  pass capability, 
which yields zero distortion filtering, equivalent to FFT  filtering. Amaseis 
does not implement that.
 
  If you should want to convert PSN files to SAC binary format,  WinQuake can 
do that. However, there is a bug in WinQuake that prevents transfer  of event 
time. I am able to patch in a correction on the files I convert. I told  
Cochrane what was wrong, but he has essentially retired  from maintaining his 
software, after a lifetime of good and dedicated  work.
 
  I am almost done with a heliplot program that filters and displays  WinDaq 
log files, including the one currently open, up to 27 hours before the  
present hour or any lesser number of hours.
 
Regards,
 
Bob





A memo to Larry Conklin and to others who might be interested.
 
Hi Larry,
 
  You mentioned that you compare your data with PAL. I do also sin= ce=20 they are only 21km from me. Have you tried looking at BINY? They are closer=20= to=20 you.
 
  WinQuake and Amaseis will open SAC binary files. You can downloa= d=20 such files from major events from almost any network station in the world ea= sily=20 from IRIS at
 
  http://www.iris.washington.edu/=
 
  Click on "Earthquakes" and "Wilber II" to get to the download=20 page. BINY belongs to the US network. Do not try to select an event usi= ng=20 the picture display. Go to the listings of events instead. After you have=20 selected your event and stations (you have to unselect almost everything), b= e=20 sure you select "SAC BINARY individual files" as the file format. Select BHZ= ,=20 BHN, and BHE as the channels. Give a user name and email address. Elect to h= ave=20 an email sent when processing is complete. If you do not ask for many files,= you=20 will have a high priority and will usually get your files prepared in a few=20 minutes after the request is entered.
 
  I occasionally download events from PAL. My correspondence with=20 their data is remarkably close, once I have processed my data using my WQFil= ter=20 utility. Something is not right with the IIR filters in WinQuake. They have=20= more=20 ringing and sharper cutoffs than they should, so I avoid using them. The FFT= =20 filters are OK. Cochrane did not write the filter code himself. He lifted th= e=20 code from a book  "Analog and Digital Filter Design using C" by Le= s=20 Thede, which is no longer in print.
 
  The filters in Amaseis are OK. Both Dr. Jones and I worked=20= out=20 the algorithms for Butterworth filters and ended up with identical results.=20= I=20 gave him my algorithm for the period extending filter, and his filter=20 performance agrees with mine. I do have an optional back-and-forth filt= er=20 pass capability, which yields zero distortion filtering, equivalent to FFT=20 filtering. Amaseis does not implement that.
 
  If you should want to convert PSN files to SAC binary forma= t,=20 WinQuake can do that. However, there is a bug in WinQuake that prevents tran= sfer=20 of event time. I am able to patch in a correction on the files I convert. I=20= told=20 Cochrane what was wrong, but he has essentially retired=20 from maintaining his software, after a lifetime of good and dedicated=20 work.
 
  I am almost done with a heliplot program that filters and displa= ys=20 WinDaq log files, including the one currently open, up to 27 hours before th= e=20 present hour or any lesser number of hours.
 
Regards,
 
Bob

[ Top ] [ Back ] [ Home Page ]