Skip to content

Error reading 1s low-latency h(t) GWF files for LIGO

The following RuntimeError occurs when attempting to read multiple seconds of data from a series of 1s files at LIGO:

>>> from gwpy.timeseries import TimeSeries
>>> print(TimeSeries.find('H1:GDS-CALIB_STRAIN_INJ1_O3Replay', 1331537687, 1331537689, frametype='H1_O3ReplayMDC_llhoft', verbose=True))
/home/duncan.macleod/git/gwpy-fork/gwpy/signal/filter_design.py:33: DeprecationWarning: Please use `LinearTimeInvariant` from the `scipy.signal` namespace, the `scipy.signal.ltisys` namespace is deprecated.
  from scipy.signal.ltisys import LinearTimeInvariant
Reading H1_O3ReplayMDC_llhoft frames: |                                         | 0/2 (  0%) ETA ?     XLAL Error - XLALFrStreamGetREAL8TimeSeries (LALFrStreamReadTS_source.c:125): End of frame stream while 16384 points remain to be read
XLAL Error - XLALFrStreamGetREAL8TimeSeries (LALFrStreamReadTS_source.c:125): I/O error
XLAL Error - XLALFrStreamReadREAL8TimeSeries (LALFrStreamReadTS_source.c:236): Internal function call failed: I/O error
Traceback (most recent call last):
  File "<string>", line 1, in <module>
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/core.py", line 537, in find
    return cls.DictClass.find(
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/core.py", line 1275, in find
    new = cls.read(cache, names, start=start, end=end, pad=pad,
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/core.py", line 897, in read
    return timeseries_reader(cls, source, *args, **kwargs)
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/io/core.py", line 50, in read
    return io_read_multi(joiner, cls, source, *args, **kwargs)
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/io/mp.py", line 90, in read_multi
    output = mp_utils.multiprocess_with_queues(
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/utils/mp.py", line 123, in multiprocess_with_queues
    return list(map(_inner, inputs))
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/utils/mp.py", line 118, in _inner
    return func(x)
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/io/mp.py", line 119, in _read_single_file
    return fobj, io_read(cls, fobj, *args, **kwargs)
  File "/home/duncan.macleod/.conda/envs/gwcelery/lib/python3.9/site-packages/astropy/io/registry/compat.py", line 43, in wrapper
    return getattr(registry, method_name)(*args, **kwargs)
  File "/home/duncan.macleod/.conda/envs/gwcelery/lib/python3.9/site-packages/astropy/io/registry/core.py", line 199, in read
    data = reader(*args, **kwargs)
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/io/gwf/__init__.py", line 405, in read_
    return reader(*args, **kwargs)
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/io/gwf/__init__.py", line 272, in read_timeseriesdict
    out.append(libread_(src, channels, start=start, end=end,
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/io/gwf/lalframe.py", line 156, in read
    _read_channel(stream, str(name), start=start, duration=duration),
  File "/home/duncan.macleod/git/gwpy-fork/gwpy/timeseries/io/gwf/lalframe.py", line 171, in _read_channel
    return reader(stream, channel, start, duration, 0)
RuntimeError: Internal function call failed: I/O error
Reading H1_O3ReplayMDC_llhoft frames: |████████████████████▌                    | 1/2 ( 50%) ETA 00:00