From =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=  Mon Jun  1 00:16:27 2009
From: =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= (=?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=)
Date: Mon, 1 Jun 2009 10:16:27 +0300
Subject: [SAC-HELP] Initialization macro
Message-ID: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>

Hi,

I am trying to run an initialization macro as sac starts up.
I have created an init.m macro file in the /sac/macro directory and I  
am able to run it from within sac by executing the following two  
commands:

SETMACRO sac/macro
MACRO init.m

I tried aliasing sac like this:

alias sac='sac sac/macro/init.m'

but it doesn't work and I don't think that it should but this is what  
it says to do in the README file for SAC so I tried.

How do I make something run as I call up sac?

Thanks

Shahar

From meixue.geo at gmail.com  Mon Jun  1 01:50:19 2009
From: meixue.geo at gmail.com (Mei Xue)
Date: Mon, 1 Jun 2009 16:50:19 +0800
Subject: [SAC-HELP] Initialization macro
In-Reply-To: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>
References: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>
Message-ID: <9cda3feb0906010150p3460929w23f2019f022364a9@mail.gmail.com>

Hi Shahar,

You may try:
alias sac 'sac m sac/macro/init.m'

see if it works.

Mei


2009/6/1 Shahar Kadmiel ??? ?????? <kadmiel at bgu.ac.il>:
> Hi,
>
> I am trying to run an initialization macro as sac starts up.
> I have created an init.m macro file in the /sac/macro directory and I am
> able to run it from within sac by executing the following two commands:
>
> SETMACRO sac/macro
> MACRO init.m
>
> I tried aliasing sac like this:
>
> alias sac='sac sac/macro/init.m'
>
> but it doesn't work and I don't think that it should but this is what it
> says to do in the README file for SAC so I tried.
>
> How do I make something run as I call up sac?
>
> Thanks
>
> Shahar
> _______________________________________________
> sac-help mailing list
> sac-help at iris.washington.edu
> http://www.iris.washington.edu/mailman/listinfo/sac-help
>


From vijay.srtmu at gmail.com  Mon Jun  1 02:30:33 2009
From: vijay.srtmu at gmail.com (Vijay Tulasingam)
Date: Mon, 1 Jun 2009 15:00:33 +0530
Subject: [SAC-HELP] Sac in Mac
Message-ID: <2e484aaf0906010230w69b9b24ocf5cc3b388a103b8@mail.gmail.com>

Hi All

I am new user to MAC, i have installed sac as in Linux. i have set the
SACAUX variable and the path. I am unable to see the graphics window with
plot command.

Thank you
Vijay Kumar
-- 
T. Vijay Kumar
Lecturer
School of Earth Sciences
SRTMU
Nanded
Maharashtra
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090601/231ec5f8/attachment.htm>

From snokeja at vt.edu  Mon Jun  1 04:58:33 2009
From: snokeja at vt.edu (Arthur Snoke)
Date: Mon, 01 Jun 2009 07:58:33 -0400
Subject: [SAC-HELP] Initialization macro
In-Reply-To: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>
References: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>
Message-ID: <4A23C269.3010609@vt.edu>

In the README file, the absolute path is set for the init.m file.  Try using the 
absolute path and see if it finds your init.m.

Do you have a directory ${SACHOME}/macro?  If so, try

alias sac '${SACHOME}/sac ${SACHOME}/macro/init.m'

I think your 'SETMACRO sac/macro' command requires that SAC is already running 
before it will be recognized.

Shahar Kadmiel ??? ?????? wrote:
> Hi,
> 
> I am trying to run an initialization macro as sac starts up.
> I have created an init.m macro file in the /sac/macro directory and I am 
> able to run it from within sac by executing the following two commands:
> 
> SETMACRO sac/macro
> MACRO init.m
> 
> I tried aliasing sac like this:
> 
> alias sac='sac sac/macro/init.m'
> 
> but it doesn't work and I don't think that it should but this is what it 
> says to do in the README file for SAC so I tried.
> 
> How do I make something run as I call up sac?
> 
> Thanks
> 
> Shahar

