A forum for SAC users to exchange ideas and ask questions about SAC usage, installation, portability, etc.
This list has moved to https://groups.google.com/a/earthscope.org/g/sac-community.
- Email address
- sac-help@lists.ds.iris.edu
- Type
-
Help
Help lists are for submitting support questions. Anyone can post a question for NSF SAGE (or another team) to answer.
- Moderated by
- IRIS Webmaster
- Related tags
- SAC Software SAC
October
September
August
July
June
May
April
March
February
January
September
August
July
June
May
April
March
February
January
December
November
October
September (1)
August
July
June (4)
May (11)
April (3)
March (5)
February (2)
January (8)
November
October
September (1)
August
July
June (4)
May (11)
April (3)
March (5)
February (2)
January (8)
Active Message Threads for October 2006
andreas wessel
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]
Replies
-
Re: Passing commands to SAC using sac<<eof - andreas wessel - 2006-10-27 21:29:11
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
Replies
-
Re: Byteswap .sgf files? - Arthur Snoke - 2006-10-27 15:28:12
Anna M. Courtier
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]
Replies
-
Re: Slowness & Headers - Renate Hartog - 2006-10-26 17:04:35
Geo .
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]
Replies
-
Re: cut command - Michael Thorne - 2006-10-25 03:11:48
Juan Gomez
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
Derek
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]
Replies
-
Re: sac on 64-bit SMP box - Brian Savage - 2006-10-18 17:29:32
-
Re: sac on 64-bit SMP box - Derek Schutt - 2006-10-18 18:21:46
-
Re: sac on 64-bit SMP box - Derek - 2006-10-18 20:15:25
-
Re: sac on 64-bit SMP box - Brian Savage - 2006-10-18 22:24:36
-
-
-
Garrett Euler
Something like this was mentioned a little earlier:
http://www.iris.washington.edu/pipermail/sac-help/2006-August/000111.html
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
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
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, sac-help-request@iris.washington.edu <
sac-help-request@iris.washington.edu> wrote:
>
> Send sac-help mailing list submissions to
> sac-help@iris.washington.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://www.iris.washington.edu/mailman/listinfo/sac-help
>…
[more]
Replies
-
Re: Transverse & Radial Components - Georgia Cua - 2006-10-12 18:20:23
Heloise.Cadet@ujf-grenoble.fr
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
Replies
-
Re: get the header in ascii - andreas wessel - 2006-10-10 20:17:59
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
Replies
-
Re: get header data - Peter Goldstein - 2006-10-10 01:02:14
-
Re: get header data - Frederik Tilmann - 2006-10-10 16:29:53
