[Antelope] overlapping wf records

Sun Oct 19 04:18:46 CDT 2008

Hello Jonathan,

Have you tried the trexcerpt(1) command? It takes a small investment  
to learn how to use it for the more complicated operations, however  
it's pretty powerful and may be able to get you through this  
extraction process more cleanly. I don't see any evidence in your  
email below that you're already using it, but you should at least  
consider it.

Best regards,
Kent

On Oct 19, 2008, at 12:09 PM, Jonathan MacCarthy wrote:

> 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
>
> _______________________________________________
> Antelope mailing list
> Antelope at brtt.net
> http://brtt.net/mailman/listinfo/antelope_brtt.net

--
Dr. Kent Lindquist                     kent at lindquistconsulting.com
Lindquist Consulting, Inc.
59 College Rd. Suite #7
Fairbanks, AK 99701              Phone/FAX 907-457-2374

http://www.lindquistconsulting.com

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://brtt.net/pipermail/antelope_brtt.net/attachments/20081019/1d0b3c64/attachment.html>