From panj151 at gmail.com Thu Nov 3 12:27:21 2011 From: panj151 at gmail.com (Jiatie Pan) Date: Thu, 3 Nov 2011 15:27:21 -0400 Subject: [SAC-HELP] ERROR: Number 901 in INITCOMLISTS , Problems opening Command List: /opt/sac/aux/clstd In-Reply-To: References: Message-ID: *Hi All..can someone help me with this problem?I can not execute the sac program.**I'm using Fedora12 64 bit. My installation as following:** [/home/pan/sac-101.4]: configure --prefix=/home/pan/sac ..... ..... Compliation Options for sac (Seismic Analysis Code 101.4) CFLAGS = LDFLAGS = LIBS = -ldl -lm -lncurses Debug = Optimization = Database = off Matlab = off X11-buffer = on X11-depth = on X_LIBS = X_CFLAGS = Readline = off Editline = on Install Base = /home/pan/sac System = linux_x86_64 (x86_64-unknown-linux-gnu) To compile sac type make [/home/pan/sac-101.4]:make ..... ..... [/home/pan/sac-101.4]:make install ..... ..... * *Then I set the environment variables just like what README file says. During this , no error messages founded. But when I execute **sac ,** [/home/pan]: sac WARNING: Number 108 /opt/sac/aux/csspickpref . readcss will not read picks. ERROR: Number 901 in INITCOMLISTS. Problems opening Command List: /opt/sac/aux/clstd [/home/pan]: which sac ~/sac/bin/sac [/home/pan]: echo $SACHOME /home/pan/sac [/home/pan]: echo "- $SACAUX -" - /home/pan/sac/aux - [/home/pan]: vi /home/pan/sac/bin/sacinit.csh ..... setenv SACHOME /home/pan/sac setenv PATH ${PATH}:${SACHOME}/bin setenv SACAUX ${SACHOME}/aux ...... **So I think the **environment variables **setting has no problem, but I still can not execute **sac**, I need your help and suggestions.Thanks!** Jiatie Pan * -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Fri Nov 4 03:29:51 2011 From: snoke at vt.edu (Arthur Snoke) Date: Fri, 4 Nov 2011 06:29:51 -0400 (EDT) Subject: [SAC-HELP] ERROR: Number 901 in INITCOMLISTS , Problems opening Command List: /opt/sac/aux/clstd In-Reply-To: References: Message-ID: Do a ln /home/pan/sac/ I had an error message like that once, and it turned out tht for some reason my sac/aux had the wrong name: sac/winaux. mv winaux to aux fixed it. On Thu, 3 Nov 2011, Jiatie Pan wrote: > > Hi All.. > can someone help me with this problem? > I can not execute the sac program. > I'm using Fedora12 64 bit. > My installation as following: > > > [/home/pan/sac-101.4]: configure --prefix=/home/pan/sac > ..... > ..... > Compliation Options for sac (Seismic Analysis Code 101.4) > CFLAGS = > LDFLAGS = > LIBS = -ldl -lm -lncurses > > > Debug = > Optimization = > Database = off > Matlab = off > X11-buffer = on > X11-depth = on > X_LIBS = > X_CFLAGS = > Readline = off > Editline = on > > > Install Base = /home/pan/sac > System = linux_x86_64 (x86_64-unknown-linux-gnu) > To compile sac type make > [/home/pan/sac-101.4]:make > ..... > ..... > [/home/pan/sac-101.4]:make install > > > ..... > ..... > Then I set the environment variables just like what README file says. > During this , no error messages founded. But when I execute sac , > > [/home/pan]: sac > WARNING: Number 108 /opt/sac/aux/csspickpref . readcss will not read picks. > ERROR: Number 901 in INITCOMLISTS. > Problems opening Command List: /opt/sac/aux/clstd > [/home/pan]: which sac > ~/sac/bin/sac > [/home/pan]: echo $SACHOME > /home/pan/sac > [/home/pan]: echo "- $SACAUX -" > - /home/pan/sac/aux - > [/home/pan]: vi /home/pan/sac/bin/sacinit.csh > ..... > setenv SACHOME /home/pan/sac > > > setenv PATH ${PATH}:${SACHOME}/bin > setenv SACAUX ${SACHOME}/aux > ...... > So I think the environment variables setting has no problem, but I still can not execute sac, I need your he > lp and suggestions.Thanks! > Jiatie Pan > > > From hlehto at mail.usf.edu Fri Nov 4 11:46:19 2011 From: hlehto at mail.usf.edu (Heather Lehto) Date: Fri, 4 Nov 2011 14:46:19 -0400 Subject: [SAC-HELP] Plotting original and synthetic waveforms for comparison In-Reply-To: References: Message-ID: I have a series of synthetic waveforms that I need to compare to the original waveform and chose the best fitting synthetic. Is there are way to plot each synthetic overlayed by the original waveform in a series of plots (one on top of the other) on a single page? -Heather -------------- next part -------------- An HTML attachment was scrubbed... URL: From mplasencia at ogs.trieste.it Fri Nov 4 12:30:51 2011 From: mplasencia at ogs.trieste.it (Milton P. Plasencia Linares) Date: Fri, 04 Nov 2011 20:30:51 +0100 Subject: [SAC-HELP] Plotting original and synthetic waveforms for comparison In-Reply-To: References: Message-ID: <20111104203051.196320btkp3oemi3@webmail.ogs.trieste.it> Hi Heather, You can follow the next commands: > read file1.sac % read the first file > color red increment % read the manual for customize, "help color" > read more file2.sac % read the second file > p2 % plot2, read "help plot2", plot a file over > the other one HTH, Cheers, 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) ********************************* Quoting "Heather Lehto" : > I have a series of synthetic waveforms that I need to compare to the > original waveform and chose the best fitting synthetic. Is there are way > to plot each synthetic overlayed by the original waveform in a series of > plots (one on top of the other) on a single page? > > -Heather > ---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program. From hrakotoarisaona at yahoo.com Sat Nov 5 14:47:50 2011 From: hrakotoarisaona at yahoo.com (Rakotoarisaona Hasiniaina) Date: Sat, 5 Nov 2011 14:47:50 -0700 (PDT) Subject: [SAC-HELP] (no subject) Message-ID: <1320529670.36652.yint-ygo-j2me@web59806.mail.ac4.yahoo.com> http://occalgym.com/wp-content/photos/misc/reg/markehjf.htm From shenxzh at gmail.com Wed Nov 9 20:15:00 2011 From: shenxzh at gmail.com (shenxzh at gmail.com) Date: Thu, 10 Nov 2011 12:15:00 +0800 Subject: [SAC-HELP] select waveform with sac Message-ID: <201111101214526857506@gmail.com> Hello everyone, I met a question to select waveform with sac. I hope I can check the waveform and save some good data in a folder or put bad data in another folder. The following is the part of my shell script. for fil in *.SAC do sac < From george.helffrich at bristol.ac.uk Wed Nov 9 21:30:57 2011 From: george.helffrich at bristol.ac.uk (George Helffrich) Date: Thu, 10 Nov 2011 05:30:57 +0000 Subject: [SAC-HELP] select waveform with sac In-Reply-To: <201111101214526857506@gmail.com> References: <201111101214526857506@gmail.com> Message-ID: Dear Xuzhang - SAC is reading input from the shell script, not your terminal. It is probably giving an end-of-file to the REPLY function. Use a review procedure structured like this: # Define sac macro called /tmp/view.m cat << 'EOF' > /tmp/view.m message "Processing $1$" read $1$ setbb resp "(reply Enter t to trash the file)" if _%resp% eq _t sc mv $1$ Trash endif 'EOF' # Create sac input file ls *.SAC | awk '{print "m /tmp/view.m",$1}' > /tmp/sac.in # Run SAC once, process all files sac qdp off ygrid on message "Many questions about files will happen now ..." m /tmp/sac.in quit On 10 Nov 2011, at 04:15, shenxzh at gmail.com wrote: > Hello everyone, > I met a question to select waveform with sac. I hope I can check the waveform and save some good data in a folder or put bad data in another folder. The following is the part of my shell script. > for fil in *.SAC > do > sac < qdp off > ygrid on > r AA.SAC > p1 > setbb resp (REPLY "Enter t to trash the file ") > if %resp eq "t" then > sc mv AA.SAC Trash > endif > q > ENDD > done > After excuting the shell script, there was no any pause and just put the *.SAC into the Trash folder. What' s problem in my script and how to solve it. > Thank you. > Cheers, Xuzhang > > 2011-11-10 > _______________________________________________ > sac-help mailing list > sac-help at iris.washington.edu > http://www.iris.washington.edu/mailman/listinfo/sac-help George Helffrich george.helffrich at bris.ac.uk From u2injesus at gmail.com Thu Nov 17 09:20:25 2011 From: u2injesus at gmail.com (jenni wang) Date: Thu, 17 Nov 2011 18:20:25 +0100 Subject: [SAC-HELP] error on spectrogram Message-ID: I am using SAC version 101.4 on Ubuntu latest version. I was running the SPECTROGRAM command when I came to the error: Window size: 200 Overlap: 100 FFT size: 512 Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-04. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-03. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-02. Instability detected in spectral estimator, increasing ridge regression factor to 1.00000e-01. Error: Could not stabilize high res. spectrogram Would anyone please help me to solve the problem? I check that there are several emails in the mailist that asked the similar question but unfortunately I can't find the answer there. I appreciate any help and suggestion. Thank you in advance. Regards, Jenni -------------- next part -------------- An HTML attachment was scrubbed... URL: From snoke at vt.edu Thu Nov 17 18:00:48 2011 From: snoke at vt.edu (Arthur Snoke) Date: Thu, 17 Nov 2011 21:00:48 -0500 (EST) Subject: [SAC-HELP] SAC version 101.5 available Message-ID: Version 101.5 of SAC is now available by request from IRIS DMC . For a summary of major changes, please refer to the CHANGES file at . ?New features described in that file include additional plot-file formats, a robust traveltime routine, and the use of annotated polezero files. For a detailed list of enhancements and bug fixes, see the HISTORY file at , There are binary distributions for Oracle Solaris, Linux (32 bit and 64 bit), and Mac OS X (PPC, Intel 32 bit, and Intel 64 bit). ?You can also request the source code and create a build of your own. SAC can be run on Windows machines running Cygwin. ?For licensing reasons IRIS cannot provide a Windows Cygwin binary distribution. Prior to v101.5 IRIS provided a separate source distribution for Cygwin, but recent releases (starting with v1.7) can build on Cygwin with the same source distribution?as the other platforms. The most up to date SAC Users Manual can be accessed online at . Brought to you by: The SAC Development Team From jingyue at mst.edu Sun Nov 20 23:52:39 2011 From: jingyue at mst.edu (Yueling Jing) Date: Mon, 21 Nov 2011 15:52:39 +0800 Subject: [SAC-HELP] question about rdseed and transfer command. Message-ID: Hi, Recently, I transferred velociy record SEED file to SAC file by rdseed, First, I use the command: "rdseed -f 200805120628045.SC.SEED -R -d -o -p 1" so the Poles and Zeros file are outputed, such as: SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 *Can I specify the file name without the last part of number?* such as ?SAC_PZs_SC_ZJG_BHZ_? . Second, I need to remove the instrument response of the data. *Can I use the transfer command as:* "r 2008*ZJG*BHZ*SAC" "Transfer from polezero subtype SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 to VEL" As a result, the output on the terminal is : ********************************** Extracting polezero response for ZJG , BHZ ... Station (ZJG ), Channel (BHZ ) Calper is not available. Cannot tell if response is normalized, so transfer function will be used without scaling. Calib is not available. Cannot tell if scale factor has been applied on input, so transfer function will be used without scaling. Waveform multiplied by 1.000000 after deconvolution. ********************************** The details of SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 is as follow: * ********************************** * NETWORK (KNETWK): SC * STATION (KSTNM): ZJG * LOCATION (KHOLE): 00 * CHANNEL (KCMPNM): BHZ * CREATED : 2011-11-21T06:54:02 * START : 2007-07-01T00:00:00 * END : 2009-07-02T01:01:00 * DESCRIPTION : Jiangyouzhongjiagou,Sichuan Province * LATITUDE : 31.792500 * LONGITUDE : 104.672900 * ELEVATION : 612.0 * DEPTH : 0.0 * DIP : 0.0 * AZIMUTH : 0.0 * SAMPLE RATE : 100.0 * INPUT UNIT : M * OUTPUT UNIT : COUNTS * INSTTYPE : Guralp CMG-3ESPC Seismometer * INSTGAIN : 1.929980e+03 (M/S) * COMMENT : N/A * SENSITIVITY : 1.214590e+09 (M/S) * A0 : 5.715080e+08 * ********************************** ZEROS 3 +0.000000e+00 +0.000000e+00 +0.000000e+00 +0.000000e+00 +0.000000e+00 +0.000000e+00 POLES 5 -7.397840e-02 +7.397840e-02 -7.397840e-02 -7.397840e-02 -1.004800e+03 +0.000000e+00 -5.024000e+02 +0.000000e+00 -1.130400e+03 +0.000000e+00 CONSTANT +6.941478e+17 Thanks. Yueling -------------- next part -------------- An HTML attachment was scrubbed... URL: From sheila at blacknest.gov.uk Mon Nov 21 00:40:54 2011 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Mon, 21 Nov 2011 08:40:54 +0000 Subject: [SAC-HELP] question about rdseed and transfer command. In-Reply-To: References: Message-ID: <4ECA0E96.1090509@blacknest.gov.uk> Dear Yueling Jing, You do not need to worry about the messages from SAC, I think. RDSEED has scaled the "constant" of your poles-and-zeros to incorporate the "sensitivity" into it, so you do not need a separate calibration factor ("calib"), i.e. constant, 6.941478e17 is equal to A0 normalisation factor 5.71508e8 x sensitivity 1.21459e9 m/s/digital count and the A0 normalisation factor was calculated on the assumption that the calibration period (calper) was one second (I have checked this with my program that finds the ratio of the polynomials). You can check that the calibration period really is one second by sending for the data sheet for the seismometer from Guralp. Send an email to caldoc at guralp.com with the serial number of the seismometer in the header. You should receive the data sheet as a Microsoft Word document automatically. The calibration factor "calib", in case you need it for other purposes, is given by 1/(sensitivity in m/s * 2pi), but in metres/count (so multiply by 10^9 if you need it in nanometres/count). I make yours 1.310359e-10 metres/count. If you want to check that the value is correct, you need to get the digitiser data sheet as well as the seismometer data sheet. If it is a Guralp digitiser then you can obtain this from caldoc at guralp.com by putting the serial number of the digitiser on to the subject line of the email (but beware that the resulting Microsoft Word file will not display the values properly in OpenOffice, only in genuine Microsoft Word). Guralp seismometer data sheets give poles and zeros in terms of hertz not radians, and in terms of velocity not displacement, so you will find that there are only two zeros, not three, in the data sheet for the seismometer, and the poles need to be multiplied by 2pi/calper to match the ones reported by rdseed. Regards, Sheila Peacock. On 11/21/2011 07:52 AM, Yueling Jing wrote: > Hi, > > Recently, I transferred velociy record SEED file to SAC file by rdseed, > > First, I use the command: > > "rdseed -f 200805120628045.SC.SEED -R -d -o -p 1" > > so the Poles and Zeros file are outputed, > > such as: SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 > > *Can I specify the file name without the last part of number?* such as ?SAC_PZs_SC_ZJG_BHZ_? . > Second, I need to remove the instrument response of the data. > *Can I use the transfer command as:* > "r 2008*ZJG*BHZ*SAC" > > "Transfer from polezero subtype SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 to VEL" > > As a result, the output on the terminal is : > ********************************** > > Extracting polezero response for ZJG, BHZ... > > Station (ZJG), Channel (BHZ) > > Calper is not available. Cannot tell if response is normalized, so transfer function will be used without scaling. > > Calib is not available. Cannot tell if scale factor has been applied on input, so transfer function will be used without scaling. > > Waveform multiplied by 1.000000 after deconvolution. > > ********************************** > > The details of SAC_PZs_SC_ZJG_BHZ_00_2007.182.00.00.00.0000_2009.182.24.60.60.99999 is as follow: > > * ********************************** > > * NETWORK(KNETWK): SC > > * STATION(KSTNM): ZJG > > * LOCATION(KHOLE): 00 > > * CHANNEL(KCMPNM): BHZ > > * CREATED: 2011-11-21T06:54:02 > > * START: 2007-07-01T00:00:00 > > * END: 2009-07-02T01:01:00 > > * DESCRIPTION: Jiangyouzhongjiagou,Sichuan Province > > * LATITUDE: 31.792500 > > * LONGITUDE: 104.672900 > > * ELEVATION: 612.0 > > * DEPTH: 0.0 > > * DIP: 0.0 > > * AZIMUTH: 0.0 > > * SAMPLE RATE: 100.0 > > * INPUT UNIT: M > > * OUTPUT UNIT: COUNTS > > * INSTTYPE: Guralp CMG-3ESPC Seismometer > > * INSTGAIN: 1.929980e+03 (M/S) > > * COMMENT: N/A > > * SENSITIVITY: 1.214590e+09 (M/S) > > * A0: 5.715080e+08 > > * ********************************** > > ZEROS3 > > +0.000000e+00+0.000000e+00 > > +0.000000e+00+0.000000e+00 > > +0.000000e+00+0.000000e+00 > > POLES5 > > -7.397840e-02+7.397840e-02 > > -7.397840e-02-7.397840e-02 > > -1.004800e+03+0.000000e+00 > > -5.024000e+02+0.000000e+00 > > -1.130400e+03+0.000000e+00 > > CONSTANT+6.941478e+17 > Thanks. > Yueling > > > > _______________________________________________ > 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 Nov 21 03:38:44 2011 From: snoke at vt.edu (Arthur Snoke) Date: Mon, 21 Nov 2011 06:38:44 -0500 (EST) Subject: [SAC-HELP] question about rdseed and transfer command Message-ID: Following up on Sheila's response ... Header variable CONSTANT is not used in the transfer polezero instrument correction, and as of RDSEED v5.1 is not put in the header. There was a bug in SAC prior to the just released v101.5 regarding CONSTANT, so best to get the new version. However, there is a bug in RDSEED v5.1 that may cause problems working with transfer polezero on some operating systems. See for further discussion. Finally, one should included the FREQLIMITIS argument in the TRANSFER call or one gets excessive correction at long periods. The effect is not as bad with a correction to velocity as it is with a correction to dislplacement. The online HELP file discusses this as well. From jingyue at mst.edu Mon Nov 21 04:01:00 2011 From: jingyue at mst.edu (Jing, Yueling) Date: Mon, 21 Nov 2011 06:01:00 -0600 Subject: [SAC-HELP] question about rdseed and transfer command References: Message-ID: <6212FB22FD17584D90840DD8370D7B3711EAD5@mst-vmail01.srv.mst.edu> Hi Arthur & Sheila, Thank you very much. I will download a new version of SAC and include the FREQlimit argument. Best regards. Yueling -----Original Message----- From: sac-help-bounces at iris.washington.edu on behalf of Arthur Snoke Sent: Mon 11/21/2011 5:38 AM To: SAC-help Listserv Subject: Re: [SAC-HELP] question about rdseed and transfer command Following up on Sheila's response ... Header variable CONSTANT is not used in the transfer polezero instrument correction, and as of RDSEED v5.1 is not put in the header. There was a bug in SAC prior to the just released v101.5 regarding CONSTANT, so best to get the new version. However, there is a bug in RDSEED v5.1 that may cause problems working with transfer polezero on some operating systems. See for further discussion. Finally, one should included the FREQLIMITIS argument in the TRANSFER call or one gets excessive correction at long periods. The effect is not as bad with a correction to velocity as it is with a correction to dislplacement. The online HELP file discusses this as well. _______________________________________________ sac-help mailing list sac-help at iris.washington.edu http://www.iris.washington.edu/mailman/listinfo/sac-help From f.darbyshire at gmail.com Wed Nov 23 12:26:11 2011 From: f.darbyshire at gmail.com (Fiona Darbyshire) Date: Wed, 23 Nov 2011 15:26:11 -0500 Subject: [SAC-HELP] Installation of SAC 101.5; operating system issues Message-ID: Hello, Today I obtained and installed SAC 101.5 for a Linux 64-bit system. (I've successfully installed previous versions on this machine). When I tried to run the program, I received the following error message: sac: /lib64/libc.so.6: version `GLIBC_2.7' not found (required by sac) sac: /lib64/libc.so.6: version `GLIBC_2.11' not found (required by sac) Not being any expert on OS libraries, I asked the system administrator for advice. He told me that my machine has operating system CentOS 5.7 which has the version GLIBC 2.5 and one cannot install a more recent version by the distribution update system (apologies if this is clumsy wording; I'm translating directly from French...). CentOS 6 uses GLIBC 2.12 but my machine is 5 years old... So my question is: is there a way to do a successful install of the sac 101.5 binaries without having to upgrade my entire operating system? I'd imagine there's a simple solution out there somewhere, but I'm not sufficiently computer-savvy to figure it out myself! The adminstrator wondered if there's a version with all libraries included that was therefore independent of the OS libraries. Thanks in advance for your advice! Fiona Darbyshire Universit? du Qu?bec ? Montr?al From sheila at blacknest.gov.uk Thu Nov 24 00:56:47 2011 From: sheila at blacknest.gov.uk (Sheila Peacock) Date: Thu, 24 Nov 2011 08:56:47 +0000 Subject: [SAC-HELP] Installation of SAC 101.5; operating system issues In-Reply-To: References: Message-ID: <4ECE06CF.6070002@blacknest.gov.uk> Dear Fiona, I got the same message from the 64-bit version of the pre-compiled SAC. The 32-bit version worked OK on an older 32-bit computer. Sorry I don't know how to get round the GLIBC version problem. I have had to compile up a working version of SAC 101.5 from the source code anyway (because I had to modify one of the subroutines to satisfy local needs here). As these things go it was straightforward. If your system has the Gnu C compiler (type "which gcc" to see if it has) then perhaps you can compile the source code with the default settings, i.e. unpack the tar file of source code, type "configure" then "make" then "make install" (you might need to be "root" to run "make install"). You will have to ask for the source code distribution from IRIS if you didn't get it along with the precompiled binaries. Actually I have simplified this too much. This was the command I issued to invoke the "configure" script: configure --prefix=/sharedprograms/sac/v101.5/64c The "--prefix=/sharedprograms....." directed "make install" to install the completed SAC program into /sharedprograms/sac/v101.5/64c which does not need me to be "root" for access. My CentOS system seemed otherwise to have everything SAC needed (in particular, "libpng.a", in the default places, fortunately. Regards, Sheila Peacock.