From mdsilveira at fc.ul.pt Mon Jan 4 09:27:13 2010 From: mdsilveira at fc.ul.pt (=?iso-8859-1?q?Gra=E7a_Silveira?=) Date: Mon, 4 Jan 2010 17:27:13 +0000 Subject: [SAC-HELP] plotinng using colors Message-ID: <201001041727.14019.mdsilveira@fc.ul.pt> Dear all, I need to plot several times the same recordsection (with ~250 files, and using SSS) in the following way: Need to make a cycle were 1) The first trace is red (for example) and all the other ones stay black Then I make the decision f I want to delete it or not 2) The second trace is red and all the other ones stay black Similarly I make the decision 3) The third trace is red and all the other ones stay black -........ ETC ETC I know everything else (decisions, SSS, etc), the only part is the color part. I have already used color increment option, but with less files and I didn't need to specify the colors like that. The point here is in order to better do the decision I would like to have all the files plotted with a specific one enhanced. Thanks Gra?a Silveira -- Gra?a M. Silveira Centro de Geof?sica da Universidade de Lisboa Edif?cio C8 Campo Grande 1749-016 LISBOA tel: +351217500812 email: mdsilveira at fc.ul.pt From fernandez.denise at gmail.com Mon Jan 11 15:05:43 2010 From: fernandez.denise at gmail.com (Denise Fernandez) Date: Tue, 12 Jan 2010 12:05:43 +1300 Subject: [SAC-HELP] how to pick amplitudes Message-ID: <895c839c1001111505h4a9cb477l2798539725dd74a1@mail.gmail.com> Dear all, I'm trying to pick the amplitudes of the s and p arrivals. What I've been done so far is picking them using the command ppk and then positioning the cursor and pressing m. That draws a little box around the portion of the graph that contains the arrival. I thought SAC will store the info in the header but I can't find any variable with information of this specific portion of the graph. I'm new using SAC so I guess I am missing something important here. Do I have to set a blackboard variable before hand? I am not looking for the depmin and depmax just for the amplitude value of the p, s and a pre-s phase. Any input will be very appreciated Cheers, Denise Fernandez -------------- next part -------------- An HTML attachment was scrubbed... URL: From okaven at usgs.gov Tue Jan 12 15:42:14 2010 From: okaven at usgs.gov (J. Ole Kaven) Date: Tue, 12 Jan 2010 15:42:14 -0800 Subject: [SAC-HELP] SAC not executing: " ERROR: Number 901 in INITCOMLISTS." Message-ID: <1263339734.5333.7.camel@kaven.wr.usgs.gov> Hi all! I'm having problems right after the installation. I received a tarball "sac-101.3b.tar.gz" and followed the instructions in the email as well as those listed in $SACHOME/README Problem: SAC isn't executing despite following the installation instructions carefully. System information: System: "uname -ris": Linux 2.6.18-164.9.1.el5 x86_64 I'm running RedHat Linux 5.4 (CentOS 5.4) SAC Location: "echo $SACHOME": /usr/local/bin Settings & ~/.bashrc: export SACHOME=/usr/local/sac export PATH=${PATH}:${SACHOME}/bin export SACAUX=${SACHOME}/aux export SAC_USE_DATABASE=0 export SAC_PPK_LARGE_CROSSHAIRS=0 When calling "./sac", this error is output: " /aux/csspickpref . readcss will not read picks. ERROR: Number 901 in INITCOMLISTS. /aux/clstdpening Command List: /usr/local/sac " This seem awefully similar to the problem reported in this thread: http://www.iris.washington.edu/pipermail/sac-help/2006-November/000188.html for which no solution was found and only a newer version helped. I have no clue what is going wrong. Any suggestions? Thanks in advance, Ole From magreent at gmail.com Wed Jan 13 09:18:39 2010 From: magreent at gmail.com (Mac Green) Date: Wed, 13 Jan 2010 11:18:39 -0600 Subject: [SAC-HELP] FFT with unevenly sampled data in SAC Message-ID: Hi all, I am using data collected with labview and saved as text files. I have to sets of data, the amplitude data set and a set of time intervals which are not even. I use readtable to import the data into sac and then write it as a sac file. If my data was evenly sampled, I would use delta in the header to give the time interval. Does anyone know how I can get the FFT of data whose delta is not constant, it I have a second file with time intervals? Thanks. -------------- next part -------------- An HTML attachment was scrubbed... URL: From ldorman at ucsd.edu Wed Jan 13 10:17:38 2010 From: ldorman at ucsd.edu (LeRoy Dorman) Date: Wed, 13 Jan 2010 10:17:38 -0800 Subject: [SAC-HELP] FFT with unevenly sampled data in SAC In-Reply-To: References: Message-ID: <4B4E0E42.5030608@ucsd.edu> This question has no simple answer, and is the subject of current research. As a starting point you might look at the following papers: I. M. Fodor and P. B. Stark, Multitaper Spectral Estimation for Time Series with Gaps, IEEE transactions on Signal Processing, Vol 48, No 12 December, 2000 and www.congrex.nl/09c24/S3_Posters/S3_P05_Smith_paper.pdf Hope this helps. LeRoy Dorman LeRoy M. Dorman Professor Scripps Institution of Oceanography, MC0220 University of California, San Diego 9500 Gilman Drive, La Jolla, CA 92093-0220 voice: +1 858 534-2406 fax: +1 858 534-6849 inet: ldorman at ucsd.edu http://www.mpl.ucsd.edu/people/ldorman/index.html Physical address: 436 Ritter Hall 8635 Charles Kennel Way (formerly Discovery Way) La Jolla, CA 92037 Mac Green wrote: > Hi all, > > I am using data collected with labview and saved as text files. I have > to sets of data, > the amplitude data set and a set of time intervals which are not even. > > I use readtable to import the data into sac and then write it as a sac > file. If my data > was evenly sampled, I would use delta in the header to give the time > interval. Does anyone > know how I can get the FFT of data whose delta is not constant, it I > have a second file with > time intervals? > > Thanks. From hectorg at miners.utep.edu Wed Jan 20 10:22:39 2010 From: hectorg at miners.utep.edu (hectorg at miners.utep.edu) Date: Wed, 20 Jan 2010 11:22:39 -0700 Subject: [SAC-HELP] SAC Zero Y-axis Message-ID: I would like to know how to draw a zero y-axis (in addition to the top and bottom axes) when using ppk. Thank you. Atte Hector Gz. Huizar -------------- next part -------------- An HTML attachment was scrubbed... URL: From N.Houlie at leeds.ac.uk Mon Jan 25 09:38:27 2010 From: N.Houlie at leeds.ac.uk (Nicolas Houlie) Date: Mon, 25 Jan 2010 17:38:27 +0000 Subject: [SAC-HELP] ppk? Message-ID: <147708352AE8E54BA8388E0BE2F6327A49BA5311C9@HERMES8.ds.leeds.ac.uk> >Hello, > >I am looking for some help to pick up times on SAC waveforms. I saw you >can use PPK but I do not find any help on your webpage (nothing in >graphics capabilities). > >is this possible to save the times/amplitudes in a file ? > >Thanks, > >nicolas Nicolas Houlie School of Earth and Environment 15-19 Hyde Terrace, Email: n.houlie at leeds.ac.uk From mplasencia at ogs.trieste.it Mon Jan 25 10:12:19 2010 From: mplasencia at ogs.trieste.it (Milton P. Plasencia Linares) Date: Mon, 25 Jan 2010 19:12:19 +0100 Subject: [SAC-HELP] ppk? In-Reply-To: <147708352AE8E54BA8388E0BE2F6327A49BA5311C9@HERMES8.ds.leeds.ac.uk> References: <147708352AE8E54BA8388E0BE2F6327A49BA5311C9@HERMES8.ds.leeds.ac.uk> Message-ID: <20100125191219.17222vb67lop8mwz@webmail.ogs.trieste.it> Hi, First, see the tutorial of Zhigang Peng (Other Basic tools - Phase Pick item) http://geophysics.eas.gatech.edu/classes/SAC/ Then, try this commands: Generate a file for the picks, use the command: touch picks sac % run SAC funcgen % generate impulse function, (here read your time series). oapf picks % Open the file where will be written the picks ppk % Pick the signal of your interest, use the commands of the Peng tutorial, x and left click for zoom, p and s, and then q on the window for quit of the ppk. capf picks % close the picks file. (important) In the file picks will be written the information of the picking including the amplitude. Read the man pages for PPK. (help ppk) You can extract the amplitude of the file picks using awk command. Cheers, ********************************** Milton P. PLASENCIA LINARES Dipartimento Centro di Ricerche Sismologiche Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS Borgo Grotta Gigante 42/C (34010) Sgonico - TRIESTE - ITALIA Tel: +39-040-2140305 Fax: +39-040-327307 E-mail: mplasencia at ogs.trieste.it ASAIN (Antarctic Seismographic Argentinean Italian Network) ********************************* Quoting "Nicolas Houlie" : >> Hello, >> >> I am looking for some help to pick up times on SAC waveforms. I saw you >> can use PPK but I do not find any help on your webpage (nothing in >> graphics capabilities). >> >> is this possible to save the times/amplitudes in a file ? >> >> Thanks, >> >> nicolas > > > Nicolas Houlie > School of Earth and Environment > 15-19 Hyde Terrace, > Email: n.houlie at leeds.ac.uk > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > > Email secured by **CeSIT** Check Point gateway > ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. From stubailo at gmail.com Tue Jan 26 12:42:26 2010 From: stubailo at gmail.com (Igor Stubailo) Date: Tue, 26 Jan 2010 20:42:26 +0000 Subject: [SAC-HELP] question about "do file wild *SAC" Message-ID: <4843440a1001261242l6521924ds8523d8fbf5dd170d@mail.gmail.com> Dear SAC users, Could you please help me to find out what's wrong with my file reading sequence? The code snippet from sac_LHZ.test: ... do file wild *SAC setbb vert $file r %vertdo file wild *SAC setbb station &1,KSTNM setbb otime &1,O evaluate to cotime %otime * (-1.0) ch ALLT %cotime IZTYPE IO setbb network &1,KNETWK setbb polezero %polezerod%%station%.%network%.PZ w %vert enddo ... The content of the current dir: PTRP.TO.SAC LVOF.TO.SAC MMMR.TO.SAC The errors: >do file wild *SAC >./tr: line 4: 26622 Segmentation fault sac <m sac_LHZ.test >quit >! The script I run: >cat tr sac < From stubailo at gmail.com Tue Jan 26 13:10:48 2010 From: stubailo at gmail.com (Igor Stubailo) Date: Tue, 26 Jan 2010 21:10:48 +0000 Subject: [SAC-HELP] question about "do file wild *SAC" In-Reply-To: <4B5F596C.8030500@lbl.gov> References: <4843440a1001261242l6521924ds8523d8fbf5dd170d@mail.gmail.com> <4B5F596C.8030500@lbl.gov> Message-ID: <4843440a1001261310m5cbd8f1cx4eca0ca463887a6@mail.gmail.com> Oops.. Thanks for pointing that out, that was copy and paste error caused on my side. This is the correct version: do file wild *SAC setbb vert $file r %vert setbb station &1,KSTNM setbb otime &1,O evaluate to cotime %otime * (-1.0) ch ALLT %cotime IZTYPE IO setbb network &1,KNETWK setbb polezero %polezerod%%station%.%network%.PZ w %vert enddo Also, the SAC version I use is 101.3b Sorry about that, Igor On Tue, Jan 26, 2010 at 9:06 PM, Katie Boyle wrote: > i'm not sure if it's causing the error, but you need to separate the read > command "r %vert" from the do loop "do file wild *SAC" > > right now they are on the same line... > > > " r %vertdo file wild *SAC" > > Igor Stubailo wrote: > >> Dear SAC users, >> >> Could you please help me to find out what's wrong with my file reading >> sequence? >> >> The code snippet from sac_LHZ.test: >> ... >> do file wild *SAC >> setbb vert $file >> r %vertdo file wild *SAC >> setbb station &1,KSTNM >> setbb otime &1,O >> evaluate to cotime %otime * (-1.0) >> ch ALLT %cotime IZTYPE IO >> setbb network &1,KNETWK >> setbb polezero %polezerod%%station%.%network%.PZ >> w %vert >> enddo >> ... >> >> The content of the current dir: >> PTRP.TO.SAC LVOF.TO.SAC MMMR.TO.SAC >> >> The errors: >> >do file wild *SAC >> >./tr: line 4: 26622 Segmentation fault sac <> >m sac_LHZ.test >> >quit >> >! >> >> The script I run: >> >cat tr >> sac <> m sac_LHZ.test >> quit >> ! >> >> >> If I comment out the lines "do file.." and "enddo" and read a single file >> like: >> setbb vert PTRP.TO.SAC >> the code continues to execute fine. >> So the problem is about using wild with *SAC. >> >> Thanks, >> Igor >> ------------------------------------------------------------------------ >> >> _______________________________________________ >> sac-help mailing list >> sac-help at iris.washington.edu >> http://www.iris.washington.edu/mailman/listinfo/sac-help >> >> > > > -- > Katie Boyle > Lawrence Berkeley National Laboratory > 510.486.4018 > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From stubailo at gmail.com Tue Jan 26 14:58:28 2010 From: stubailo at gmail.com (Igor Stubailo) Date: Tue, 26 Jan 2010 22:58:28 +0000 Subject: [SAC-HELP] question about "do file wild *SAC" In-Reply-To: <4843440a1001261444x54f732dbmd063123d1aeb71f9@mail.gmail.com> References: <4843440a1001261242l6521924ds8523d8fbf5dd170d@mail.gmail.com> <4843440a1001261333q1e614c4ag6997abe5689845d7@mail.gmail.com> <4B5F61D9.3010906@lbl.gov> <4843440a1001261351j41610bdhaba0669dd21950b5@mail.gmail.com> <4B5F647B.3020501@lbl.gov> <4843440a1001261403h2192fd75lb2ec5d0dfead5ad3@mail.gmail.com> <4B5F676F.10908@lbl.gov> <4843440a1001261409l6d95bc52neee3fb39509f7ef5@mail.gmail.com> <4B5F6CD3.3070206@lbl.gov> <4843440a1001261444x54f732dbmd063123d1aeb71f9@mail.gmail.com> Message-ID: <4843440a1001261458p47ad5e31q5425ad1456dde0b1@mail.gmail.com> This is an update to my original message about "do file wild *SAC" problem. Katie Boyle helped me to figure out what was happening. So, the problem is with the sac memory limit. The script works fine with a few files but when I increase the number of them to 150, sac fails with the segfault. Apparently, "do file wild *SAC" wants to fit all *SAC data into the buffer and cannot. The data files are 1 sps, 2 hour long LHZs. Each file is 29K (29436) long. My sac version is the latest and greatest ver. 101.3b. Is there a way to increase the total data size limit that can be used in sac? The workaround in my case is to read the sac files one by one through the loop in the tshell outside sac. Thanks, Igor -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Tue Jan 26 16:04:59 2010 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 26 Jan 2010 19:04:59 -0500 (EST) Subject: [SAC-HELP] question about "do file wild *SAC" In-Reply-To: <4843440a1001261458p47ad5e31q5425ad1456dde0b1@mail.gmail.com> References: <4843440a1001261242l6521924ds8523d8fbf5dd170d@mail.gmail.com> <4843440a1001261333q1e614c4ag6997abe5689845d7@mail.gmail.com> <4B5F61D9.3010906@lbl.gov> <4843440a1001261351j41610bdhaba0669dd21950b5@mail.gmail.com> <4B5F647B.3020501@lbl.gov> <4843440a1001261403h2192fd75lb2ec5d0dfead5ad3@mail.gmail.com> <4B5F676F.10908@lbl.gov> <4843440a1001261409l6d95bc52neee3fb39509f7ef5@mail.gmail.com> <4B5F6CD3.3070206@lbl.gov> <4843440a1001261444x54f732dbmd063123d1aeb71f9@mail.gmail.com> <4843440a1001261458p47ad5e31q5425ad1456dde0b1@mail.gmail.com> Message-ID: In the HISTORY file that accompanies the version of SAC you are using, it says ... - Known Bugs/deficiencies - SSS traveltimes command (endian issue) - do file wild segmentation fault with a large number of files This bug will be corrected in the next release. When I first found this bug some years ago, I found it was not the number of files but the number of cumulative characters in the filespecs. The way we have gotten around it for now is to break up the list into several pieces: do file wild [A-J]*bhz enddo do file wild [K-N]*bhz enddo I think this has been discussed before in sac-help exchanges. On Tue, 26 Jan 2010, Igor Stubailo wrote: > This is an update to my original message about "do file wild *SAC" problem. > Katie Boyle helped me to figure out what was happening. > > So, the problem is with the?sac memory limit. The? script works fine with a > few files but when I increase the number of them to 150, sac fails with the > segfault. Apparently, "do file wild *SAC" wants to fit all *SAC data into > the buffer and cannot. > > The data files are 1 sps, 2 hour long LHZs. Each file is 29K (29436) long. > > My sac version is the latest and greatest ver. 101.3b. > Is there a way to increase the total data size limit that can be used in > sac? > > The workaround in my case is to read the sac files one by one through the > loop in the tshell outside sac. > > Thanks, > Igor > > From stubailo at gmail.com Tue Jan 26 16:53:57 2010 From: stubailo at gmail.com (Igor Stubailo) Date: Wed, 27 Jan 2010 00:53:57 +0000 Subject: [SAC-HELP] question about "do file wild *SAC" In-Reply-To: References: <4843440a1001261242l6521924ds8523d8fbf5dd170d@mail.gmail.com> <4843440a1001261351j41610bdhaba0669dd21950b5@mail.gmail.com> <4B5F647B.3020501@lbl.gov> <4843440a1001261403h2192fd75lb2ec5d0dfead5ad3@mail.gmail.com> <4B5F676F.10908@lbl.gov> <4843440a1001261409l6d95bc52neee3fb39509f7ef5@mail.gmail.com> <4B5F6CD3.3070206@lbl.gov> <4843440a1001261444x54f732dbmd063123d1aeb71f9@mail.gmail.com> <4843440a1001261458p47ad5e31q5425ad1456dde0b1@mail.gmail.com> Message-ID: <4843440a1001261653k67722e1ehe146db1608bda583@mail.gmail.com> Arthur, Thanks for clarification and one more way to work around this problem. -Igor On Wed, Jan 27, 2010 at 12:04 AM, Arthur Snoke wrote: > In the HISTORY file that accompanies the version of SAC you are using, it > says ... > > > - Known Bugs/deficiencies > - SSS traveltimes command (endian issue) > - do file wild segmentation fault with a large number of files > > This bug will be corrected in the next release. > > When I first found this bug some years ago, I found it was not the number > of files but the number of cumulative characters in the filespecs. The way > we have gotten around it for now is to break up the list into several > pieces: > do file wild [A-J]*bhz > enddo > do file wild [K-N]*bhz > enddo > > I think this has been discussed before in sac-help exchanges. > > > On Tue, 26 Jan 2010, Igor Stubailo wrote: > > This is an update to my original message about "do file wild *SAC" >> problem. >> Katie Boyle helped me to figure out what was happening. >> >> So, the problem is with the sac memory limit. The script works fine with >> a >> few files but when I increase the number of them to 150, sac fails with >> the >> segfault. Apparently, "do file wild *SAC" wants to fit all *SAC data into >> the buffer and cannot. >> The data files are 1 sps, 2 hour long LHZs. Each file is 29K (29436) long. >> >> My sac version is the latest and greatest ver. 101.3b. >> Is there a way to increase the total data size limit that can be used in >> sac? >> >> The workaround in my case is to read the sac files one by one through the >> loop in the tshell outside sac. >> >> Thanks, >> Igor >> >> -------------- next part -------------- An HTML attachment was scrubbed... URL: From sheila at blacknest.gov.uk Fri Jan 29 06:40:53 2010 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Fri, 29 Jan 2010 14:40:53 +0000 Subject: [SAC-HELP] readGSE lats/longs Message-ID: <4B62F375.4020200@blacknest.gov.uk> Does the latest SAC "readgse" command pick up the station latitudes/longitudes/elevations from the STA2 header line of GSE2.0 files? If not, are there plans to incorporate this into "readgse" in upcoming versions? Regards, Sheila Peacock.