[Antelope] overlapping wf records
Jonathan MacCarthy
jkmacc at nmt.edu
Sat Oct 18 20:09:54 CDT 2008
Hello Antelope users,
I'm new to Antelope, and I've been looking for a command/trick that
might solve a problem I'm having, but with no luck so far. I'm
hoping someone might have a suggestion; any help would be much
appreciated...
With my origin, site, and wfdisc tables, I'm doing the canonical
"extract parrival-30sec to parrival+60sec from wfdisc holdings", but
I'm running into overlap problems. My wfdisc is populated by
downloading P-60sec to P+120sec from IRIS for select events. The
problem is due to a few of those origins being close to each other in
time & space such that my " parrival-30sec to parrival+60sec" call for
the earlier event also grabs from wfdisc records intended to be
associated with the later event. Here's an example for database name
"TA":
The events...
>> dbsubset TA.origin 'orid == 16 || orid == 17' | dbselect - orid
lat lon depth 'strtime(time)'
16 36.1780 141.5450 19.0000 5/07/2008 16:02:02.600
17 36.2260 141.7660 35.0000 5/07/2008 16:03:40.430
The extract call:
>> dbjoin TA.origin TA.site | dbsubset - 'orid==16' | dbjoin - :sta
':parrival()-30::parrival()+60#time::endtime' wfdisc | dbselect - orid
sta chan 'strtime(wfdisc.time)' 'strtime(parrival())'
'strtime(wfdisc.endtime)' 'dir . "/" . dfile' | head
16 SRU BHZ 5/07/2008 16:13:00.965 5/07/2008
16:14:08.829 5/07/2008 16:16:20.715 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_02_02/UU.SRU.BHZ.sac
16 SRU BHE 5/07/2008 16:13:04.090 5/07/2008
16:14:08.829 5/07/2008 16:16:14.715 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_02_02/UU.SRU.BHE.sac
16 SRU BHN 5/07/2008 16:13:05.390 5/07/2008
16:14:08.829 5/07/2008 16:16:19.690 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_02_02/UU.SRU.BHN.sac
16 SRU BHN 5/07/2008 16:14:31.040 5/07/2008
16:14:08.829 5/07/2008 16:17:44.940 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_03_40/UU.SRU.BHN.sac
16 SRU BHZ 5/07/2008 16:14:37.690 5/07/2008
16:14:08.829 5/07/2008 16:18:00.690 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_03_40/UU.SRU.BHZ.sac
16 SRU BHE 5/07/2008 16:14:43.640 5/07/2008
16:14:08.829 5/07/2008 16:17:52.715 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_03_40/UU.SRU.BHE.sac
16 ISCO BHN 5/07/2008 16:10:51.224 5/07/2008
16:14:21.957 5/07/2008 16:18:30.024 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_02_02/US.ISCO.BHN.sac
16 ISCO BHE 5/07/2008 16:11:19.049 5/07/2008
16:14:21.957 5/07/2008 16:18:55.199 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_02_02/US.ISCO.BHE.sac
16 ISCO BHZ 5/07/2008 16:13:22.874 5/07/2008
16:14:21.957 5/07/2008 16:18:08.124 /fs/scratch/jkmacc/Data/CREST/TA/
P/wf/Event_2008_128_16_03_40/US.ISCO.BHZ.sac
...
I've output 3 times: parrival-30sec, parrival, and parrival+60sec, so
that the overlapping wfdisc times are apparent. I've also output the
the "dir" field to demonstrate that the files are from different event
directories.
Is there some way to combine overlapping wfdisc records so that this
doesn't happen, or to otherwise distinguish between the records I want
and the ones I don't?
Thanks for your help!
Best,
Jonathan MacCarthy
Jonathan MacCarthy
PhD Candidate in Seismology
New Mexico Institute of Mining and Technology
Dept. of Earth and Environmental Science, MSEC208
801 Leroy Place
Socorro, NM 87801
ph: 505.835.5418
web: www.ees.nmt.edu/~jkmacc
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://brtt.net/pipermail/antelope_brtt.net/attachments/20081018/8f82e461/attachment.html>
More information about the Antelope
mailing list