From nolet at geoazur.unice.fr  Mon Jun  1 05:41:59 2009
From: nolet at geoazur.unice.fr (Guust Nolet)
Date: Mon, 1 Jun 2009 14:41:59 +0200
Subject: [SAC-HELP] Initialization macro
In-Reply-To: <4A23C269.3010609@vt.edu>
References: <FB18493B-C9AA-4F81-9754-B7831886AB70@bgu.ac.il>
	<4A23C269.3010609@vt.edu>
Message-ID: <D1130F0B-DD1F-4AC6-9A0F-48C93D284396@geoazur.unice.fr>

In my case the alias works, but I take care to specify the *complete*  
address of the macro. I use tcsh:

alias sac '/usr/local/sac/bin/sac /Users/guust/.sacrc'

Guust Nolet
nolet at geoazur.unice.fr

Geosciences Azur
250, Rue Albert Einstein
06560 Sophia Antipolis
France
+33.4.92.94.26.32

http://geoweb.princeton.edu/people/faculty/nolet/



On 1 juin 09, at 13:58, Arthur Snoke wrote:

> In the README file, the absolute path is set for the init.m file.   
> Try using the absolute path and see if it finds your init.m.
>
> Do you have a directory ${SACHOME}/macro?  If so, try
>
> alias sac '${SACHOME}/sac ${SACHOME}/macro/init.m'
>
> I think your 'SETMACRO sac/macro' command requires that SAC is  
> already running before it will be recognized.
>
> Shahar Kadmiel ??? ?????? wrote:
>> Hi,
>> I am trying to run an initialization macro as sac starts up.
>> I have created an init.m macro file in the /sac/macro directory and  
>> I am able to run it from within sac by executing the following two  
>> commands:
>> SETMACRO sac/macro
>> MACRO init.m
>> I tried aliasing sac like this:
>> alias sac='sac sac/macro/init.m'
>> but it doesn't work and I don't think that it should but this is  
>> what it says to do in the README file for SAC so I tried.
>> How do I make something run as I call up sac?
>> Thanks
>> Shahar
> _______________________________________________
> sac-help mailing list
> sac-help at iris.washington.edu
> http://www.iris.washington.edu/mailman/listinfo/sac-help


From =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=  Tue Jun  9 06:14:06 2009
From: =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= (=?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=)
Date: Tue, 9 Jun 2009 16:14:06 +0300
Subject: [SAC-HELP] READSP, WRITESP and spectral ratio
Message-ID: <BD31ACF6-49EA-485B-8EFF-1D4645127475@bgu.ac.il>

  I am trying to compute spectral ratios of between a hefty number of  
files and I try it like this:

read /path_to_files/filenames*.sac
fft
writesp over

read /path_to_other_files/otherfiles*.sac
fft
writesp over

To this point every thing works well and I can even view with plotsp  
and see what I am doing.

I now want to divide filenames by other_files hence get the spectral  
ratio between them.

the * stands for sequential numbers and the number of file is the same  
in the nominator and denominator.

how do you suggest I go at this?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090609/ca853714/attachment.htm>

From liszews2 at tcnj.edu  Tue Jun  9 12:03:00 2009
From: liszews2 at tcnj.edu (Kelly M Liszewski)
Date: Tue, 9 Jun 2009 15:03:00 -0400 (EDT)
Subject: [SAC-HELP] ppk zoom error
Message-ID: <745776526.1562861244574180404.JavaMail.root@zcs.TCNJ.EDU>

 I am using SAC version 101.2 on intel Mac OS X version 10.5.6, and I am running ppk to pick p waves for several stations at a time.  I type ppk p 4 to view four stations at a time, but I can only zoom properly in the first subwindow (at the top of the first four traces plotted).  As I scroll through all of the stations to pick them, when I try to zoom, the zoom function zooms in on the wrong portion of the trace.  Has anyone else had this problem, and does anyone have any advice? 
        Thank you,
        Kelly Liszewski
        The College of New Jersey

From =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=  Wed Jun 10 06:56:45 2009
From: =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= (=?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=)
Date: Wed, 10 Jun 2009 16:56:45 +0300
Subject: [SAC-HELP] Spectral data
Message-ID: <7351640C-B9C2-4AD2-857D-95A0EAA554F9@bgu.ac.il>

How can I divide a spectral data file by another spectral data file?

this doesn't work:

read sac1.sac sac2.sac
fft
writesp over

readsp sac1.sac
binoperr npts ignore delta ignore
divf sac2.sac

I get the following error:

ERROR 1307: Illegal operation on spectral file sac1.sac


From Margaret_Doll at brown.edu  Thu Jun 11 12:51:00 2009
From: Margaret_Doll at brown.edu (Margaret Doll)
Date: Thu, 11 Jun 2009 15:51:00 -0400
Subject: [SAC-HELP] floating point exception
Message-ID: <3743D2A7-8234-4441-AB3F-7BA17CE9B466@brown.edu>

We have just installed sac 101.2 on a linux 64 bit computer running   
RedHat.   The SAC distribution was not a source distribution.

We are getting a "Floating point error" when we run the program.  I  
assume that our 64 bit libraries are causing the problem in that SAC  
is a 32 bit program.

How can we fix the problem?

From savage at uri.edu  Thu Jun 11 13:21:01 2009
From: savage at uri.edu (Brian Savage)
Date: Thu, 11 Jun 2009 16:21:01 -0400
Subject: [SAC-HELP] floating point exception
In-Reply-To: <3743D2A7-8234-4441-AB3F-7BA17CE9B466@brown.edu>
References: <3743D2A7-8234-4441-AB3F-7BA17CE9B466@brown.edu>
Message-ID: <D868A44F-AB0A-41F9-AAFE-F62BCE62B493@uri.edu>

Margaret,

Does the linux machine have the 32 bit libraries installed ?
You can check which libraries sac is using by issuing the command

% ldd /usr/local/sac/bin/sac

Assuming sac is installed in the default location.
My system has 32 and 64 bit libraries installed at
32 bit library
	/lib/ld-linux.so.2
64 bit library
	/lib64/ld-linux-x86-64.so.2

The sac binary should be using for the first one.

By the way, the 64 bit version of sac will be coming out soon.

Brian Savage
savage at uri.edu


On Jun 11, 2009, at  3:51 PM , Margaret Doll wrote:

> We have just installed sac 101.2 on a linux 64 bit computer  
> running  RedHat.   The SAC distribution was not a source distribution.
>
> We are getting a "Floating point error" when we run the program.  I  
> assume that our 64 bit libraries are causing the problem in that  
> SAC is a 32 bit program.
>
> How can we fix the problem?
> _______________________________________________
> sac-help mailing list
> sac-help at iris.washington.edu
> http://www.iris.washington.edu/mailman/listinfo/sac-help


From =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=  Sun Jun 14 08:16:08 2009
From: =?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?= (=?UTF-8?Q?Shahar_Kadmiel_=D7=A9=D7=97=D7=A8_=D7=A7=D7=93=D7=9E?=)
Date: Sun, 14 Jun 2009 18:16:08 +0300
Subject: [SAC-HELP] how should I cite SAC in a paper?
Message-ID: <E697BC67-B697-4E0D-9561-828EB1336436@bgu.ac.il>



From peterg at llnl.gov  Mon Jun 15 11:30:11 2009
From: peterg at llnl.gov (Peter Goldstein)
Date: Mon, 15 Jun 2009 11:30:11 -0700
Subject: [SAC-HELP] how should I cite SAC in a paper?
In-Reply-To: <E697BC67-B697-4E0D-9561-828EB1336436@bgu.ac.il>
References: <E697BC67-B697-4E0D-9561-828EB1336436@bgu.ac.il>
Message-ID: <7be0l7$1eb03v@nspiron-1.llnl.gov>




Dear Shahar Kadamiel,

I would use the following references.

Cheers,

Peter

Goldstein, P., A. Snoke, (2005), ?SAC 
Availability for the IRIS Community?, 
Incorporated Institutions for Seismology Data 
Management Center Electronic 
Newsletter.  http://www.iris.edu/news/newsletter/vol7no1/page1.htm


Goldstein, P., D. Dodge, M. Firpo, Lee Minner 
(2003) ?SAC2000: Signal processing  and analysis 
tools for seismologists and engineers, Invited 
contribution to ?The IASPEI International 
Handbook of Earthquake and Engineering 
Seismology?, Edited by WHK Lee, H. Kanamori, P.C. 
Jennings, and C. Kisslinger, Academic Press, London.


