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
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