SAC Help

A forum for SAC users to exchange ideas and ask questions about SAC usage, installation, portability, etc.

Email address
Moderated by
IRIS Webmaster
Related tags
SAC Software SAC

Active Message Threads for October 2006

andreas wessel
2006-10-11 06:30:39 - 2006-10-27 21:29:11
Good evening, I'm trying to pass some commands to SAC with << But all I get is the SAC startup lines: SEISMIC ANALYSIS CODE [06/06/2005 (Version 100.1)] Copyright 1995 Regents of the University of California .. and a blinking cursor a few lines below. All I can do then is kill the script with Ctrl+C. I am running a minimal example. I get the "echo start", it's only the SAC part that causes trouble. #!/bin/sh echo start sac<<eof r test.e p q eof exit 0 I also tried csh instead of sh, but … [more]
andreas wessel
2006-10-27 20:53:01
Hello, For accessing .sgf files that were created on another machine I probably have to byteswap them. With SAC data files thats done easily by just 2 lines in SAC: r file w file Now I'm wondering, if there is a similar (easy) way to do this with sgf files. The sgftops script only works for sgf files that I created myself, with other ones I end up with a "Segmentation fault" Thanks, Andreas
Anna M. Courtier
2006-10-26 18:29:41
Hi, I'm trying to store a slowness value into one of the SAC header user fields to be used later for binning. I was hoping to use the value calculated in TauP for the event distance/depth specific to each trace and then store it in the header. I saw the man page on 'traveltime' but am not sure if I can extend this to slowness or how to implement it (preferably in the macro I am using to pick arrival times and cut traces, etc). Can this be done? If so, what is the best approach? I am us… [more]
Geo .
2006-10-25 07:00:10
Hi, I've been having a little trouble with SAC cut command lately. Following is an example script that I use to cut data files, cut on r 2006* cut 75 300 cut off w f1 f2 f3 r f1 f2 f3 w over When I read these cut files in & plot them later, the "cut" doesn't come in to effect. In fact I never could get the files to cut for some reason. It keeps on plotting the original file instead of the cut file. xlim is much more effective in this sense but you can't write files on to the disk with that e… [more]
Juan Gomez
2006-10-20 01:14:03
Hi there, The SAC binary work well also on a AMD-64 PC, however, there is another related problem: the library sacio.a, it can not be used on local compiled programs. It is, I have several fortran programs made by my self that work very well on a 32-bits PC, however, I can nnot compile then on my AMD-64 bits PC. An help, comment or suggestion? Thanks, JuanM At 02:01 p.m. 19/10/2006, you wrote: >Send sac-help mailing list submissions … [more]
No replies
2006-10-18 00:46:46 - 2006-10-18 22:24:36
Hi, I'm hoping to install sac on a linux box running 64bit SMP linux >uname -a Linux localhost.localdomain 2.6.15-1.1833_FC4smp #1 SMP Wed Mar 1 23:55:52 EST 2006 x86_64 x86_64 x86_64 GNU/Linux but, in reading the archives of this mail list, it seems SAC won't work in 64-bit mode, and indeed, when I compile the source version I have (05Mar01.sac.tgz), it fails. Although, it seems to fail because it can't find or can't make a proper X11 library: … [more]
Garrett Euler
2006-10-13 06:58:08
Something like this was mentioned a little earlier: I didn't see a resolution or much further discussion so I wanted to point out that while the older version of sac (00.56a) would throw an error on merging if two files overlapped or were less than a sample interval apart (ERROR 1802: Time overlap:), the newer version (100.1) does not spit any error (on both solaris and linux), does nothing with the files in memor… [more]
No replies
Garrett Euler
2006-10-12 20:37:34
I encountered an error (ERROR 2118: No transfer function applied.) when attempting to remove the instrument response with the following polezeros: ZEROS 3 POLES 6 -0.0048 0.0000 -0.0729 0.0000 -18.4725 18.8458 -18.4725 -18.8458 -33.9245 0.0000 -118.8514 0.0000 CONSTANT 4.439785e+15 I dug into it a little and found this only occured when attempted on linux (worked fine on solaris). Both running 100.1. Then I kept digging and found that it will work if you change the CONSTANT value by an… [more]
No replies
Geo .
2006-10-12 05:42:14 - 2006-10-12 18:20:23
Hi, I want to convert horizontal components in to transverse and radial components in SAC. Would anyone have a code/instructions to do this? You feedback is greatly appreciated. Januka. On 10/10/06, <> wrote: > > Send sac-help mailing list submissions to > > > To subscribe or unsubscribe via the World Wide Web, visit > >… [more]
2006-10-09 13:59:26 - 2006-10-10 20:17:59
Good mornning I'd like to know how to get the header of sac files into an ascii file? thank you Regards Héloïse
2006-10-09 23:03:04 - 2006-10-10 16:29:53
Good morning, I'd like to get some information in the header but the subroutines getfhv doesn't seem to work for depmin, depmax and dist in a fortran code: call GETFHV ('DEPMIN', depmin, nerr) I got this error message: ERROR: Undefined header field value. DEPMIN but it works for: call GETNHV ('NZYEAR', nyear, nerr) is there a bug? or an other way to extract data from the header? Best regards Héloïse CETE Nice, LGIT Grenoble
20:15:04 v.22510d55