At 08:16 AM 6/14/2009, Shahar Kadmiel ?????? ?????? wrote:

>_______________________________________________
>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)
  Lawrence Livermore National Security
  L-103, PO Box 808        (925) 422-5844 (fax)
  Livermore, CA 94551      peterg at llnl.gov (email)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090615/7255363d/attachment.htm>

From courtiam at jmu.edu  Tue Jun 16 08:31:10 2009
From: courtiam at jmu.edu (Anna Courtier)
Date: Tue, 16 Jun 2009 11:31:10 -0400
Subject: [SAC-HELP] Header Calculations
Message-ID: <F19B8694-A121-464A-9F2B-6553D755987C@jmu.edu>

Hello,

I have some SAC files and need az, baz, dist, and gcarc in the  
headers. The event and station locations are already in the headers  
and LCALDA = TRUE.  With LCALDA=TRUE, I think the az, baz, etc should  
be automatically calculated and put into the headers, but this is only  
the case for about half of the data files. Switching LCALDA to false  
does not change anything, and resetting the LCALDA  to true does not  
fix the problem either.  Does anyone have a solution for this?  There  
is no obvious pattern for which files work and which do not.

I'm working with SEISMIC ANALYSIS CODE [02/05/2008 (Version 101.1)] on  
Mac OS X.

Thanks,
ANNA







From Onur.Tan at mam.gov.tr  Tue Jun 16 11:12:44 2009
From: Onur.Tan at mam.gov.tr (Onur Tan)
Date: Tue, 16 Jun 2009 21:12:44 +0300
Subject: [SAC-HELP] ascii to binary
Message-ID: <9E88105FE4AC8C47BD6D00EB9831482F2ABFE40C@posta.mam.gov.tr>

Hi,

I need a simple program that convert  ascii file to binary SAC file format.
Can you help me?

best regards
Onur

 
Dr. Onur TAN
----------------------------------------------   40.7866N 29.4500E ---------
 T?B?TAK Marmara Ara?t?rma Merkezi,  Yer ve Deniz Bilimleri Enstit?s?
 TUBITAK Marmara Research Center, Earth and Marine Sciences Institute
 Gebze - Kocaeli - TURKEY 

From a1981117 at hotmail.com  Thu Jun 25 15:51:28 2009
From: a1981117 at hotmail.com (guorongwen)
Date: Fri, 26 Jun 2009 06:51:28 +0800
Subject: [SAC-HELP] FW: Does not work
Message-ID: <SNT119-W154442F09E97333E21A8BB98340@phx.gbl>


Hi, All
        I have installed the sac in the ubuntu. But it seems does not work. I just input the 'sac' in the command line and it is said 'bash:sac: command not found''. But I use echo $SACAUX command, it showed /usr/local/sac/aux. what is the problem? I am a newie, and here nobody know about it. Thanks very much!
        Best wishes!
 
Your truly
Wen



??????????MSN??? ?????

_________________________________________________________________
????????????MClub?????????
http://club.msn.cn/?from=10
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090626/dcdfc7ad/attachment.htm>

From snoke at vt.edu  Thu Jun 25 16:02:46 2009
From: snoke at vt.edu (Arthur Snoke)
Date: Thu, 25 Jun 2009 19:02:46 -0400 (EDT)
Subject: [SAC-HELP] FW: Does not work
In-Reply-To: <SNT119-W154442F09E97333E21A8BB98340@phx.gbl>
References: <SNT119-W154442F09E97333E21A8BB98340@phx.gbl>
Message-ID: <alpine.OSX.2.00.0906251902290.29631@macarthur-2.local>

Is the executable sac in your path?

On Fri, 26 Jun 2009, guorongwen wrote:

> Hi, All
> ????????I have installed the sac in the ubuntu. But it seems does not work.
> I just input the 'sac' in the command line and it is said 'bash:sac: command
> not found''.?But I?use?echo $SACAUX?command,?it showed /usr/local/sac/aux.
> what?is the problem? I?am a newie, and here nobody know about it. Thanks
> very much!
> ??????? Best wishes!
> ?
> Your truly
> Wen
> 
> ____________________________________________________________________________
> ??????????MSN??? ?????
> 
> ____________________________________________________________________________
> ?MSN???????????????? ?????
>

