From re.rezaei at gmail.com Mon Dec 2 12:40:27 2019 From: re.rezaei at gmail.com (Reza Rezaei) Date: Mon, 2 Dec 2019 15:10:27 +0330 Subject: [obspy-users] Variables in header Message-ID: Hello every one, In SAC, we can write couple of variables in the header (T0~T10). Is there any similar approach to write and retrieve variables in Stream or Trace? Best regards, Reza -------------- next part -------------- An HTML attachment was scrubbed... URL: From megies at geophysik.uni-muenchen.de Tue Dec 3 14:17:43 2019 From: megies at geophysik.uni-muenchen.de (Tobias Megies) Date: Tue, 3 Dec 2019 14:17:43 +0100 Subject: [obspy-users] PPSD workflow In-Reply-To: <1616240347.49079.1574741556391@mail.yahoo.com> References: <1616240347.49079.1574741556391.ref@mail.yahoo.com> <1616240347.49079.1574741556391@mail.yahoo.com> Message-ID: Hi Januka, the approach very closely follows the processing workflows described in the references papers / technical reports by McNamara et al. http://docs.obspy.org/citations.html?highlight=mcnamara By default PPSD works on 1-hour slices (like PQLX), but that can be changed. For details please see the source code: https://github.com/obspy/obspy/blob/4c83dfc6e4a6ae4c93c8681349661d67d48505bc/obspy/signal/spectral_estimation.py#L903-L1001 As you can see in there the psd is reused from matplotlib: https://matplotlib.org/3.1.1/api/mlab_api.html#matplotlib.mlab.psd cheers, T On 26/11/2019 05:12, Januka Attanayake wrote: > Hi,? > Does the PPSD class follow the same workflow as that described in the > PQLX manual? i.e., is the day long file truncated to the next lowest > power of two?? > > "First, to significantly improve the Fast Fourier Transform (FFT) speed > ratio, by reducing the number of operations, the number of samples in > the time series, N, is truncated to the next lowest power of two, 217, > leaving N=131,072, thereby reducing the series length such that > Th=3276.8s. Second, in order to further reduce the variance of the final > PSD estimates, each roughly 1-hour time series record is divided into 13 > segments, overlapping by 75 percent, where the length of each new time > series segment is now, Tr=Th/4=819.2s with N=32,768=215." > > Thanks! > > Januka > From megies at geophysik.uni-muenchen.de Tue Dec 3 14:28:38 2019 From: megies at geophysik.uni-muenchen.de (Tobias Megies) Date: Tue, 3 Dec 2019 14:28:38 +0100 Subject: [obspy-users] ObsPy at AGU 2019 Message-ID: Hi all, I'll be at AGU and there will be a lightning talk in an Open Source software session (NS21A) on Tuesday morning: https://agu.confex.com/agu/fm19/meetingapp.cgi/Paper/577013 Also following this oral session there are a few posters from this session (NS21B) so we'll be gathering there and you can get in touch with me there on Tuesday around 10am to noon: https://agu.confex.com/agu/fm19/meetingapp.cgi/Session/78079 We'll also organize an informal gathering for food/drinks in the evening on one day. It's not yet fixed when, but it'll be posted on our Twitter and (hopefully I'll remember) also here. https://twitter.com/obspy Also you can send a direct message on Twitter if you wanna get in touch during AGU. cheers, Tobias From jkmacc at lanl.gov Tue Dec 3 16:10:43 2019 From: jkmacc at lanl.gov (MacCarthy, Jonathan K) Date: Tue, 3 Dec 2019 15:10:43 +0000 Subject: [obspy-users] Variables in header In-Reply-To: References: Message-ID: Hi Reza, Trace objects have a .meta attribute that is essentially just a dictionary, and can hold any keys and values you add to it. Just be careful not to overwrite the core keys, like ?starttime?, ?sampling_rate?, etc. Best, Jon ------------------------------------------------------------------------------------ Jonathan K. MacCarthy Seismoacoustics Team Los Alamos National Laboratory Phone: 505-665-4448 Cell: 505-695-5044 ------------------------------------------------------------------------------------ On Dec 2, 2019, at 4:40 AM, Reza Rezaei > wrote: Hello every one, In SAC, we can write couple of variables in the header (T0~T10). Is there any similar approach to write and retrieve variables in Stream or Trace? Best regards, Reza _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://lists.swapbytes.de/mailman/listinfo/obspy-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From jattanayake at yahoo.com Wed Dec 4 09:17:51 2019 From: jattanayake at yahoo.com (Januka Attanayake) Date: Wed, 4 Dec 2019 08:17:51 +0000 (UTC) Subject: [obspy-users] PPSD workflow In-Reply-To: References: <1616240347.49079.1574741556391.ref@mail.yahoo.com> <1616240347.49079.1574741556391@mail.yahoo.com> Message-ID: <551085015.3127481.1575447471030@mail.yahoo.com> Thanks as always Tobias! Januka Attanayake Research Fellow | Earthquake Seismology Homepage:?https://sites.google.com/site/janukaattanayake School of Earth Sciences | McCoy Bldg. 200 University of Melbourne | Parkville 3010 VIC Australia On Wednesday, 4 December 2019, 2:05:51 am AEDT, Tobias Megies wrote: Hi Januka, the approach very closely follows the processing workflows described in the references papers / technical reports by McNamara et al. http://docs.obspy.org/citations.html?highlight=mcnamara By default PPSD works on 1-hour slices (like PQLX), but that can be changed. For details please see the source code: https://github.com/obspy/obspy/blob/4c83dfc6e4a6ae4c93c8681349661d67d48505bc/obspy/signal/spectral_estimation.py#L903-L1001 As you can see in there the psd is reused from matplotlib: https://matplotlib.org/3.1.1/api/mlab_api.html#matplotlib.mlab.psd cheers, T On 26/11/2019 05:12, Januka Attanayake wrote: > Hi,? > Does the PPSD class follow the same workflow as that described in the > PQLX manual? i.e., is the day long file truncated to the next lowest > power of two?? > > "First, to significantly improve the Fast Fourier Transform (FFT) speed > ratio, by reducing the number of operations, the number of samples in > the time series, N, is truncated to the next lowest power of two, 217, > leaving N=131,072, thereby reducing the series length such that > Th=3276.8s. Second, in order to further reduce the variance of the final > PSD estimates, each roughly 1-hour time series record is divided into 13 > segments, overlapping by 75 percent, where the length of each new time > series segment is now, Tr=Th/4=819.2s with N=32,768=215." > > Thanks! > > Januka > _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://lists.swapbytes.de/mailman/listinfo/obspy-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From vivekgbabu at gmail.com Thu Dec 5 19:25:17 2019 From: vivekgbabu at gmail.com (Vivek.G.Babu) Date: Thu, 5 Dec 2019 23:55:17 +0530 Subject: [obspy-users] obspy.clients : Help needed Message-ID: Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- *Vivek G Bab* -------------- next part -------------- An HTML attachment was scrubbed... URL: From calum.chamberlain at vuw.ac.nz Sun Dec 8 03:15:02 2019 From: calum.chamberlain at vuw.ac.nz (Calum Chamberlain) Date: Sun, 8 Dec 2019 02:15:02 +0000 Subject: [obspy-users] obspy.clients : Help needed In-Reply-To: References: Message-ID: Hi Vivek, You should just be able to cat.write("catalog.out", format="NORDIC") This is the general format for writing Catalog objects from obspy. See the docs page here https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html Note there have been a lot of improvements for NORDIC support for the upcoming version 1.2.0 that are in the current obspy master, so if you find that the files are missing something you need, then let us know and we can help you install the current master branch. Calum. Calum Chamberlain Postdoctoral Research Fellow in Seismology and Fault Mechanics E: calum.chamberlain at vuw.ac.nz P: +64 420 463 6353 On Dec 5 2019, at 10:25 am, Vivek.G.Babu wrote: Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- Vivek G Bab _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 -------------- next part -------------- An HTML attachment was scrubbed... URL: From calum.chamberlain at vuw.ac.nz Sun Dec 8 03:15:02 2019 From: calum.chamberlain at vuw.ac.nz (Calum Chamberlain) Date: Sun, 8 Dec 2019 02:15:02 +0000 Subject: [obspy-users] obspy.clients : Help needed In-Reply-To: References: Message-ID: Hi Vivek, You should just be able to cat.write("catalog.out", format="NORDIC") This is the general format for writing Catalog objects from obspy. See the docs page here https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html Note there have been a lot of improvements for NORDIC support for the upcoming version 1.2.0 that are in the current obspy master, so if you find that the files are missing something you need, then let us know and we can help you install the current master branch. Calum. Calum Chamberlain Postdoctoral Research Fellow in Seismology and Fault Mechanics E: calum.chamberlain at vuw.ac.nz P: +64 420 463 6353 On Dec 5 2019, at 10:25 am, Vivek.G.Babu wrote: Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- Vivek G Bab _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 -------------- next part -------------- An HTML attachment was scrubbed... URL: From vivekgbabu at gmail.com Sun Dec 8 10:37:53 2019 From: vivekgbabu at gmail.com (Vivek.G.Babu) Date: Sun, 8 Dec 2019 15:07:53 +0530 Subject: [obspy-users] obspy-users Digest, Vol 109, Issue 3 In-Reply-To: References: Message-ID: Dear Calum, Thank you so much for your suggestion. But I already tried that. I am getting an error like: Traceback (most recent call last): File "", line 1, in File "E:\WIHG\Programming\python\others\obspy\core\event\catalog.py", line 509, in write return write_format(self, filename, **kwargs) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 931, in write_select string_io=select, high_accuracy=high_accuracy) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 1172, in _write_nordic newpicks = '\n'.join(nordpick(event, high_accuracy=high_accuracy)) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 1470, in nordpick polarity.rjust(1) + ' ') TypeError: can only concatenate str (not "NoneType") to str "E:\WIHG\Programming\python\others" is my current working directory, On Sun, 8 Dec 2019 at 08:05, wrote: > Send obspy-users mailing list submissions to > obspy-users at lists.swapbytes.de > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.swapbytes.de/mailman/listinfo/obspy-users > or, via email, send a message with subject or body 'help' to > obspy-users-request at lists.swapbytes.de > > You can reach the person managing the list at > obspy-users-owner at lists.swapbytes.de > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of obspy-users digest..." > > > Today's Topics: > > 1. obspy.clients : Help needed (Vivek.G.Babu) > 2. Re: obspy.clients : Help needed (Calum Chamberlain) > 3. Re: obspy.clients : Help needed (Calum Chamberlain) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 5 Dec 2019 23:55:17 +0530 > From: "Vivek.G.Babu" > To: obspy-users at lists.swapbytes.de > Subject: [obspy-users] obspy.clients : Help needed > Message-ID: > HXw at mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Hi, > I am very new to obspy package, I am not sure this is the right place to > ask this question. > I am trying to download events from ISC (includearrivals =True), And I am > getting an object, But I don't know how exactly I could make an output file > with the arrival phases in NORDIC format from this.!?!? > > Any suggestions? > > > eg: > > client=Client("ISC") > > cat = client.get_events(eventid=867982,includearrivals=True) > > how to make an output file in NORDIC (or text format) which contains event > parameters and the phase arrivals? > > > THANK YOU > > > > > -- > *Vivek G Bab* > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://lists.swapbytes.de/pipermail/obspy-users/attachments/20191205/b3dbd747/attachment-0001.html > > > > ------------------------------ > > Message: 2 > Date: Sun, 8 Dec 2019 02:15:02 +0000 > From: Calum Chamberlain > To: "users at obspy.org" > Cc: "obspy-users at lists.swapbytes.de" > Subject: Re: [obspy-users] obspy.clients : Help needed > Message-ID: > Content-Type: text/plain; charset="utf-8" > > Hi Vivek, > > You should just be able to cat.write("catalog.out", format="NORDIC") This > is the general format for writing Catalog objects from obspy. See the docs > page here > https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html > Note there have been a lot of improvements for NORDIC support for the > upcoming version 1.2.0 that are in the current obspy master, so if you find > that the files are missing something you need, then let us know and we can > help you install the current master branch. > > Calum. > > Calum Chamberlain > Postdoctoral Research Fellow in Seismology and Fault Mechanics > E: calum.chamberlain at vuw.ac.nz > P: +64 420 463 6353 > On Dec 5 2019, at 10:25 am, Vivek.G.Babu wrote: > Hi, > I am very new to obspy package, I am not sure this is the right place to > ask this question. > I am trying to download events from ISC (includearrivals =True), And I am > getting an object, But I don't know how exactly I could make an output file > with the arrival phases in NORDIC format from this.!?!? > > Any suggestions? > > > eg: > client=Client("ISC") > cat = client.get_events(eventid=867982,includearrivals=True) > how to make an output file in NORDIC (or text format) which contains event > parameters and the phase arrivals? > > THANK YOU > > > > -- > Vivek G Bab > > > _______________________________________________ > obspy-users mailing list > obspy-users at lists.swapbytes.de > > https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://lists.swapbytes.de/pipermail/obspy-users/attachments/20191208/79fab50f/attachment.html > > > > ------------------------------ > > Message: 3 > Date: Sun, 8 Dec 2019 02:15:02 +0000 > From: Calum Chamberlain > To: "users at obspy.org" > Cc: "obspy-users at lists.swapbytes.de" > Subject: Re: [obspy-users] obspy.clients : Help needed > Message-ID: > Content-Type: text/plain; charset="utf-8" > > Hi Vivek, > > You should just be able to cat.write("catalog.out", format="NORDIC") This > is the general format for writing Catalog objects from obspy. See the docs > page here > https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html > Note there have been a lot of improvements for NORDIC support for the > upcoming version 1.2.0 that are in the current obspy master, so if you find > that the files are missing something you need, then let us know and we can > help you install the current master branch. > > Calum. > > Calum Chamberlain > Postdoctoral Research Fellow in Seismology and Fault Mechanics > E: calum.chamberlain at vuw.ac.nz > P: +64 420 463 6353 > On Dec 5 2019, at 10:25 am, Vivek.G.Babu wrote: > Hi, > I am very new to obspy package, I am not sure this is the right place to > ask this question. > I am trying to download events from ISC (includearrivals =True), And I am > getting an object, But I don't know how exactly I could make an output file > with the arrival phases in NORDIC format from this.!?!? > > Any suggestions? > > > eg: > client=Client("ISC") > cat = client.get_events(eventid=867982,includearrivals=True) > how to make an output file in NORDIC (or text format) which contains event > parameters and the phase arrivals? > > THANK YOU > > > > -- > Vivek G Bab > > > _______________________________________________ > obspy-users mailing list > obspy-users at lists.swapbytes.de > > https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://lists.swapbytes.de/pipermail/obspy-users/attachments/20191208/79fab50f/attachment-0001.html > > > > ------------------------------ > > _______________________________________________ > obspy-users mailing list > obspy-users at lists.swapbytes.de > https://lists.swapbytes.de/mailman/listinfo/obspy-users > > End of obspy-users Digest, Vol 109, Issue 3 > ******************************************* > -- *Vivek G Babu* *Research Fellow* *Geophysics group* *Wadia Institute of Himalayan Geology* *Dehradun - 248001* *Mob : +919446270718* *ORCID : 0000-0001-7068-0302* -------------- next part -------------- An HTML attachment was scrubbed... URL: From megies at geophysik.uni-muenchen.de Tue Dec 10 15:30:06 2019 From: megies at geophysik.uni-muenchen.de (Tobias Megies) Date: Tue, 10 Dec 2019 06:30:06 -0800 Subject: [obspy-users] ObsPy at AGU 2019 In-Reply-To: References: Message-ID: <69c18d9e-6673-e6cf-76fa-da1a0fa12611@geophysik.uni-muenchen.de> Hi all, just a reminder/follow-up: You can meet me at the end of the oral session NS21A today at around 10am, or after that at the accompanying poster session NS21B around 10-12am. In the evening you can join us for drinks/dinner, let's meet at 6:15pm at the poster pickup booth in Moscone South (when you walk in from the street at the left side of the entrance hall). If you wanna, please drop me a private email so that we know whom to wait for and maybe I can plan ahead and reserve some table. best, Tobias On 03/12/2019 05:28, Tobias Megies wrote: > Hi all, > > I'll be at AGU and there will be a lightning talk in an Open Source > software session (NS21A) on Tuesday morning: > https://agu.confex.com/agu/fm19/meetingapp.cgi/Paper/577013 > > Also following this oral session there are a few posters from this > session (NS21B) so we'll be gathering there and you can get in touch > with me there on Tuesday around 10am to noon: > https://agu.confex.com/agu/fm19/meetingapp.cgi/Session/78079 > > We'll also organize an informal gathering for food/drinks in the evening > on one day. It's not yet fixed when, but it'll be posted on our Twitter > and (hopefully I'll remember) also here. > > https://twitter.com/obspy > > Also you can send a direct message on Twitter if you wanna get in touch > during AGU. > > cheers, > Tobias > _______________________________________________ > obspy-users mailing list > obspy-users at lists.swapbytes.de > https://lists.swapbytes.de/mailman/listinfo/obspy-users > -- EMail: tobias.megies at lmu.de WWW: https://www.geophysik.uni-muenchen.de Geophysikalisches Observatorium Ludwigsh?he 8 82256 F?rstenfeldbruck Tel: +49 (0) 89 2180-73981 Fax: +49 (0) 89 2180-73970 Ludwig-Maximilians-Universit?t Department f?r Geo- und Umweltwissenschaften Sektion Geophysik Theresienstrasse 41/IV 80333 M?nchen Tel: +49 (0) 89 2180-4326 Fax: +49 (0) 89 2180-4205 From calum.chamberlain at vuw.ac.nz Wed Dec 11 02:46:05 2019 From: calum.chamberlain at vuw.ac.nz (Calum Chamberlain) Date: Wed, 11 Dec 2019 01:46:05 +0000 Subject: [obspy-users] obspy-users Digest, Vol 109, Issue 3 In-Reply-To: References: , Message-ID: Hi Viveck, Would you be able to send me the catalog (maybe as a quakeml). It may be something that is fixed in the current master, but it's worth me checking before making you go to the effort of installing that. C CJ Chamberlain, out of office ________________________________ From: obspy-users-bounces at lists.swapbytes.de on behalf of Vivek.G.Babu Sent: Sunday, December 8, 2019 1:37:53 AM To: obspy-users at lists.swapbytes.de Subject: Re: [obspy-users] obspy-users Digest, Vol 109, Issue 3 Dear Calum, Thank you so much for your suggestion. But I already tried that. I am getting an error like: Traceback (most recent call last): File "", line 1, in File "E:\WIHG\Programming\python\others\obspy\core\event\catalog.py", line 509, in write return write_format(self, filename, **kwargs) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 931, in write_select string_io=select, high_accuracy=high_accuracy) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 1172, in _write_nordic newpicks = '\n'.join(nordpick(event, high_accuracy=high_accuracy)) File "E:\WIHG\Programming\python\others\obspy\io\nordic\core.py", line 1470, in nordpick polarity.rjust(1) + ' ') TypeError: can only concatenate str (not "NoneType") to str "E:\WIHG\Programming\python\others" is my current working directory, On Sun, 8 Dec 2019 at 08:05, > wrote: Send obspy-users mailing list submissions to obspy-users at lists.swapbytes.de To subscribe or unsubscribe via the World Wide Web, visit https://lists.swapbytes.de/mailman/listinfo/obspy-users or, via email, send a message with subject or body 'help' to obspy-users-request at lists.swapbytes.de You can reach the person managing the list at obspy-users-owner at lists.swapbytes.de When replying, please edit your Subject line so it is more specific than "Re: Contents of obspy-users digest..." Today's Topics: 1. obspy.clients : Help needed (Vivek.G.Babu) 2. Re: obspy.clients : Help needed (Calum Chamberlain) 3. Re: obspy.clients : Help needed (Calum Chamberlain) ---------------------------------------------------------------------- Message: 1 Date: Thu, 5 Dec 2019 23:55:17 +0530 From: "Vivek.G.Babu" > To: obspy-users at lists.swapbytes.de Subject: [obspy-users] obspy.clients : Help needed Message-ID: > Content-Type: text/plain; charset="utf-8" Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- *Vivek G Bab* -------------- next part -------------- An HTML attachment was scrubbed... URL: > ------------------------------ Message: 2 Date: Sun, 8 Dec 2019 02:15:02 +0000 From: Calum Chamberlain > To: "users at obspy.org" > Cc: "obspy-users at lists.swapbytes.de" > Subject: Re: [obspy-users] obspy.clients : Help needed Message-ID: > Content-Type: text/plain; charset="utf-8" Hi Vivek, You should just be able to cat.write("catalog.out", format="NORDIC") This is the general format for writing Catalog objects from obspy. See the docs page here https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html Note there have been a lot of improvements for NORDIC support for the upcoming version 1.2.0 that are in the current obspy master, so if you find that the files are missing something you need, then let us know and we can help you install the current master branch. Calum. Calum Chamberlain Postdoctoral Research Fellow in Seismology and Fault Mechanics E: calum.chamberlain at vuw.ac.nz> P: +64 420 463 6353 On Dec 5 2019, at 10:25 am, Vivek.G.Babu > wrote: Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- Vivek G Bab _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 -------------- next part -------------- An HTML attachment was scrubbed... URL: > ------------------------------ Message: 3 Date: Sun, 8 Dec 2019 02:15:02 +0000 From: Calum Chamberlain > To: "users at obspy.org" > Cc: "obspy-users at lists.swapbytes.de" > Subject: Re: [obspy-users] obspy.clients : Help needed Message-ID: > Content-Type: text/plain; charset="utf-8" Hi Vivek, You should just be able to cat.write("catalog.out", format="NORDIC") This is the general format for writing Catalog objects from obspy. See the docs page here https://docs.obspy.org/packages/autogen/obspy.core.event.Catalog.write.html Note there have been a lot of improvements for NORDIC support for the upcoming version 1.2.0 that are in the current obspy master, so if you find that the files are missing something you need, then let us know and we can help you install the current master branch. Calum. Calum Chamberlain Postdoctoral Research Fellow in Seismology and Fault Mechanics E: calum.chamberlain at vuw.ac.nz> P: +64 420 463 6353 On Dec 5 2019, at 10:25 am, Vivek.G.Babu > wrote: Hi, I am very new to obspy package, I am not sure this is the right place to ask this question. I am trying to download events from ISC (includearrivals =True), And I am getting an object, But I don't know how exactly I could make an output file with the arrival phases in NORDIC format from this.!?!? Any suggestions? eg: client=Client("ISC") cat = client.get_events(eventid=867982,includearrivals=True) how to make an output file in NORDIC (or text format) which contains event parameters and the phase arrivals? THANK YOU -- Vivek G Bab _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://apc01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.swapbytes.de%2Fmailman%2Flistinfo%2Fobspy-users&data=02%7C01%7Ccalum.chamberlain%40vuw.ac.nz%7C1ef087d2616e45b29d7c08d77b83609b%7Ccfe63e236951427e8683bb84dcf1d20c%7C0%7C0%7C637113676578869585&sdata=ki3OH34DYnXu%2FWMkkmrHP9aDRg8HvoF7awrCF1z%2BkIY%3D&reserved=0 -------------- next part -------------- An HTML attachment was scrubbed... URL: > ------------------------------ _______________________________________________ obspy-users mailing list obspy-users at lists.swapbytes.de https://lists.swapbytes.de/mailman/listinfo/obspy-users End of obspy-users Digest, Vol 109, Issue 3 ******************************************* -- Vivek G Babu Research Fellow Geophysics group Wadia Institute of Himalayan Geology Dehradun - 248001 Mob : +919446270718 ORCID : 0000-0001-7068-0302 -------------- next part -------------- An HTML attachment was scrubbed... URL: From racine at sed.ethz.ch Mon Dec 16 14:53:37 2019 From: racine at sed.ethz.ch (Roman Racine) Date: Mon, 16 Dec 2019 14:53:37 +0100 Subject: [obspy-users] Resource leak in RoutingClient Message-ID: <72502763-d64e-14be-6ece-a374d3553b60@sed.ethz.ch> Hi all, I am trying to replace arclink requests by fdsn requests in my software. I have tried using obspy.clients.fdsn.RoutingClient. However, when I make a large number of requests, what I usually do, lsof shows more and more resources allocated by my script which are never freed again. Here a trivial example: ****** import obspy.clients.fdsn client = obspy.clients.fdsn.RoutingClient('eida-routing') starttime = obspy.core.UTCDateTime()-86400 endtime = starttime + 600 for i in xrange(1,1500): ??? st = client.get_waveforms(network='CH',station='BALST',location='*',channel='HHZ',starttime=starttime,endtime=endtime) ??? print st ******** lsof|grep |wc -l shows a rapidly growing number. I'm currently using the ubuntu package for bionic: python-obspy????????????????????????????????? 1.1.1-1~bionic Kind regards Roman Racine -- Schweizerischer Erdbebendienst (SED) ETH Z?rich Roman Racine Software Engineer NO, FO 69.6 Sonneggstrasse 5 8092 Z?rich Telefon: +41 44 633 24 69 E-Mail: racine at sed.ethz.ch www.seismo.ethz.ch From racine at sed.ethz.ch Tue Dec 17 14:35:29 2019 From: racine at sed.ethz.ch (Roman Racine) Date: Tue, 17 Dec 2019 14:35:29 +0100 Subject: [obspy-users] Resource leak in RoutingClient In-Reply-To: <1576589585.1555.0@mail.univie.ac.at> References: <72502763-d64e-14be-6ece-a374d3553b60@sed.ethz.ch> <1576589585.1555.0@mail.univie.ac.at> Message-ID: <479096ba-7995-ce78-8fb2-b80851eea033@sed.ethz.ch> Hi Florian, this is exactly the same problem. I wonder why I did not find it when searching for any clues yesterday. Regards Roman On 17/12/2019 14.33, Florian Fuchs wrote: > Hallo Roman, > > I think this issues was already once discussed here: > https://github.com/obspy/obspy/issues/2342 > > Does that look like your situation? > > Best regards, > Florian > _____________________________________ Dr. Florian Fuchs Department of > Meteorology and Geophysics University of Vienna UZA II, 2D508 > Althanstra?e 14 1090 Vienna Austria Phone: +43 1 4277 53726 Mobile: > +43 664 60277 53726 www: http://imgw.univie.ac.at > _____________________________________ > > Am Mo, 16. Dez, 2019 um 2:53 NACHMITTAGS schrieb Roman Racine > : >> Hi all, I am trying to replace arclink requests by fdsn requests in >> my software. I have tried using obspy.clients.fdsn.RoutingClient. >> However, when I make a large number of requests, what I usually do, >> lsof shows more and more resources allocated by my script which are >> never freed again. Here a trivial example: ****** import >> obspy.clients.fdsn client = >> obspy.clients.fdsn.RoutingClient('eida-routing') starttime = >> obspy.core.UTCDateTime()-86400 endtime = starttime + 600 for i in >> xrange(1,1500): ??? st = >> client.get_waveforms(network='CH',station='BALST',location='*',channel='HHZ',starttime=starttime,endtime=endtime) >> ??? print st ******** lsof|grep |wc -l shows a rapidly growing >> number. I'm currently using the ubuntu package for bionic: >> python-obspy????????????????????????????????? 1.1.1-1~bionic Kind >> regards Roman Racine >> -- >> Schweizerischer Erdbebendienst (SED) ETH Z?rich Roman Racine Software >> Engineer NO, FO 69.6 Sonneggstrasse 5 8092 Z?rich Telefon: +41 44 633 >> 24 69 E-Mail: racine at sed.ethz.ch >> www.seismo.ethz.ch >> _______________________________________________ obspy-users mailing >> list obspy-users at lists.swapbytes.de >> >> https://lists.swapbytes.de/mailman/listinfo/obspy-users -- Schweizerischer Erdbebendienst (SED) ETH Z?rich Roman Racine Software Engineer NO, FO 69.6 Sonneggstrasse 5 8092 Z?rich Telefon: +41 44 633 24 69 E-Mail: racine at sed.ethz.ch www.seismo.ethz.ch From florian.fuchs at univie.ac.at Tue Dec 17 14:33:05 2019 From: florian.fuchs at univie.ac.at (Florian Fuchs) Date: Tue, 17 Dec 2019 14:33:05 +0100 Subject: [obspy-users] Resource leak in RoutingClient In-Reply-To: <72502763-d64e-14be-6ece-a374d3553b60@sed.ethz.ch> References: <72502763-d64e-14be-6ece-a374d3553b60@sed.ethz.ch> Message-ID: <1576589585.1555.0@mail.univie.ac.at> Hallo Roman, I think this issues was already once discussed here: https://github.com/obspy/obspy/issues/2342 Does that look like your situation? Best regards, Florian _____________________________________ Dr. Florian Fuchs Department of Meteorology and Geophysics University of Vienna UZA II, 2D508 Althanstra?e 14 1090 Vienna Austria Phone: +43 1 4277 53726 Mobile: +43 664 60277 53726 www: http://imgw.univie.ac.at _____________________________________ Am Mo, 16. Dez, 2019 um 2:53 NACHMITTAGS schrieb Roman Racine : > Hi all, > > I am trying to replace arclink requests by fdsn requests in my > software. I have tried using obspy.clients.fdsn.RoutingClient. > However, when I make a large number of requests, what I usually do, > lsof shows more and more resources allocated by my script which are > never freed again. > > Here a trivial example: > ****** > import obspy.clients.fdsn > > client = obspy.clients.fdsn.RoutingClient('eida-routing') > > starttime = obspy.core.UTCDateTime()-86400 > endtime = starttime + 600 > for i in xrange(1,1500): > st = > client.get_waveforms(network='CH',station='BALST',location='*',channel='HHZ',starttime=starttime,endtime=endtime) > print st > ******** > > lsof|grep |wc -l > > shows a rapidly growing number. I'm currently using the ubuntu > package for bionic: > > python-obspy 1.1.1-1~bionic > > Kind regards > > Roman Racine > > -- > Schweizerischer Erdbebendienst (SED) > ETH Z?rich > > Roman Racine > Software Engineer > > NO, FO 69.6 > Sonneggstrasse 5 > 8092 Z?rich > Telefon: +41 44 633 24 69 > E-Mail: racine at sed.ethz.ch > www.seismo.ethz.ch > > _______________________________________________ > obspy-users mailing list > obspy-users at lists.swapbytes.de > https://lists.swapbytes.de/mailman/listinfo/obspy-users -------------- next part -------------- An HTML attachment was scrubbed... URL: From rfkartal at gmail.com Thu Dec 26 15:20:54 2019 From: rfkartal at gmail.com (Recai F. KARTAL) Date: Thu, 26 Dec 2019 17:20:54 +0300 Subject: [obspy-users] Connectin seedlink Message-ID: *Dear Authorized,* *I have a question related to seedlink connection. * *My connection code is below:* from obspy.clients.seedlink.easyseedlink import * def handle_data(trace): print(trace.stats.station) print(trace.stats.endtime) client = create_client('remote.org:18000', handle_data) station = ?ABC? net = ?CD? client.select_stream(net, station, 'EL?') client.run() *My goal is to connect to the seedlink server and do station query. I'm linking with the code above. But I can't get a return from handle_data. * *If the station is coming to the seedlink server, I want to do the same query for another station after reading the end time and receiving this information.* *How can do it? can you help me?* *Best regards,* *Recai.* -- *Recai F. KARTAL* *Jeofizik Y?k. M?hendisi / Geophysical Engineer (M. Sc.)* **************************************************************************************************************************************** T. C. ??i?leri Bakanl??? Afet ve Acil Durum Y?netimi Ba?kanl??? / Ministry of Interior Disaster and Emergency Management Presidency Deprem Dairesi / Earthquake Department ?niversiteler Mah. Dumlup?nar Bul. No:159 06800 ?ankaya-ANKARA-TURKEY Phone: +90 312 258 23 23-2153 **************************************************************************************************************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: From gatti.flpp at gmail.com Sat Dec 28 18:22:13 2019 From: gatti.flpp at gmail.com (Filippo) Date: Sat, 28 Dec 2019 18:22:13 +0100 Subject: [obspy-users] TauP resizing Message-ID: Dear users does any of you know why the taup package provides me the following warning? Resizing a TauP array inplace failed due to the existence of other references to the array, creating a new array when asking for arrivals for this model: model = TauPyModel(model="ak135f_no_mud") ?phase_list=["P","S","pP","sS","p","s"] Thanks Filippo -------------- next part -------------- An HTML attachment was scrubbed... URL: