[Antelope] Fwd: overlapping wf records

Sun Oct 19 04:20:09 CDT 2008

Begin forwarded message:

> From: "Dr. Kent Lindquist" <kent at lindquistconsulting.com>
> Date: October 19, 2008 8:18:46 PM GMT+11:00
> To: Jonathan MacCarthy <jkmacc at nmt.edu>
> Cc: antelope at brtt.net
> Subject: Re: [Antelope] overlapping wf records
>
> 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
>

--
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/54d45851/attachment-0001.html>