Inconsistent kztime in SAC Headers affecting Moment Tensor Inversion in MTUQ
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:
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.
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