mass_downloader ignoring minimum_length
mikehagerty opened this issue · comments
mikehagerty commented
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
Tobias Megies commented
Do you have self-contained example code?