From a1981117 at hotmail.com  Thu Jun 25 21:14:29 2009
From: a1981117 at hotmail.com (guorongwen)
Date: Fri, 26 Jun 2009 12:14:29 +0800
Subject: [SAC-HELP] Thanks
Message-ID: <SNT119-W53BC1CF0889D0BF998A36F98350@phx.gbl>


Hi, All,

        Thanks for your help. Now it works.There is another sac in the system which caused the problem(As  Arthur recommended). I really appreciate it for the help of Arthur and Charlotte. 

           Best wishes!

 

Your truly

Wen

_________________________________________________________________
?Windows Live ??????????Messenger?
http://www.windowslive.cn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090626/797410df/attachment.htm>

From wangwtustc at gmail.com  Mon Jun 29 06:37:40 2009
From: wangwtustc at gmail.com (weitao wang)
Date: Mon, 29 Jun 2009 21:37:40 +0800
Subject: [SAC-HELP] how to rotate R T comp into N E ?
Message-ID: <406b435e0906290637w34dcc633ie98c09bf4b668980@mail.gmail.com>

Hi,All
       I have a question while rotating horizontal component with sac. We
can use "rotate to GCP" to rotate N E into  R T comp,but how to rotate back
into N E components? Is there an easy way to do this?
      thanks for your help.
      Best wishes

weitao
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090629/cb9e4f49/attachment.htm>

From snoke at vt.edu  Mon Jun 29 07:23:05 2009
From: snoke at vt.edu (Arthur Snoke)
Date: Mon, 29 Jun 2009 10:23:05 -0400 (EDT)
Subject: [SAC-HELP] how to rotate R T comp into N E ?
In-Reply-To: <406b435e0906290637w34dcc633ie98c09bf4b668980@mail.gmail.com>
References: <406b435e0906290637w34dcc633ie98c09bf4b668980@mail.gmail.com>
Message-ID: <alpine.OSX.2.00.0906291018200.3393@macarthur-2.local>

In the help command for rotate, there are options for rotating "to" and 
"through" a fixed angle.  It would be easy to write a macro that could 
give you back the N and E components.  I would do test cases to make sure 
it works correctly for you.

Here is a macro I use to do the "normal" rotation:

r $1$.n $1$.e
rtr
rotate to gcp
write $1$.r $1$.t
r $1$.r
chnhdr kcmpnm "radial"
write over
r $1$.t
chnhdr kcmpnm "trnsvers"
write over

If it is named rotaate.m, one calls it for  abc.n and abc.e using

SAC> m rotate.m abc

On Mon, 29 Jun 2009, weitao wang wrote:

> Hi,All
> I have a question while rotating horizontal component with sac. We can 
> use "rotate to GCP" to rotate N E into? R T comp,but how to rotate back 
> into N E components? Is there an easy way to do this?

From bjlemon at 163.com  Mon Jun 29 21:25:50 2009
From: bjlemon at 163.com (bjlemon)
Date: Tue, 30 Jun 2009 12:25:50 +0800 (CST)
Subject: [SAC-HELP] Fw:Fwd:  how to rotate R T comp into N E ?
Message-ID: <19790075.682291246335950793.JavaMail.coremail@bj163app84.163.com>

Hi all,
Thank you for your letter. 
In your letter, you have written how to rotate the N E comp to R T. But my problem is how to rotate the R T comp to N E. I did see the help of sac command "rotate", but there is no instruction abut it? So I wrote a letter again.


Best wishes,

Ning Liu and Weitao Wang
 
 
 





---------- Forwarded message ----------
From: Arthur Snoke<snoke at vt.edu>
Date: 2009/6/29
Subject: Re: [SAC-HELP] how to rotate R T comp into N E ?
To: weitao wang <wangwtustc at gmail.com>
Cc: sac-help at iris.washington.edu


In the help command for rotate, there are options for rotating "to" and "through" a fixed angle.  It would be easy to write a macro that could give you back the N and E components.  I would do test cases to make sure it works correctly for you.

Here is a macro I use to do the "normal" rotation:

