obspy / obspy

ObsPy: A Python Toolbox for seismology/seismological observatories.

Home Page:https://www.obspy.org

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

mass_downloader ignoring minimum_length

mikehagerty opened this issue · comments

Avoid duplicates

  • I searched existing issues

Bug Summary

I'm running mass_downloader with minimum_length, e.g:
2024-02-02 12:17:30,960 [ INFO] retrieve_mass_download:retrieve_mass_download: _restrictions={'minimum_length': 0.95, 'channel_priorities': ['BH[ZNE12]',

And asking for a time window that should be 14400 pnts for dt=.025 s
Nevertheless, several downloaded channels have npts < 10000
which is clearly fewer than 95% of 14400.

Thanks!

Code to Reproduce

My code is hitting the following fdsnws waveform servers:
  waveform_servers = ["https://service.ncedc.org", "https://service.scedc.caltech.edu"]

To request data for the following event:
      <origin publicID="quakeml:nc.anss.org/Origin/NC/14320789" catalog:datasource="nc" catalog:dataid="nc73997351" catalog:eventsource="nc" catalog:eventid="73997351">
        <time> 2024-02-01T17:32:54.300000Z
        <latitude> 36.6861667
        <longitude> -121.3251667

It is requesting 6 mins of data: starttime:2024-02-01T17:31:54.300000Z endtime:2024-02-01T17:37:54.300000Z

Here is where mass_downloader is called:

    minradius = min_dist_km / 111.19
    maxradius = max_dist_km / 111.19
    domain = CircularDomain(latitude=evla, longitude=evlo, minradius=minradius, maxradius=maxradius)

    _restrictions = Restrictions(starttime=starttime, endtime=endtime, **restrictions)

    logger.info("retrieve_mass_download: _restrictions=%s" % restrictions)

    providers = []
    for server in waveform_servers:
        providers.append(Client(base_url=server))

    mdl = MassDownloader(providers=providers)
    mdl.download(domain, _restrictions, mseed_storage=mseed_storage, stationxml_storage=stationxml_storage)

Here is some log info:

2024-02-02 11:51:34,438 [   INFO] retrieve_mass_download:retrieve_mass_download: starttime:2024-02-01T17:31:54.300000Z endtime:2024-02-01T17:37:54.300000Z evla:36.6861667 evlo:-121.3251667 min_dist:60 max_dist:900
2024-02-02 11:51:34,438 [   INFO] retrieve_mass_download:retrieve_mass_download: _restrictions={'minimum_length': 0.95, 'channel_priorities': ['BH[ZNE12]', 'HH[ZNE12]'], 'reject_channels_with_gaps': True, 'sanitize': True, 'exclude_stations': ...}

Error Traceback

No response

ObsPy Version?

1.4

Operating System?

OSX

Python Version?

3.10.11

Installation Method?

conda

Do you have self-contained example code?