petermakus / seismic Goto Github PK

Problem

Lately, in particular, scipy and numpy have the habit of changing their syntax.

Proposed solution

Fix version codes in the various dependency files.

Problem

Currently, CorrDB does not differ between data from various locations.

Proposed solution

The probably easiest way to implement this into the current framework would be to create an additional subgroup in the hdf5 files, in the sense that the current hdf5 tree would be updated to net0-net1/stat0-stat1/loc0-loc1/cha0-cha1.

Important is to allow for backward compatibility with the old hdf5 structure. probably only possible by allowing for both file structures in the params.yaml.

When remove_response=True, Correlator.pxcorr sometimes acts up, when processing using mpi. This seems to be caused by a process that opens an hdf5 file at CorrelationDataBase._write() and freezes while doing so. Then, the next process wanting to write to that file crashes as files cannot be opened by several processes.

Hello, I collected some data by myself with a three-component seismograph, the original format of the data is DLD, the data can be converted into miniseed format using sololite, but there is no network information in the data, may I ask how to edit the network when I edit and sort the files? Can i define one myself?

Avoid duplicates

• I searched existing issues

Bug Summary

In the examples, if I run the tutorial "Plot a time-dependent spectrogram"

from seismic.plot.plot_spectrum import plot_spct_series
from matplotlib import pyplot as plt

f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
fig = plt.figure(figsize=(9, 7))
plot_spct_series(S, f, t, log_scale=True)

the module sc is not found

NameError Traceback (most recent call last)
Cell In[1], line 4
1 from seismic.plot.plot_spectrum import plot_spct_series
2 from matplotlib import pyplot as plt
----> 4 f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
5 fig = plt.figure(figsize=(9, 7))
6 plot_spct_series(S, f, t, log_scale=True)

NameError: name 'sc' is not defined

Code to Reproduce

from seismic.plot.plot_spectrum import plot_spct_series
from matplotlib import pyplot as plt

f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
fig = plt.figure(figsize=(9, 7))
plot_spct_series(S, f, t, log_scale=True)

Error Traceback

NameError                                 Traceback (most recent call last)
Cell In[1], line 4
      1 from seismic.plot.plot_spectrum import plot_spct_series
      2 from matplotlib import pyplot as plt
----> 4 f, t, S = sc.compute_spectrogram('X9', 'IR1', 'HHE', starttime, endtime, 7200, freq_max=12.5)
      5 fig = plt.figure(figsize=(9, 7))
      6 plot_spct_series(S, f, t, log_scale=True)

NameError: name 'sc' is not defined

SeisMIC Version?

last

Operating System?

OS/X Ventura

Python Version?

Python 3.11.5

Installation Method?

conda

When attempting to load any correlation from julday 360 (which year doesn't matter), the following error occurs:
ValueError: unconverted data remains: 1

This is caused by an issue with obspy's UTCDateTime and should be fixed with their next release (1.23). However, until then it should be possible to implement a dirty fix.