Mitchell Liddell <liddell<at>ualberta.ca>
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to
make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to
make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
-
Hi Mitch,
Based on the manual, you have to instruct what SAC should do when you have
gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you
can use AVERAGE. It's a decision you will have to make. In your case, using
AVERAGE might be better, although the decision should be based on the data
you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to
make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
-
Hello Januka,
Thank you very much for the response. Unfortunately, if I run commands in
that way I get an error because SAC thinks I'm trying to merge the file
"GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC GS.OK030.00.HHZ.D.2015.015.212306.SAC
GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC GS.OK030.00.HHZ.D.2015.015.213111.SAC
GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options
rather than the file list, as there is no example on the manual page (
https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how to
control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
--
Based on the manual, you have to instruct what SAC should do when you have
gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you
can use AVERAGE. It's a decision you will have to make. In your case, using
AVERAGE might be better, although the decision should be based on the data
you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
----------------------
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to
make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
-
Try
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
Thank you very much for the response. Unfortunately, if I run commands in
that way I get an error because SAC thinks I'm trying to merge the file
"GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC
GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC
GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options
rather than the file list, as there is no example on the manual page (
https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how
to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
--
Based on the manual, you have to instruct what SAC should do when you
have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you
can use AVERAGE. It's a decision you will have to make. In your case, using
AVERAGE might be better, although the decision should be based on the data
you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
----------------------
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to
make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
-
Hello Januka,
I don't think the way that I'm reading in the files is the issue. I have
tried reading in all the files at one time and then attempting merge, and
also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how to
use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com> wrote:
Try
--
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Hello Januka,
----------------------
Thank you very much for the response. Unfortunately, if I run commands in
that way I get an error because SAC thinks I'm trying to merge the file
"GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC
GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC
GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options
rather than the file list, as there is no example on the manual page (
https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how
to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
--
Based on the manual, you have to instruct what SAC should do when you
have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE
you can use AVERAGE. It's a decision you will have to make. In your case,
using AVERAGE might be better, although the decision should be based on the
data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
----------------------
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My
commands are working if there is a gap in the data (fills in with 0s) but
if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option
to make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
-
What computer platform are you using? Apparently MERGE does not work right on Macs. If you search the web you can find lots of links about MERGE problems.
Cliff Thurber
On Jan 10, 2019, at 3:01 AM, Mitchell Liddell <liddell<at>ualberta.ca<liddell<at>ualberta.ca>> wrote:
Hello Januka,
I don't think the way that I'm reading in the files is the issue. I have tried reading in all the files at one time and then attempting merge, and also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how to use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com<jattanayake<at>gmail.com>> wrote:
Try
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
[https://docs.google.com/uc?export=download&id=18UprUc6tr2T7pShyWEVxhJNmiEGEFySy&revid=0B2dyBa1wSDlRV3Z5a1Z4UEl5WHhHaExLc1M4S0s1blFLc3RrPQ]
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca<liddell<at>ualberta.ca>> wrote:
Hello Januka,
Thank you very much for the response. Unfortunately, if I run commands in that way I get an error because SAC thinks I'm trying to merge the file "GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options rather than the file list, as there is no example on the manual page (https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com<jattanayake<at>gmail.com>> wrote:
Hi Mitch,
Based on the manual, you have to instruct what SAC should do when you have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you can use AVERAGE. It's a decision you will have to make. In your case, using AVERAGE might be better, although the decision should be based on the data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
[https://docs.google.com/uc?export=download&id=18UprUc6tr2T7pShyWEVxhJNmiEGEFySy&revid=0B2dyBa1wSDlRV3Z5a1Z4UEl5WHhHaExLc1M4S0s1blFLc3RrPQ]
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca<liddell<at>ualberta.ca>> wrote:
Mitchell Liddell <liddell<at>ualberta.ca<liddell<at>ualberta.ca>>
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My commands are working if there is a gap in the data (fills in with 0s) but if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (ds.iris.edu/files/sac-manual/commands/merge.htmlhttp://ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to make data overlaps average.
I could not find the syntax for applying the options for data gaps or overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
-
Hello,
I downloaded 10 days continuous waveform data in seed format from BREQ_FAST
Request. And try to convert it into SAC format. I applied rdseed -df
(seed file name) command and got the following error. However, for
one-day data,
I had no error. In this regards, Could you please give me some suggestions?
Thank you in advance!
My error is:
Writing XQ.NA010..HH2, 19999039 samples (binary), starting 2015,229
03:46:38.6700 UT
WARNING (process_data): too much data (burp!).
Maximum data length exceeded for station NA010, channel HH2.
Due to indigestion, this seismogram will be split into sections.
Try increasing the Sample Buffer Length
Writing XQ.NA010..HH2, 19998908 samples (binary), starting 2015,230
07:33:13.8650 UT
WARNING (process_data): too much data (burp!).
Maximum data length exceeded for station NA010, channel HH2.
Due to indigestion, this seismogram will be split into sections.
Try increasing the Sample Buffer Length
Sanjev Dhakal
On Fri, Jan 11, 2019 at 1:30 AM Cliff Thurber <cthurber<at>wisc.edu> wrote:
What computer platform are you using? Apparently MERGE does not work
--
right on Macs. If you search the web you can find lots of links about
MERGE problems.
Cliff Thurber
On Jan 10, 2019, at 3:01 AM, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
I don't think the way that I'm reading in the files is the issue. I have
tried reading in all the files at one time and then attempting merge, and
also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how
to use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com> wrote:
Try
--
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Hello Januka,
----------------------
Thank you very much for the response. Unfortunately, if I run commands
in that way I get an error because SAC thinks I'm trying to merge the file
"GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC
GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC
GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options
rather than the file list, as there is no example on the manual page (
https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows
how to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
--
Based on the manual, you have to instruct what SAC should do when you
have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE
you can use AVERAGE. It's a decision you will have to make. In your case,
using AVERAGE might be better, although the decision should be based on the
data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
----------------------
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments.
My commands are working if there is a gap in the data (fills in with 0s)
but if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option
to make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Thank you!
Best Regards
*------------------------------------------------------------*
*Sanjev Dhakal*
Institute of Tibetan Plateau Research, Chinese Academy of Sciences
(ITP-CAS).
*Address:* No.16, Lincui Rd., Chaoyang District, Beijing 100101, China
*Cell:* *+86- <+86%20130%205144%203523>18811333057 (China) and
+977-9849335053 (Nepal)*
Alternative E-mail:
*Sanjev<at>itpcas.ac.cn <Sanjev<at>itpcas.ac.cn>*
-
-
Dear All -
This could be down to a difference between IRIS/SAC and BRIS/SAC. The syntax of the MERGE command changed incompatibly at some point in the evolution of IRIS’ SAC that the BRIS version does not implement.
Suggest that you handle the overlap by using CUTIM to trim the overlap from the partially merged in-memory merged trace before doing the next merge.
As Januka says, merging data without knowing why the overlap arises is worrisome.
On 10 Jan 2019, at 18:00, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
George Helffrich
I don't think the way that I'm reading in the files is the issue. I have tried reading in all the files at one time and then attempting merge, and also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how to use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com> wrote:
Try
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
Thank you very much for the response. Unfortunately, if I run commands in that way I get an error because SAC thinks I'm trying to merge the file "GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options rather than the file list, as there is no example on the manual page (https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
Based on the manual, you have to instruct what SAC should do when you have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you can use AVERAGE. It's a decision you will have to make. In your case, using AVERAGE might be better, although the decision should be based on the data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My commands are working if there is a gap in the data (fills in with 0s) but if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to make data overlaps average.
I could not find the syntax for applying the options for data gaps or overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
george<at>elsi.jp
-
Mitchell,
Did you find a solution to your problem?
As George Helffrich suggests the problem may be differences between IRIS/SAC and BRIS/SAC, or possibly an older version of IRIS/SAC. Which version are you running? The copyright shown when sac starts should identify it rather easily.
Januka’s suggestion is based on v101.6a of IRIS/SAC.
Brian
On Jan 10, 2019, at 4:00 AM, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
I don't think the way that I'm reading in the files is the issue. I have tried reading in all the files at one time and then attempting merge, and also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how to use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com <jattanayake<at>gmail.com>> wrote:
Try
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca <liddell<at>ualberta.ca>> wrote:
Hello Januka,
Thank you very much for the response. Unfortunately, if I run commands in that way I get an error because SAC thinks I'm trying to merge the file "GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options rather than the file list, as there is no example on the manual page (https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows how to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com <jattanayake<at>gmail.com>> wrote:
Hi Mitch,
Based on the manual, you have to instruct what SAC should do when you have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE you can use AVERAGE. It's a decision you will have to make. In your case, using AVERAGE might be better, although the decision should be based on the data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca <liddell<at>ualberta.ca>> wrote:
Mitchell Liddell <liddell<at>ualberta.ca <liddell<at>ualberta.ca>>
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments. My commands are working if there is a gap in the data (fills in with 0s) but if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (ds.iris.edu/files/sac-manual/commands/merge.html http://ds.iris.edu/files/sac-manual/commands/merge.html) there is an option to make data overlaps average.
I could not find the syntax for applying the options for data gaps or overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
-
Hello Brian,
I suspect the issue has to do with the IRIS/SAC vs BRIS/SAC differences as
George Hellfrich proposed. In my case, I had very situations that required
an overlap so it was just easier to trim them instead of working out a
solution to the overlap itself.
The data gap filling function of merge has been fine, although I haven't
been able to find the syntax to change the command from fill with 0s to
interpolate.
Mitch
On Mon, Jan 14, 2019 at 7:48 AM Brian Savage <savage<at>uri.edu> wrote:
Mitchell,
--
Did you find a solution to your problem?
As George Helffrich suggests the problem may be differences between
IRIS/SAC and BRIS/SAC, or possibly an older version of IRIS/SAC. Which
version are you running? The copyright shown when sac starts should
identify it rather easily.
Januka’s suggestion is based on v101.6a of IRIS/SAC.
Brian
On Jan 10, 2019, at 4:00 AM, Mitchell Liddell <liddell<at>ualberta.ca> wrote:
Hello Januka,
I don't think the way that I'm reading in the files is the issue. I have
tried reading in all the files at one time and then attempting merge, and
also merging one by one like:
r file1
merge file2
merge file3
.
.
.
write merged.sac
This works as long as there are no overlaps because I can't find out how
to use the option in the merge command.
Does that make sense?
Thanks!
Mitch
On Wed, Jan 9, 2019 at 10:54 PM Geo . <jattanayake<at>gmail.com> wrote:
Try
--
READ *GS.OK030.00.HHZ.D.2015.015*.SAC
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 15:16, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Hello Januka,
----------------------
Thank you very much for the response. Unfortunately, if I run commands
in that way I get an error because SAC thinks I'm trying to merge the file
"GAP" to "ZERO" etc:
SAC> READ GS.OK030.00.HHZ.D.2015.015*SAC
GS.OK030.00.HHZ.D.2015.015.181259.SAC
GS.OK030.00.HHZ.D.2015.015.212306.SAC GS.OK030.00.HHZ.D.2015.015.212213.SAC
GS.OK030.00.HHZ.D.2015.015.180849.SAC
GS.OK030.00.HHZ.D.2015.015.213111.SAC GS.OK030.00.HHZ.D.2015.015.000000.SAC
SAC> MERGE GAP ZERO OVERLAP AVERAGE
ERROR 108: File does not exist: GAP
I don't know the syntax to make SAC know that I am providing the options
rather than the file list, as there is no example on the manual page (
https://ds.iris.edu/files/sac-manual/commands/merge.html) that shows
how to control these. Do you have any ideas?
Best,
Mitch
On Wed, Jan 9, 2019 at 3:27 PM Geo . <jattanayake<at>gmail.com> wrote:
Hi Mitch,
--
Based on the manual, you have to instruct what SAC should do when you
have gaps or overlapping time series.
e.g.
READ file1 file2
MERGE GAP ZERO OVERLAP COMPARE
Note that in place of ZERO you can use INTERP and in place of COMPARE
you can use AVERAGE. It's a decision you will have to make. In your case,
using AVERAGE might be better, although the decision should be based on the
data you have.
Januka
Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: *https://sites.google.com/site/janukaattanayake
https://sites.google.com/site/janukaattanayake*
School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia
On Thu, 10 Jan 2019 at 04:12, Mitchell Liddell <liddell<at>ualberta.ca>
wrote:
Mitchell Liddell <liddell<at>ualberta.ca>
----------------------
10:33 PM (1 minute ago)
to sac-help
Hello,
I am attempting to merge several sac files into single day segments.
My commands are working if there is a gap in the data (fills in with 0s)
but if there is an overlap I get, for example, this error:
SAC> r GS.OK030.00.HHZ.D.2015.327.000000.SAC
SAC> merge GS.OK030.00.HHZ.D.2015.327.062124.SAC
ERROR 1802: Time overlap: GS.OK030.00.HHZ.D.2015.327.000000.SAC
GS.OK030.00.HHZ.D.2015.327.062124.SAC
END1: 2015 327 6 24 50 580
BEG2: 2015 327 6 21 24 890
On the manual page for MERGE (
ds.iris.edu/files/sac-manual/commands/merge.html) there is an option
to make data overlaps average.
I could not find the syntax for applying the options for data gaps or
overlaps in MERGE. Can you help me with this?
Thank you!
Mitch
--
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
----------------------
SAC Help
Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu
Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
Update subscription preferences at http://ds.iris.edu/account/profile/
Mitch V. Liddell, PhD
Postdoctoral Researcher
University of Alberta
587-372-5060
-
-
-
-
-