From nadav.wetzler at gmail.com Thu Jul 2 10:49:10 2009 From: nadav.wetzler at gmail.com (Nadav Wetzler) Date: Thu, 2 Jul 2009 20:49:10 +0300 Subject: [SAC-HELP] Problem in installing SAC v101.2 - Mac Message-ID: > Thank you for sending me the Seismic Analysis Code software, Mac > version v101.2 > > I read the attached installation txt file. However, I have a problem > with the installation. When I type "make install", I get: > Making install in libedit > make[2]: Nothing to be done for `install-exec-am'. > make[2]: Nothing to be done for `install-data-am'. > Making install in utils > test -z "/usr/local/sac/bin" || /Users/nadavwetzler/documents/ > university/phd/softwares/sac-101.2/config/install-sh -d "/usr/local/ > sac/bin" > mkdir: /usr/local/sac: Permission denied > make[2]: *** [install-binPROGRAMS] Error 1 > make[1]: *** [install-am] Error 2 > make: *** [install-recursive] Error 1 > > I set the Permission to "Read & Write" for the sac-101.2 folder. > > What am I missing?? > > I'm using a MacBook 2.4, with Mac OS X version 10.5.7 > > Thanks From pbragato at inogs.it Fri Jul 3 04:31:18 2009 From: pbragato at inogs.it (pier luigi bragato) Date: Fri, 03 Jul 2009 13:31:18 +0200 Subject: [SAC-HELP] floating point exception Message-ID: <4A4DEC06.8030403@inogs.it> I've the same problem for a 64 bit computer running Slackware. Here the description of the system: % uname -a Linux antelope-wks 2.6.16.29 #2 SMP Tue Oct 31 13:41:37 CET 2006 x86_64 x86_64 x86_64 GNU/Linux the involved libraries are: % ldd /opt/sac/bin/sac linux-gate.so.1 => (0xffffe000) libSM.so.6 => /usr/lib/libSM.so.6 (0xf7fa6000) libICE.so.6 => /usr/lib/libICE.so.6 (0xf7f8e000) libX11.so.6 => /usr/lib/libX11.so.6 (0xf7ec1000) libdl.so.2 => /lib/libdl.so.2 (0xf7ebd000) libm.so.6 => /lib/libm.so.6 (0xf7e9a000) libncurses.so.5 => /lib/libncurses.so.5 (0xf7e58000) libc.so.6 => /lib/libc.so.6 (0xf7d28000) /lib/ld-linux.so.2 (0xf7fce000) All the libraries exist and are 32 bit. Here it is their description using "file -L": /usr/lib/libSM.so.6: ELF 32-bit LSB shared object, Intel 80386, version1 (SYSV), not stripped /usr/lib/libICE.so.6: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped /usr/lib/libX11.so.6: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped /lib/libdl.so.2: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped /lib/libm.so.6: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped /lib/libncurses.so.5: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), stripped /lib/libc.so.6: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped /lib/ld-linux.so.2: ELF 32-bit LSB shared object, Intel 80386, version 1 (SYSV), not stripped Any suggestion ? Pier Luigi Bragato Email secured by **CeSIT** Check Point gateway From re.rezaei at gmail.com Sun Jul 12 05:36:41 2009 From: re.rezaei at gmail.com (reza rezaei) Date: Sun, 12 Jul 2009 16:06:41 +0330 Subject: [SAC-HELP] a question about transfer command in Sac Message-ID: <5feb258a0907120536v234d2c86gc1b47b734941bd5f@mail.gmail.com> Hi I have a problem with transfer command in Sac. I have removed an instrument response and applied Wood Anderson to a wave in two different ways but I saw some problems. First, I removed SS1 response with: transfer from polezero subtype $resp freq 1 10, then I applied WA response with this command: transfer to WA Furthermore, I did both with this command: transfer from polezero subtype $resp to WA freq 1 10 >From the above you might got idea that the result should be similar, but the results from the first and the second method are not similar. I have attached a file. Please run and look at the script in it. You can change parameters and test it. I am very grateful if you help me. Thanks, Reza -------------- next part -------------- A non-text attachment was scrubbed... Name: transfer.tar.gz Type: application/x-gzip Size: 22675 bytes Desc: not available URL: From snoke at vt.edu Sun Jul 12 13:15:32 2009 From: snoke at vt.edu (Arthur Snoke) Date: Sun, 12 Jul 2009 16:15:32 -0400 (EDT) Subject: [SAC-HELP] a question about transfer command in Sac In-Reply-To: <5feb258a0907120536v234d2c86gc1b47b734941bd5f@mail.gmail.com> References: <5feb258a0907120536v234d2c86gc1b47b734941bd5f@mail.gmail.com> Message-ID: I have been working on the transfer function, so I looked at your example. A few comments: It is not clear which version of sac you are using, as your command is sac2000. It probably does not matter. The command sequence rmean rtr could be replaced by rtr only. Try both ways and compare the lh output. The command FREQ (short for FREQLIMITS) should have 4 arguments. I am not sure what it means if you put in only two. I am guessing that you want the filter response flat between 1 and 10 Hz. Here is a run with settings I think do what you want. I did not do all the tests. Recall that a second call to transfer with no FREQ will probably use the values from the first time. Not sure that is what you want. SAC> r 2000-08-03-0620-19S.TBZ___001_SP_E_TBZ___SP_E_SAC SAC> rtr Slope and standard deviation are: -0.00089445 0.14512 Intercept and standard deviation are: 0.13525 22.458 Data standard deviation is: 1300.1 Data correlation coefficient is: 5.3244e-05 SAC> transfer from polezero subtype ss1rd3_2.pz to WA FREQ 0.1 1.0 10.0 20.0 Extracting polezero response for TBZ SP , SP E ... Station (TBZ SP ), Channel (SP E ) Waveform multiplied by 1.000000 after deconvolution. SAC> p1 SAC> ppk SAC> write wa.sac SAC> quit Let me know if you (aor anyone else) have questions about what I am doing or claiming, as I want the next writeup on transfer to be both complete and accurate. Arthur Snoke On Sun, 12 Jul 2009, reza rezaei wrote: > Hi > I have a problem with transfer command in Sac. > I have removed an instrument response and applied Wood Anderson to a > wave in two different ways but I saw some problems. > > First, I removed SS1 response with: transfer from polezero subtype > $resp freq 1 10, then I applied WA response with this command: > transfer to WA > Furthermore, I did both with this command: transfer from polezero > subtype $resp to WA freq 1 10 > >> From the above you might got idea that the result should be similar, > but the results from the first and the second method are not similar. > > I have attached a file. Please run and look at the script in it. You > can change parameters and test it. I am very grateful if you help me. > > Thanks, > Reza > From cquinteros at funvisis.gob.ve Sun Jul 12 19:26:50 2009 From: cquinteros at funvisis.gob.ve (CLAUDIA QUINTEROS) Date: Sun, 12 Jul 2009 23:26:50 -0300 Subject: [SAC-HELP] Sac files UNIX to Mac OS Message-ID: <20090713021148.M55210@funvisis.gob.ve> Hi, I used to use SAC in UNIX, but currently I have SAC for Mac OS. I would like to read the files in SAC format but I don't know how to byteswap files. Thanks. Claudia Quinteros. -- Claudia Quinteros Ingeniera Geof?sica FUNVISIS - Fundaci?n Venezolana de Investigaciones Sismol?gicas Tel?fonos: +58212 257-5153/7672/9346 Ext. 271 Fax: +58212 257-9977 Celular: +58412 737-6381 http://www.funvisis.gob.ve -- Este mensaje ha sido analizado por MailScanner en busca de virus y otros contenidos peligrosos, y se considera que est? limpio. Fundaci?n Venezolana de Investigaciones Sismol?gicas (FUNVISIS) From m.tchelo at gmail.com Mon Jul 13 05:49:23 2009 From: m.tchelo at gmail.com (Marcelo Bianchi) Date: Mon, 13 Jul 2009 09:49:23 -0300 Subject: [SAC-HELP] Sac files UNIX to Mac OS In-Reply-To: <20090713021148.M55210@funvisis.gob.ve> References: <20090713021148.M55210@funvisis.gob.ve> Message-ID: <89b338c10907130549m4ee37cb8y633a6360a9e4f277@mail.gmail.com> Claudia If you want I write a tool that do this job but I think that the new versions of the sac software could read both endianess. The software that I write can be obtained from: http://my.opera.com/m.tchelo/blog/2009/05/18/sacswapbyte-version-1-0 best regards, marcelo 2009/7/12 CLAUDIA QUINTEROS > > Hi, > > I used to use SAC in UNIX, but currently I have SAC for Mac OS. > I would like to read the files in SAC format but I don't know how to > byteswap > files. > > Thanks. > Claudia Quinteros. > > -- > > > Claudia Quinteros > Ingeniera Geof?sica > FUNVISIS - Fundaci?n Venezolana de Investigaciones Sismol?gicas > Tel?fonos: +58212 257-5153/7672/9346 Ext. 271 > Fax: +58212 257-9977 > Celular: +58412 737-6381 > http://www.funvisis.gob.ve > > > -- > Este mensaje ha sido analizado por MailScanner > en busca de virus y otros contenidos peligrosos, > y se considera que est? limpio. > Fundaci?n Venezolana de Investigaciones Sismol?gicas (FUNVISIS) > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > -- :: "This is all about life !" :: :: Find me @ http://my.opera.com/m.tchelo :: -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Mon Jul 13 07:21:52 2009 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 13 Jul 2009 10:21:52 -0400 (EDT) Subject: [SAC-HELP] Sac files UNIX to Mac OS In-Reply-To: <89b338c10907130549m4ee37cb8y633a6360a9e4f277@mail.gmail.com> References: <20090713021148.M55210@funvisis.gob.ve> <89b338c10907130549m4ee37cb8y633a6360a9e4f277@mail.gmail.com> Message-ID: The LLNL/IRIS distribution of SAC can read SAC files in either endian. In addition, the distribution includes a utility program sacswap that can switch endians. The next distribution of SAC will also include a program that can swap the byte order of SGF files (plot files produced by SAC). From snoke at vt.edu Mon Jul 13 14:56:28 2009 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 13 Jul 2009 17:56:28 -0400 (EDT) Subject: [SAC-HELP] a question about transfer command in Sac In-Reply-To: <5feb258a0907120536v234d2c86gc1b47b734941bd5f@mail.gmail.com> References: <5feb258a0907120536v234d2c86gc1b47b734941bd5f@mail.gmail.com> Message-ID: When running the SAC TRANSFER a second time within a SAC session, it remembers the settings for things like FROM and TO and FREQ from the first run. Hence in Reza's example, SAC> r 2000-08-03-0620-19S.TBZ___001_SP_E_TBZ___SP_E_SAC SAC> rtr SAC> transfer from polezero subtype ss1rd3_2.pz freq 0.1 1 10 12 followed by SAC> transfer to WA actually was translated by SAC as the complete command transfer from polezero subtype ss1rd3_2.pz to WA freq 0.1 1 10 12 which is what he did in his second run, but the waveform has been transformed by the first transform. One could either follow the first transfer with a SAC r SAC> transfer to WA or one could write out the waveform after the first transfer SAC> write try0.sac SAC> quit Restart SAC SAC> r try0.sac SAC> transfer to WA SAC> write try4.sac Shown in the attached plot is a plot of this try4 and try2, which was written applying the two-stage transfer to the original waveform. Arthur -------------- next part -------------- A non-text attachment was scrubbed... Name: try2-try4.pdf Type: application/pdf Size: 19318 bytes Desc: URL: From snoke at vt.edu Tue Jul 14 05:04:19 2009 From: snoke at vt.edu (Arthur Snoke) Date: Tue, 14 Jul 2009 08:04:19 -0400 (EDT) Subject: [SAC-HELP] clarifying calls to TRANSFER Message-ID: Inspired by the recent discussion about TRANSFER on this list, I am putting in the following paragraph in the TRANSFER HELP file: Care must be taken when calling TRANSFER a second time within a single SAC session, because in the second call TRANSFER will use the same arguments for FROM, TO, FREQ, etc. as in the first call unless an alternative argument is explicitly provided. Comments? From cquinteros at funvisis.gob.ve Tue Jul 14 10:40:43 2009 From: cquinteros at funvisis.gob.ve (Claudia Quinteros) Date: Tue, 14 Jul 2009 13:10:43 -0430 Subject: [SAC-HELP] Sac files UNIX to Mac OS Message-ID: <4A5CC31B.4000206@funvisis.gob.ve> Hi, Certainly, the LLNL/IRIS distribution of SAC can read SAC files in either endian. I was trying to read a bad data set and for that reason SAC couldn't recognize my SAC files. Now, I'm using good data and I don't have to byteswap the file to read it with SAC in Mac. Thank you very much.. -- ________________________________________________________________ Claudia Quinteros Ingeniera Geof?sica FUNVISIS - Fundaci?n Venezolana de Investigaciones Sismol?gicas /Tel?fonos: (212) 2575153 / 7672 / 9346. Ext. 271./ /Fax: (212) 2579977/ Celular: (0412) 7376381 / (0424) 1489930 /http://www.funvisis.gob.ve/ -- Este mensaje ha sido analizado por MailScanner en busca de virus y otros contenidos peligrosos, y se considera que est? limpio. Fundaci?n Venezolana de Investigaciones Sismol?gicas (FUNVISIS) -------------- next part -------------- An HTML attachment was scrubbed... URL: From stubailo at moho.ess.ucla.edu Tue Jul 14 14:30:46 2009 From: stubailo at moho.ess.ucla.edu (Igor Stubailo) Date: Tue, 14 Jul 2009 14:30:46 -0700 (PDT) Subject: [SAC-HELP] rotated plot question Message-ID: Dear SAC users, I am trying to get ps images in my script but all plots become rotated by -90 degrees at the end: SAC>r 2009010129/* SAC>p1 --> plots 6 stations where t=0 is on the LEFT SAC>SGF PREFIX data NUMBER 1 DIRECTORY /home/test/data SIZE NORMAL SAC>bd sgf; p1; ed sgf SAC>sgftops data001.sgf data001.ps >gv data001.ps --> plots 6 stations where t=0 is on the BOTTOM Why this rotation is happenig? Is this normal? How do I rotate it back? My version: SEISMIC ANALYSIS CODE [01/29/2008 (Version 101.1)] Regards, Igor From stubailo at moho.ess.ucla.edu Wed Jul 15 14:59:12 2009 From: stubailo at moho.ess.ucla.edu (Igor Stubailo) Date: Wed, 15 Jul 2009 14:59:12 -0700 (PDT) Subject: [SAC-HELP] number formatting question In-Reply-To: References: <20090713021148.M55210@funvisis.gob.ve> <89b338c10907130549m4ee37cb8y633a6360a9e4f277@mail.gmail.com> Message-ID: Dear SAC community, I would like to display a BAZ in the title of the plot but the number is somewhat hard to read and it takes too many characters to display: SAC> read /home/data/test.HHZ.sac SAC> ch evla 18 evlo 12 evdp 10 mag 8 SAC> setbb baz &1,BAZ SAC> getbb baz baz = 6.559583e+01 How do I display the number as 65.59 without going through the shell processing? How would you do from within the SAC with a shell help? Thanks, Igor From hmahardi at ictp.it Sat Jul 25 05:35:19 2009 From: hmahardi at ictp.it (Mahardika Harry) Date: Sat, 25 Jul 2009 14:35:19 +0200 (CEST) Subject: [SAC-HELP] How to convert sac to ascii format Message-ID: <4056.10.57.18.10.1248525319.squirrel@webmail.ictp.it> Greetings, Do anybody know how to convert a sac file into ascii format. I'm doing a thesis research in cross-correlation of seismic noise , and I want to use several inputs in ascii, while my noise data are in sac format. Can anyone help me? Harry Mahardika Earth System Physics Diploma Programme The Abdus Salam ICTP From rm438 at cam.ac.uk Sat Jul 25 06:28:47 2009 From: rm438 at cam.ac.uk (Bob Myhill) Date: Sat, 25 Jul 2009 14:28:47 +0100 Subject: [SAC-HELP] How to convert sac to ascii format In-Reply-To: <4056.10.57.18.10.1248525319.squirrel@webmail.ictp.it> References: <4056.10.57.18.10.1248525319.squirrel@webmail.ictp.it> Message-ID: Dear Harry, Use the CONVERT function (which can be run within a macro if necessary). Details can be found in the SAC Manual at http://www.iris.edu/manuals/sac/manual.html. Bob -- __________________________________ Bob Myhill Department of Earth Sciences University of Cambridge Bullard Laboratories Madingley Rise Cambridge CB3 0EZ United Kingdom Email: rm438 at cam.ac.uk __________________________________ -------------- next part -------------- An HTML attachment was scrubbed... URL: From jollya2k at bk.ru Sat Jul 25 11:04:50 2009 From: jollya2k at bk.ru (=?koi8-r?Q?=E1=CC=C5=CB=D3=C1=CE=C4=D2_=F6?=) Date: Sat, 25 Jul 2009 22:04:50 +0400 Subject: [SAC-HELP] =?koi8-r?b?SG93IHRvIGNvbnZlcnQgc2FjIHRvIGFzY2lpIGZv?= =?koi8-r?b?cm1hdA==?= In-Reply-To: <4056.10.57.18.10.1248525319.squirrel@webmail.ictp.it> References: <4056.10.57.18.10.1248525319.squirrel@webmail.ictp.it> Message-ID: -----Original Message----- From: "Mahardika Harry" To: sac-help at iris.washington.edu Date: Sat, 25 Jul 2009 14:35:19 +0200 (CEST) Subject: [SAC-HELP] How to convert sac to ascii format > Greetings, > > Do anybody know how to convert a sac file into ascii format. > > I'm doing a thesis research in cross-correlation of seismic noise , and I > want to use several inputs in ascii, while my noise data are in sac > format. > > Can anyone help me? > > Harry Mahardika > Earth System Physics > Diploma Programme > The Abdus Salam ICTP > > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > Dear Harry, Working with cross-corralation functions I used command WRITE ALPHA for getting ascii file. To open such ascii file in SAC it is possible with command READ ALPHA. Alexander S. Zholondz Saint-Petersburg State University Russia