[SAC-HELP] Transverse & Radial Components

Georgia Cua georgia.cua at sed.ethz.ch
Thu Oct 12 02:20:23 PDT 2006


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 at 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
-----------------------------


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20061012/f2f6f959/attachment.html>


More information about the sac-help mailing list