Thread: Transverse & Radial Components

I think you want the "rotate" command in sac.

within sac, type "help rotate" to get more info about this command

Georgia

On Oct 12, 2006, at 4:42 AM, Januka Attanayake wrote:

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<at>iris.washington.edu < sac-help-
request<at>iris.washington.edu> wrote:
Send sac-help mailing list submissions to
sac-help<at>iris.washington.edu

To subscribe or unsubscribe via the World Wide Web, visit
http://www.iris.washington.edu/mailman/listinfo/sac-help
or, via email, send a message with subject or body 'help' to
sac-help-request<at>iris.washington.edu

You can reach the person managing the list at
sac-help-owner<at>iris.washington.edu

When replying, please edit your Subject line so it is more specific
than "Re: Contents of sac-help digest..."


Today's Topics:

1. Re: get the header in ascii (andreas wessel)
2. Re: get header data (Peter Goldstein)
3. Re: get header data (Frederik Tilmann)
4. Passing commands to SAC using sac<<eof (andreas wessel)


----------------------------------------------------------------------

Message: 1
Date: Tue, 10 Oct 2006 13:17:59 +1300
From: "andreas wessel" <awbochum<at>gmail.com>
Subject: Re: [SAC-HELP] get the header in ascii
To: sac-help<at>iris.washington.edu
Message-ID:
<2afddd310610091717i60eeb6e8r555d3574b922a3c<at>mail.gmail.com >
Content-Type: text/plain; charset="iso-8859-1"

Hello,

There is a way to write ascii, but it's far from perfect, as it is
kind of
tricky to see which values belong to which variables.
If you have the file in SAC memory, you can write it as ascii with

SAC> WRITE ALPHA filename
and then look at the values in your terminal, the header is at the
beginning
(it shows only the values, but not the names)

$ head filename

Please let me know if you find a better way to do this.


Andreas



On 10/9/06, Heloise.Cadet<at>ujf-grenoble.fr < Heloise.Cadet@ujf-
grenoble.fr>
wrote:

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
_______________________________________________
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: http://www.iris.washington.edu/pipermail/sac-help/attachments/
20061009/9d65a158/attachment-0001.html

------------------------------

Message: 2
Date: Mon, 9 Oct 2006 18:02:14 -0700
From: Peter Goldstein < peterg<at>llnl.gov>
Subject: Re: [SAC-HELP] get header data
To: <Heloise.Cadet<at>ujf-grenoble.fr>
Cc: sac-help<at>iris.washington.edu
Message-ID: <p06110409c1509d5fdfa8@[192.168.1.5]>
Content-Type: text/plain; charset="iso-8859-1" ; format="flowed"


As Art suggests, saclst will hopefully provide
the header list capabilities you need.
In the short term you can hopefully get by using listheader within
a script and
redirect the output to the file you want.

Something like the following script could be invoked as
sac_lh.csh "your_filename"
Note that the blank lines are needed to get all the headers and you'll
still have to edit the results if you want a particular format or
to get rid of
extraneous information and blank lines.


#!/bin/csh -xvf
# sac_lh.csh: script to list sac header listing to standard out
sac<<eof
r $argv[1]
lh inclusive on




q
eof
# end of sac_lh.csh

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



_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help

--

Peter Goldstein, Ph.D. (925) 423-1231 (office)
L-103, PO Box 808 (925) 422-5844 (fax)
Livermore, CA 94551 peterg<at>llnl.gov (email)
web page: http://www.llnl.gov/sac



------------------------------

Message: 3
Date: Tue, 10 Oct 2006 09:29:53 +0100
From: Frederik Tilmann <tilmann<at>esc.cam.ac.uk>
Subject: Re: [SAC-HELP] get header data
To: Heloise.Cadet<at>ujf-grenoble.fr
Cc: sac-help<at>iris.washington.edu
Message-ID: < 452B5A01.80609<at>esc.cam.ac.uk>
Content-Type: text/plain; charset=us-ascii; format=flowed

Heloise

these values all need to be calculated by sac so when the sac file is
generated externally, e.g. by rdseed they are generally left
undefined.

Try reading into sac and writing again, i.e.

sac<<eof
r file
w over
q
eof

Of course, for dist evla,evlo,stla,stlo all need to be set and lcalda
must be true.

Frederik


Heloise.Cadet<at>ujf-grenoble.fr wrote:

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



_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help

--
Frederik Tilmann
Bullard Laboratories Tel. +44 1223 765545
Department of Earth Sciences Fax. +44 1223 360779
University of Cambridge email: tilmann<at>esc.cam.ac.uk
Madingley Road
Cambridge CB3 0EZ
UK






------------------------------

Message: 4
Date: Tue, 10 Oct 2006 23:30:39 +1300
From: "andreas wessel" <awbochum<at>gmail.com>
Subject: [SAC-HELP] Passing commands to SAC using sac<<eof
To: sac-help<at>iris.washington.edu
Message-ID:
<2afddd310610100330t118cd90cg717e7c22a91d94b<at>mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"

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 that didn't help.
Are there any environmental things besides SACAUX I need to consider?

Thanks,
Andreas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.iris.washington.edu/pipermail/sac-help/attachments/
20061010/3ef772f6/attachment-0001.html

------------------------------

_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help


End of sac-help Digest, Vol 16, Issue 2
***************************************



_______________________________________________
sac-help mailing list
sac-help<at>iris.washington.edu
http://www.iris.washington.edu/mailman/listinfo/sac-help

-----------------------------
Dr. Georgia Cua
Swiss Seismological Service
Institute of Geophysics
ETH Hoenggerberg HPP P
CH-8093 Zurich
Switzerland

Tel: +41-044-633-7574
Fax: +41-044-633-1065

georgia.cua<at>sed.ethz.ch
-----------------------------