pysep icon indicating copy to clipboard operation
pysep copied to clipboard

Inconsistent kztime in SAC Headers affecting Moment Tensor Inversion in MTUQ

Open Aayushi-10 opened this issue 1 year ago • 1 comments

I am encountering an issue while downloading data with PySEP for the networks 8H, 7B, and 5C. The kztime in the SAC headers for all the station differs, which is causing difficulties in performing moment tensor inversion using the downloaded data.

I utilized the following file for data extraction: https://github.com/Aayushi-10/Test/blob/main/2016-12-12T220444_HINDU_KUSH.yaml

And the following file for conducting the MTUQ inversion: https://github.com/Aayushi-10/Test/blob/main/HinduKush_Sample.py

This is the error which I’m facing:

Screenshot (447)

Aayushi-10 avatar Sep 26 '24 15:09 Aayushi-10

Hi @Aayushi-10, is this related to #146?

kztime is not explicitly set in PySEP and is built from tr.stats.starttime - tr.stats.sac.b (from here: https://github.com/obspy/obspy/issues/163). PySEP sets b here (https://github.com/adjtomo/pysep/blob/master/pysep/utils/cap_sac.py#L188), so kztime should just be the event origin time.

Can you provide log files for your PySEP data download for one station? That might help diagnose the issue.

bch0w avatar Oct 01 '24 19:10 bch0w

Hi @bch0w,

Thank you for responding.

Yes, this issue is related to #146. Below is the link to the log file generated when I was using network 8H and station EP02:

https://github.com/Aayushi-10/Test/blob/main/pysep.log

Aayushi-10 avatar Oct 15 '24 07:10 Aayushi-10