r $1$.n $1$.e
rtr
rotate to gcp
write $1$.r $1$.t
r $1$.r
chnhdr kcmpnm "radial"
write over
r $1$.t
chnhdr kcmpnm "trnsvers"
write over

If it is named rotaate.m, one calls it for  abc.n and abc.e using

SAC> m rotate.m abc


On Mon, 29 Jun 2009, weitao wang wrote:

Hi,All
I have a question while rotating horizontal component with sac. We can use "rotate to GCP" to rotate N E into  R T comp,but how to rotate back into N E components? Is there an easy way to do this?






200????,????,?????


200????,????,?????
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20090630/6a40a954/attachment.htm>

From renate at ess.washington.edu  Mon Jun 29 22:09:15 2009
From: renate at ess.washington.edu (Renate Hartog)
Date: Mon, 29 Jun 2009 22:09:15 -0700 (PDT)
Subject: [SAC-HELP] Fw:Fwd:  how to rotate R T comp into N E ?
In-Reply-To: <19790075.682291246335950793.JavaMail.coremail@bj163app84.163.com>
References: <19790075.682291246335950793.JavaMail.coremail@bj163app84.163.com>
Message-ID: <Pine.LNX.4.64.0906292154040.27095@argon.ess.washington.edu>

Hello,

Below is a simple example to show how to use rotate to go from any 
orientation to another (orthogonal orientation). "rotate to number" 
changes the azimuth of the first listed component to "number" degrees 
clockwise from N. The second component will be number+90 from N.

You can check the orientation of your two files by listing the value of 
the cmpaz header value (clockwise from N).

-Renate

SAC> r
06100802482.D04A.BHE 06100802482.D04A.BHN
SAC> lh cmpaz


   FILE: 06100802482.D04A.BHE - 1
  --------------------------

      cmpaz = 9.000000e+01

   FILE: 06100802482.D04A.BHN - 2
  --------------------------

      cmpaz = 0.000000e+00
SAC> rotate to gcp
SAC> lh


   FILE: 06100802482.D04A.BHE - 1
  --------------------------

      cmpaz = 2.871565e+02

   FILE: 06100802482.D04A.BHN - 2
  --------------------------

      cmpaz = 1.715649e+01
SAC> rotate to 0.
SAC> lh


   FILE: 06100802482.D04A.BHE - 1
  --------------------------

      cmpaz = 0.000000e+00

   FILE: 06100802482.D04A.BHN - 2
  --------------------------

      cmpaz = 9.000000e+01
SAC>



On Tue, 30 Jun 2009, bjlemon wrote:

> Hi all,
> Thank you for your letter.
> In your letter, you have written how to rotate the N E comp to R T. But my problem is how to rotate the R T comp to N E. I did see the help of sac command "rotate", but there is no instruction abut it? So I wrote a letter again.
>
>
> Best wishes,
>
> Ning Liu and Weitao Wang
>
>
>
>
>
>
>
>
> ---------- Forwarded message ----------
> From: Arthur Snoke<snoke at vt.edu>
> Date: 2009/6/29
> Subject: Re: [SAC-HELP] how to rotate R T comp into N E ?
> To: weitao wang <wangwtustc at gmail.com>
> Cc: sac-help at iris.washington.edu
>
>
> In the help command for rotate, there are options for rotating "to" and "through" a fixed angle.  It would be easy to write a macro that could give you back the N and E components.  I would do test cases to make sure it works correctly for you.
>
> Here is a macro I use to do the "normal" rotation:
>
> r $1$.n $1$.e
> rtr
> rotate to gcp
> write $1$.r $1$.t
> r $1$.r
> chnhdr kcmpnm "radial"
> write over
> r $1$.t
> chnhdr kcmpnm "trnsvers"
> write over
>
> If it is named rotaate.m, one calls it for  abc.n and abc.e using
>
> SAC> m rotate.m abc
>
>
> On Mon, 29 Jun 2009, weitao wang wrote:
>
> Hi,All
> I have a question while rotating horizontal component with sac. We can use "rotate to GCP" to rotate N E into  R T comp,but how to rotate back into N E components? Is there an easy way to do this?
>
>
>
>
>
>
> 200????????,????????,??????????
>
>
> 200????????,????????,??????????