From savage at uri.edu Thu Apr 1 07:28:02 2010 From: savage at uri.edu (Brian Savage) Date: Thu, 1 Apr 2010 10:28:02 -0400 Subject: [SAC-HELP] matlab In-Reply-To: References: Message-ID: <0A34657C-79FC-48AB-95C5-3E853A76D373@uri.edu> Dear Laura, Are you trying to use the SAC's Matlab interface or do you want to manipulate data in Matlab. If you would like to manipulate data in matlab there are routines to do this in (v. 101.2) sac/utils/getsacdata.m sac/utils/padcat.m sac/utils/readsac.m sac/utils/sacstruct.m sac/utils/writesac.m If you would like to use SAC's Matlab interface (running Matlab from in the SAC program) this is more difficult. First, make sure the SAC was compiled with Matlab support (not the default) and then make sure you have the libraries required (libeng and libmex). Brian savage at uri.edu On Mar 29, 2010, at 3:53 AM , Laura Parisi wrote: > Hi, > I need to use matlab within sac but I've some problems. > > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2- > linux86$ ./bin/sacinit.sh > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export SACHOME=/media/ACERDATA/software/sac101.2-linux86 > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2- > linux86$ export PATH=${PATH}:${SACHOME}/bin > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export SACAUX=${SACHOME}/aux > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export SAC_USE_DATABASE=0 > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export SAC_PPK_LARGE_CROSSHAIRS=1 > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export LIBRARY_PATH=$LIBRARY_PATH:/media/ACERDATA/software/ > matlab2008/bin/glnx86/libeng.so > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ > export LIBRARY1_PATH=$LIBRARY1_PATH:/media/ACERDATA/software/ > matlab2008/bin/glnx86/libmex.so > lauretta at lauretta-laptop:/media/ACERDATA/software/sac101.2-linux86$ / > bin/sac SEISMIC ANALYSIS CODE [03/07/2007 (Version 101.0)] > Copyright 1995 Regents of the University of California > > SAC> R /home/lauretta/20070113_1512A/20070113_151222.IV.AGST. > 00.SHZ.sac > SAC> PLOT > SAC> MAT > ERROR 8002: Cannot link to MATLAB shared object: > SAC> SETMAT /media/ACERDATA/software/matlab2008/bin/glnx86/MATLAB > SAC> MAT > ERROR 8002: Cannot link to MATLAB shared object: > > what's wrong? > > help me, please! > > > > Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. Sign > up now._______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= Thu Apr 1 11:43:29 2010 From: =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= (=?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=) Date: Thu, 1 Apr 2010 11:43:29 -0700 Subject: [SAC-HELP] Power spectra Message-ID: Hi, I am using MacSac 10.6f-grh104 and trying to compute a power spectra from 10 seconds synthetic seismograms. Is there a simple way of doing that in sac or do I have to take the fft and multiply it by it's conjugate? I can't get the PDS (in the SPE sub-process) command to work properly. I might be using bad widow length. Is there a rule of thumb as to what the length should be? Thanks From snoke at vt.edu Fri Apr 2 07:59:23 2010 From: snoke at vt.edu (Arthur Snoke) Date: Fri, 2 Apr 2010 10:59:23 -0400 (EDT) Subject: [SAC-HELP] traveltime in SAC Message-ID: If you have no interested in the sac > sss> traveltime command or in being able to add phase-arrival times to SAC files, you can stop reading now. Although TRAVELTIME is part of the SSS SAC Subprocess, it should be able to work on a single data file. I have not succeeded in getting SYNTAX 1 and SYNTAX 3 to work. (I am referring to terms used in the SAC> help traveltime command.) Unless one of you can point out errors in my tests below and/or easy fixes, nothing will be done in the soon to be released 101.4, but we hope it (and other calls in the SAC subprocess) can be corrected in the next release. Before I detail my tests with traveltime, I tell you about a program that does work very well for adding phase-arrival times to SAC waveforms: the taup toolkit taup_setsac command. The toolkit can be downloaded for most operating systems from . For my waveform file test.z (which has event distance and depth information added to the header and actually has observable phase arrivals for these phases): taup_setsac -evdpkm -ph P-1,PcP-2,pP-3,PP-5,sP-6 test.z produced the correct arrival times in test.z. This uses the iasp91 model, but could use any of a number of other models as well. At several institutuions, phase-arrival times are added using the iaspei-tau Fortran program ttimes or interpolation of JB tables. Program taup_setsac is much more flexible -- and powerful. Now for my tests with traveltime. For Syntax 1, on either a PPC Mac or Sun Solaris, using the commands SAC> r test.z SAC> sss Signal Stacking Subprocess. SAC/SSS> traveltime model iasp91 phase P pP PP sP PcP picks 1 SAC/SSS> quitsub SAC> lh gives the wrong answers for the TNMARKER set. This result is the same with the currant version or (on the Sun) for SEISMIC ANALYSIS CODE [07/21/2003 (Version 00.59.46)]. When the same sequence is run on an i686 Mac or Linux, the traveltime command leads to a segmentation fault. Part of the reason for failure for SYNTAX 1 is that it reads two binary files that are in sac/aux/tables. These files, iasp91.hed and iasp91.tbl, are written using the big-endian byte order.When these are replaced on the little-endian machines with the appropriate files, SAC bombs with a "bad sizes" message. For SYNTAX 3, I downloaded the taup package and ran taup_curve -h 611.0 -ph P,PP,pP,sP,PcP SAC> r test2.z SAC> sss Signal Stacking Subprocess. SAC/SSS> traveltime picks 1 taup taup_curve.gmt SAC/SSS> quitsub SAC> lh All the arrival times added to the header were at the same time -- incorrect for any of the arrivals. From santanubaruah2007 at rediffmail.com Thu Apr 8 01:46:41 2010 From: santanubaruah2007 at rediffmail.com (Santanu Baruah) Date: 8 Apr 2010 08:46:41 -0000 Subject: [SAC-HELP] =?utf-8?q?problemin_reading_sac_header?= Message-ID: <20100408084641.61059.qmail@f4mail202.rediffmail.com> Dear Sir, Please find a sac file. When I try to read its header information by typing "lh", the programme donot allow me to see the header information. It gives "ERROR 1301: No data files read in." , though I use the command "r" for reading the sac file before typing "lh". Interestingly the header information of the same file can be seen in an another system where Centos OS is there. In my case, I am using FEDORA 12. Is it a problem related to OS? Please kindly see the problem and let me know how to overcome this problem. I shall be highly obliged for your help. With Best Regards Santanu Baruah Santanu Baruah Geoscience Division North East Institute of Science and Technology (Formerly RRL) Jorhat-785006,Assam,India. Mobile: +91 9435514805(M) -------------- next part -------------- An HTML attachment was scrubbed... URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: 1993.192.13.45.59.ATD.sac.z Type: application/octet-stream Size: 144632 bytes Desc: not available URL: From toyokuni.genchi at nipr.ac.jp Mon Apr 12 04:00:18 2010 From: toyokuni.genchi at nipr.ac.jp (Genti Toyokuni) Date: Mon, 12 Apr 2010 20:00:18 +0900 Subject: [SAC-HELP] About definition of Fourier transform Message-ID: <20100412200018.d38a8d1d.toyokuni.genchi@nipr.ac.jp> Hello, I want to ask you definition of Fourier transform used in the fft command. Does the command perform the FFT with \exp(-i\omega t), or \exp(+i\omega t)? I tried to find it out from the SAC manual, but I could not. Regards, Genti Toyokuni From snoke at vt.edu Mon Apr 12 06:18:31 2010 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 12 Apr 2010 09:18:31 -0400 (EDT) Subject: [SAC-HELP] About definition of Fourier transform In-Reply-To: <20100412200018.d38a8d1d.toyokuni.genchi@nipr.ac.jp> References: <20100412200018.d38a8d1d.toyokuni.genchi@nipr.ac.jp> Message-ID: The forward transform uses your first option: exp(-i\omega t). I thought it was documented but did not see it with a quick check. if it is there somewhere, please let me know where. If not, we will put it in the next release -- probably in the help file for FFT, if no one has a better idea. On Mon, 12 Apr 2010, Genti Toyokuni wrote: > Hello, > > I want to ask you definition of Fourier transform used in the fft command. > Does the command perform the FFT with \exp(-i\omega t), or \exp(+i\omega t)? > > I tried to find it out from the SAC manual, but I could not. > > Regards, > > Genti Toyokuni > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > From snoke at vt.edu Mon Apr 12 17:09:51 2010 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 12 Apr 2010 20:09:51 -0400 (EDT) Subject: [SAC-HELP] About definition of Fourier transform In-Reply-To: References: <20100412200018.d38a8d1d.toyokuni.genchi@nipr.ac.jp> Message-ID: If the forward fft has a +i, the unwrapped phase will increase with frequency, for a -i (SAC convention) the phase will decrease. An easy way to check is the following sequence: SAC> fg seis SAC> unwrap DC level after DFT is -0.98547 SAC> plotsp Waiting SAC> The plot displayed after "Waiting" is the unwrapped phase, which clearly decreases with increasing frequency. I have not used unwrap, and the HELP file says "The maximum length of data that can currently be transformed is 4096." The help file was last updated in January 8, 1983, If that limitation has been relaxed, the unwrap command might be useful for some: "The data is filled with zeros to the next power of two before being transformed. You may specify a larger number of zeros by using the FILL option." My non-SAC program for working with time series pads time series by a factor of 4 before the FFT. On Mon, 12 Apr 2010, Arthur Snoke wrote: > The forward transform uses your first option: exp(-i\omega t). > > I thought it was documented but did not see it with a quick check. if it is > there somewhere, please let me know where. If not, we will put it in the > next release -- probably in the help file for FFT, if no one has a better > idea. > > On Mon, 12 Apr 2010, Genti Toyokuni wrote: > >> Hello, >> >> I want to ask you definition of Fourier transform used in the fft command. >> Does the command perform the FFT with \exp(-i\omega t), or \exp(+i\omega >> t)? >> >> I tried to find it out from the SAC manual, but I could not. >> >> Regards, >> >> Genti Toyokuni >> _______________________________________________ >> sac-help mailing list >> sac-help at iris.washington.edu >> http://www.iris.washington.edu/mailman/listinfo/sac-help >> > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > From kboyle at lbl.gov Tue Apr 13 14:08:35 2010 From: kboyle at lbl.gov (Katie Boyle) Date: Tue, 13 Apr 2010 14:08:35 -0700 Subject: [SAC-HELP] SAC 64-bit Installation Message-ID: Has anyone been able to successfully compile the SAC source code on a 64-bit system? I've got both versions of libncurses, and my X11 libraries are present. It seems like the problem is with linking. I'm stuck! When I use gcc -m32, I get: dfr.c:(.text+0x81c): undefined reference to `cmplxdiv' dfr.c:(.text+0x888): undefined reference to `cmplxcj' dfr.c:(.text+0x8a5): undefined reference to `cmplxmul' dfr.c:(.text+0x8af): undefined reference to `cmplxtof' dfr.c:(.text+0x905): undefined reference to `cmplxtof' dfr.c:(.text+0x91c): undefined reference to `aimag' /usr/bin/ld: final link failed: Invalid operation collect2: ld returned 1 exit status make[1]: *** [SAC] Error 1 make[1]: Leaving directory `/home/hutchings/Desktop/SYNTHESIS/PROGRAMS/sac/SAC/src' make: *** [all] Error 2 When I use gcc -m64, I get: dfr.c:(.text+0x7b6): undefined reference to `cmplxdiv' dfr.c:(.text+0x7c0): undefined reference to `cmplxtof' dfr.c:(.text+0x7ff): undefined reference to `cmplxmul' dfr.c:(.text+0x81c): undefined reference to `cmplxdiv' dfr.c:(.text+0x888): undefined reference to `cmplxcj' dfr.c:(.text+0x8a5): undefined reference to `cmplxmul' dfr.c:(.text+0x8af): undefined reference to `cmplxtof' dfr.c:(.text+0x905): undefined reference to `cmplxtof' dfr.c:(.text+0x91c): undefined reference to `aimag' collect2: ld returned 1 exit status make[1]: *** [SAC] Error 1 make[1]: Leaving directory `/home/hutchings/Desktop/SYNTHESIS/PROGRAMS/sac/SAC/src' make: *** [all] Error 2 Thanks! Katie -------------- next part -------------- An HTML attachment was scrubbed... URL: From giantz at telkom.net Wed Apr 14 05:58:01 2010 From: giantz at telkom.net (Gian) Date: Wed, 14 Apr 2010 19:58:01 +0700 Subject: [SAC-HELP] ERROR: Number 901 in INITCOMLISTS Message-ID: Hi All.. can someone help me with this problem? I can not execute the sac program.. I'm using centos 5 desktop 32 bit.. here's what I get : [ctuser at localhost ~]$ echo $SACHOME /home/ctuser/sac [ctuser at localhost ~]$ echo $SACAUX /home/ctuser/sac/aux [ctuser at localhost ~]$ echo $PATH /usr/kerberos/bin:/usr/local/bin:/usr/bin:/bin:/usr/X11R6/bin:/home/ctuser/bin:/home/ctuser/sac/bin [ctuser at localhost bin]$ ./sac /csspickpref . readcss will not read picks. ERROR: Number 901 in INITCOMLISTS. /clstdms opening Command List: /home/ctuser/sac/aux [ctuser at localhost bin]$ it's the same in shell promt.. Thx before... ============================================================================================= Segera nikmati Free Trial 60 hari Protector Postpaid layanan keamanan online bagiPelanggan Speedy. Info lebih lanjut hubungi 147 atau http://protector.telkomspeedy.com ============================================================================================= From savage at uri.edu Wed Apr 14 07:38:59 2010 From: savage at uri.edu (Brian Savage) Date: Wed, 14 Apr 2010 10:38:59 -0400 Subject: [SAC-HELP] ERROR: Number 901 in INITCOMLISTS In-Reply-To: References: Message-ID: The messages you are getting from SAC look like the SACAUX environment variable is not set correctly. If it were, your should see something like /home/ctuser/sac/aux/csspickprefs and /home/ctuser/sac/aux/clstd You may have a return or extra character in the environment variable. You can check this by doing % echo "- $SACAUX -" Brian Savage savage at uri.edu On Apr 14, 2010, at 8:58 AM , Gian wrote: > Hi All.. > can someone help me with this problem? > I can not execute the sac program.. > I'm using centos 5 desktop 32 bit.. > > here's what I get : > > [ctuser at localhost ~]$ echo $SACHOME > /home/ctuser/sac > > [ctuser at localhost ~]$ echo $SACAUX > /home/ctuser/sac/aux > > [ctuser at localhost ~]$ echo $PATH > /usr/kerberos/bin:/usr/local/bin:/usr/bin:/bin:/usr/X11R6/bin:/home/ > ctuser/bin:/home/ctuser/sac/bin > > [ctuser at localhost bin]$ ./sac > /csspickpref . readcss will not read picks. > ERROR: Number 901 in INITCOMLISTS. > /clstdms opening Command List: /home/ctuser/sac/aux > [ctuser at localhost bin]$ > > it's the same in shell promt.. > > Thx before... > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > ====================================================================== > > Segera nikmati Free Trial 60 hari Protector Postpaid layanan > keamanan online bagiPelanggan Speedy. > Info lebih lanjut hubungi 147 atau http://protector.telkomspeedy.com > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > ====================================================================== > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From pku132 at 163.com Thu Apr 15 03:34:08 2010 From: pku132 at 163.com (Lee) Date: Thu, 15 Apr 2010 18:34:8 +0800 Subject: [SAC-HELP] Who has the zero-polar file for Guralp CMG-3ESPC? Message-ID: <4BC6EB88.02B343.30629@m50-133.163.com> Hi, all, does anyone know the zero-polar points for Guralp CMG-3ESPC? Thx. Lee 2010-04-15 From sheila at blacknest.gov.uk Thu Apr 15 03:57:25 2010 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Thu, 15 Apr 2010 11:57:25 +0100 Subject: [SAC-HELP] Who has the zero-polar file for Guralp CMG-3ESPC? In-Reply-To: <4BC6EB88.02B343.30629@m50-133.163.com> References: <4BC6EB88.02B343.30629@m50-133.163.com> Message-ID: <4BC6F115.7070007@blacknest.gov.uk> You can get the data sheets for your particular Guralp seismometer by emailing "caldoc at guralp.com" with the serial number of the seismometer in the subject line. It also works for digitisers, but for both, only if they were recently manufactured. The reply consists of Microsoft Word documents. Beware that these DO NOT OPEN PROPERLY in OpenOffice - always open them in a genuine copy of Microsoft Word. also beware that the poles and zeros are for velocity in Hz not displacement in radians, and have to be converted for SAC. Sheila Peacock Blacknest From dkementz at auth.gr Thu Apr 15 04:43:35 2010 From: dkementz at auth.gr (Despina Kementzetzidou) Date: Thu, 15 Apr 2010 14:43:35 +0300 Subject: [SAC-HELP] Who has the zero-polar file for Guralp CMG-3ESPC? In-Reply-To: <4BC6EB88.02B343.30629@m50-133.163.com> References: <4BC6EB88.02B343.30629@m50-133.163.com> Message-ID: <4BC6FBE7.6060900@auth.gr> For the model of 60sec I have the following prices (all in rads, velocity)/ Kony SENSOR Model SENSOR SENSITIVITY (V/m/sec) A0 (velocity) DIGITIZER SENSITIVITY (cnts/V) SAC CONSTANT (displacement, add 1 more zero) ZEROS POLES ZERO 1r ZERO 1im ZERO 2r ZERO 2im POLE 1r POLE 1im POLE 2r POLE 2im POLE 3r POLE 3im POLE 4r POLE 4im POLE 5r POLE 5im CMG-3ESPC/60sec 1950 5.71508E+08 400000 4.46E+17 2 5 0 0 0 0 -502.655 0 -1005.31 0 -1130.97 0 -0.07402 0.074016 -0.07402 -0.07402 Lee wrote: > Hi, all, > > does anyone know the zero-polar points for Guralp CMG-3ESPC? > > Thx. > > Lee > 2010-04-15 > > > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > > -------------- next part -------------- An HTML attachment was scrubbed... URL: From Onur.Tan at mam.gov.tr Thu Apr 15 04:10:50 2010 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Thu, 15 Apr 2010 14:10:50 +0300 Subject: [SAC-HELP] Who has the zero-polar file for Guralp CMG-3ESPC? In-Reply-To: <4BC6EB88.02B343.30629@m50-133.163.com> Message-ID: Dear Lee, I attached a 3ESPC calibration document from Guralp. Generaly they have the same Pole-Zero values. I recommend that you should send a mail to caldoc at guralp.com. write the sensor and digitizer serial numbers (i.e. T3W49, A203, T3W50, B201) in to the subject part. The system sends back the doc files. If you can not get files write to support team at Guralp. They will send you manually. The main point is that you must convert the poles-zeros in the document in to radian (poles*2*pi*f) for SAC, and add 1 extra zero for integration. and convert A0 value to radian and displacement (see the pdf file) type. best regards onur Dr. Onur TAN ---------------------------------------------- 40.7866N 29.4500E --------- T?B?TAK Marmara Ara?t?rma Merkezi, Yer ve Deniz Bilimleri Enstit?s? TUBITAK Marmara Research Center, Earth and Marine Sciences Institute Gebze - Kocaeli - TURKEY On 15.04.2010 13:34, "Lee" wrote: Hi, all, does anyone know the zero-polar points for Guralp CMG-3ESPC? Thx. Lee 2010-04-15 _______________________________________________ sac-help mailing list sac-help at iris.washington.edu http://www.iris.washington.edu/mailman/listinfo/sac-help -------------- next part -------------- A non-text attachment was scrubbed... Name: T3W49.doc Type: application/octet-stream Size: 55296 bytes Desc: T3W49.doc URL: -------------- next part -------------- A non-text attachment was scrubbed... Name: seisuk_data_manage.pdf Type: application/octet-stream Size: 572850 bytes Desc: seisuk_data_manage.pdf URL: From sheila at blacknest.gov.uk Thu Apr 15 08:33:19 2010 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Thu, 15 Apr 2010 16:33:19 +0100 Subject: [SAC-HELP] Who has the zero-polar file for Guralp CMG-3ESPC? In-Reply-To: References: Message-ID: <4BC731BF.3070604@blacknest.gov.uk> The 3ESPC calibration documents that I have seen are specified in velocity so the instructions below are probably correct. Read the calibration document from Guralp carefully: some Guralp seismometers have poles and zeros specified in acceleration, so you have to add TWO zeros to integrate to displacement for SAC. Check that your calibration document says "velocity response output" above the list of poles. Sheila P. Onur Tan wrote: > Dear Lee, > > I attached a 3ESPC calibration document from Guralp. Generaly they have the same Pole-Zero values. I recommend that you should send a mail to > caldoc at guralp.com. write the sensor and digitizer serial numbers (i.e. T3W49, A203, T3W50, B201) in to the subject part. The system sends back the doc files. If you can not get files write to support team at Guralp. They will send you manually. > > The main point is that you must convert the poles-zeros in the document in to radian (poles*2*pi*f) for SAC, and add 1 extra zero for integration. and convert A0 value to radian and displacement (see the pdf file) type. > > best regards > onur > > Dr. Onur TAN > ---------------------------------------------- 40.7866N 29.4500E --------- > T?B?TAK Marmara Ara?t?rma Merkezi, Yer ve Deniz Bilimleri Enstit?s? > TUBITAK Marmara Research Center, Earth and Marine Sciences Institute > Gebze - Kocaeli - TURKEY > > > > On 15.04.2010 13:34, "Lee" wrote: > > Hi, all, > > does anyone know the zero-polar points for Guralp CMG-3ESPC? > > Thx. > > Lee > 2010-04-15 > > > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help > > > > > ------------------------------------------------------------------------ > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From mplasencia at ogs.trieste.it Thu Apr 15 12:32:23 2010 From: mplasencia at ogs.trieste.it (Milton P. Plasencia Linares) Date: Thu, 15 Apr 2010 21:32:23 +0200 Subject: [SAC-HELP] About definition of Fourier transform Message-ID: <20100415213223.187348jmfi2cvzk7@webmail.ogs.trieste.it> Hi, I follow the commands of Artur Snoke and i obtain: SEISMIC ANALYSIS CODE [08/25/2009 (Version 101.3b)] Copyright 1995 Regents of the University of California SAC> fg seis SAC> unwrap remove this function: src/co/zmemad() The same results occur in Linux box and MacOSX, 64 bits sac binary version 101.3b. Any suggestions? Milton ********************************** 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-2140136 Fax: +39-040-327307 E-mail: mplasencia at ogs.trieste.it ASAIN (Antarctic Seismographic Argentinean Italian Network) ********************************* ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. From snoke at vt.edu Thu Apr 15 13:41:33 2010 From: snoke at vt.edu (Arthur Snoke) Date: Thu, 15 Apr 2010 16:41:33 -0400 (EDT) Subject: [SAC-HELP] About definition of Fourier transform In-Reply-To: <20100415213223.187348jmfi2cvzk7@webmail.ogs.trieste.it> References: <20100415213223.187348jmfi2cvzk7@webmail.ogs.trieste.it> Message-ID: I reproduced the results you got below for sac version 101.3b. However, it works in SEISMIC ANALYSIS CODE [09/04/2008 (Version 101.2)]. This does not help you in the immediate future, but "unwrap" works (again) in 101.4 On Thu, 15 Apr 2010, Milton P. Plasencia Linares wrote: > > Hi, > I follow the commands of Artur Snoke and i obtain: > > SEISMIC ANALYSIS CODE [08/25/2009 (Version 101.3b)] > Copyright 1995 Regents of the University of California > > SAC> fg seis > SAC> unwrap > remove this function: src/co/zmemad() > > The same results occur in Linux box and MacOSX, 64 bits > sac binary version 101.3b. > > Any suggestions? > > Milton > > ********************************** > 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-2140136 > Fax: +39-040-327307 > > E-mail: mplasencia at ogs.trieste.it From snoke at vt.edu Thu Apr 15 14:23:30 2010 From: snoke at vt.edu (Arthur Snoke) Date: Thu, 15 Apr 2010 17:23:30 -0400 (EDT) Subject: [SAC-HELP] running an older version of SAC Message-ID: If you have kept the SAC tree (or a binary tarball) for an older version of SAC that you may want to run, the following script may be useful: cat /usr/local/jas/bin/newsac #!/bin/tcsh setenv SACHOME "$1" setenv SACAUX "${SACHOME}/aux" setenv PATH "${SACHOME}/bin:${PATH}" ${SACHOME}/bin/sac /usr/local/jas/macros/init.m To test if the unwrap command worked in 101.2, I did the following: newsac /usr/local/jas/sac_archive/sac-080904 SEISMIC ANALYSIS CODE [09/04/2008 (Version 101.2)] Copyright 1995 Regents of the University of California SAC> fg seis SAC> unwrap DC level after DFT is -0.98547 SAC> quit From magreent at gmail.com Sun Apr 18 20:04:12 2010 From: magreent at gmail.com (Mac Green) Date: Sun, 18 Apr 2010 21:04:12 -0600 Subject: [SAC-HELP] "Cut" command in SAC Message-ID: How exactly do I implement the CUT command in SAC for a case like this: I have a file called data.sac The number of points=N I want only the first M points of the file. It is the option given in the help files that Uses B and N, but I can not implement it. Thanks in advance. -------------- next part -------------- An HTML attachment was scrubbed... URL: From makra at itsak.gr Mon Apr 19 00:10:11 2010 From: makra at itsak.gr (Konstantia Makra) Date: Mon, 19 Apr 2010 10:10:11 +0300 Subject: [SAC-HELP] "Cut" command in SAC In-Reply-To: References: Message-ID: <903EF9F5-4E51-43DE-BB81-450E9E94863E@itsak.gr> CUT command is used as follows cut t1 t2 r data.sac w over (if you want to keep the same filename) cut off t1 =N*dt t2 = M*dt To see the cut file r data.sac p1 Konstantia MAKRA Dr Civil Engineer Soil Dynamics Division Institute of Engineering Seismology & Earthquake Engineering (ITSAK) e-mail: makra at itsak.gr On 19 ??? 2010, at 6:04 ??, Mac Green wrote: > How exactly do I implement the CUT command in SAC for a case like > this: > > I have a file called data.sac > The number of points=N > I want only the first M points of the file. > > It is the option given in the help files that Uses B and N, but I > can not implement it. > > Thanks in advance. > > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help From makra at itsak.gr Mon Apr 19 00:12:14 2010 From: makra at itsak.gr (Konstantia Makra) Date: Mon, 19 Apr 2010 10:12:14 +0300 Subject: [SAC-HELP] "Cut" command in SAC In-Reply-To: <903EF9F5-4E51-43DE-BB81-450E9E94863E@itsak.gr> References: <903EF9F5-4E51-43DE-BB81-450E9E94863E@itsak.gr> Message-ID: <9AB4B404-6401-4E42-9E26-FFE0023B2486@itsak.gr> sorry, N is not the number of points in data.sac file but the first point you want to cut the file and M the last one. Konstantia MAKRA Dr Civil Engineer Soil Dynamics Division Institute of Engineering Seismology & Earthquake Engineering (ITSAK) e-mail: makra at itsak.gr On 19 ??? 2010, at 10:10 ??, Konstantia Makra wrote: > > CUT command is used as follows > > cut t1 t2 > r data.sac > w over (if you want to keep the same filename) > cut off > > t1 =N*dt > t2 = M*dt > > To see the cut file > > r data.sac > p1 > > > > > > Konstantia MAKRA > Dr Civil Engineer > Soil Dynamics Division > Institute of Engineering Seismology & Earthquake Engineering (ITSAK) > e-mail: makra at itsak.gr > > > > > On 19 ?????? 2010, at 6:04 ? ??, Mac Green wrote: > >> How exactly do I implement the CUT command in SAC for a case like >> this: >> >> I have a file called data.sac >> The number of points=N >> I want only the first M points of the file. >> >> It is the option given in the help files that Uses B and N, but I >> can not implement it. >> >> Thanks in advance. >> >> _______________________________________________ >> sac-help mailing list >> sac-help at iris.washington.edu >> http://www.iris.washington.edu/mailman/listinfo/sac-help > From Onur.Tan at mam.gov.tr Mon Apr 26 05:05:20 2010 From: Onur.Tan at mam.gov.tr (Onur Tan) Date: Mon, 26 Apr 2010 15:05:20 +0300 Subject: [SAC-HELP] SAC to SEG2 Message-ID: Hello, How can I convert SAC file to SEG2 format ? best regards Onur for Brian: The manual of the ?write? command on the web and its help file in SAC are different. There is no information about SEGY format in the web page. However their revision date are the same. Dr. Onur TAN ---------------------------------------------- 40.7866N 29.4500E --------- T?B?TAK Marmara Ara?t?rma Merkezi, Yer ve Deniz Bilimleri Enstit?s? TUBITAK Marmara Research Center, Earth and Marine Sciences Institute Gebze - Kocaeli - TURKEY