INGV station instrument response error

[helle-pedersen]

Dear INGV colleagues,

I am presently doing an analysis using many seismic stations across Europe, using the vertical component.
I am working in obspy, downloading inventory, then downloading data, then simple preprocessing.

All works well for INGV stations**, except for station CSOB.** The program fails when deconvolving with the instrument response.

Here is an example that doesn't work. I am working on this stream:
IV.CSOB..HHZ | 2022-02-01T22:43:07.300000Z - 2022-02-01T22:49:47.300000Z | 100.0 Hz, 40001 samples

There may be errors also on the other components and for other times.

The processing is simple : :
def preproc(st,fmin,fmax,inv) :
st_plot_filt=st.copy()
st_plot_filt=st_plot_filt.detrend("linear")
st_plot_filt.taper(max_percentage=0.04)
st_plot_filt.remove_response(output="VEL", water_level=10, inventory=inv)
st_plot_filt.filter("bandpass", freqmin=fmin, freqmax=fmax)
st_plot_filt.normalize()
return st_plot_filt

It fails at the line
----> 5 st_plot_filt_ingv=preproc(st_final_ingv,fmin,fmax, inventory_ingv)
and looking through the errors, comes with
ValueError: check_channel: Illegal RESP format

Would it be able to fix this issue? It makes it quite a lot more complex to download data from a lot of INGV stations into a single stream, then deconvolving.

I will of course make an exception in my present program, but I'd be grateful if you could drop me a note if/when the issue is fixed.

Just contact me if you have any questions

Best regards

Helle Pedersen

[helle-pedersen]

Btw, for the same time window, BNI must have a wrong instrument response. The deconvolution goes ahead, but emits a warning, and the amplitudes are off scale, by 10000 or so (again, horizontal component)

[petrrr]

Btw, for the same time window, BNI must have a wrong instrument response. The deconvolution goes ahead, but emits a warning, and the amplitudes are off scale, by 10000 or so (again, horizontal component)

@helle-pedersen: We were not able to understand this problem. Are you refering to channels HH?