Thread: RE: Re: sac-help Digest, Vol 33, Issue 4

Hi Tuna,

exactly, my data length is 1800 sec.
if you look the sample again, cutting limits is 600-700 sec.
the problem is the first variable ( %a )
when echo on, the second variable can not writen to the screen !
SAC> cut %a %b
cut %a %b
==> cut 600


Do you use the version 101.0?
I have no problem with old version.

onur


-----Original Message-----
From: sac-help-bounces<at>iris.washington.edu
[sac-help-bounces<at>iris.washington.edu] On Behalf Of Tuna Eken
Sent: Sunday, March 09, 2008 12:04 AM
To: sac-help<at>iris.washington.edu
Subject: [SAC-HELP] Re: sac-help Digest, Vol 33, Issue 4





To message 3:

Hi Onur,

when I tried the same commands with the same parameters onto my example
sac file it did worked out ... It actually gave the same error messae
when i used different time values (whis is not included by the data) ...

So problem may be caused by applying wrong values for your data. have
you ever checked starting end ending points of the data ?

regards

tuna eken


Message: 3
Date: Sat, 8 Mar 2008 13:29:51 +0200
From: "Onur Tan" <Onur.Tan<at>mam.gov.tr>
Subject: [SAC-HELP] Blackboard Variables and CUT
To: <sac-help<at>iris.washington.edu>
Message-ID:
<110C3E0C4E9ADB4C882229C2A7A1540A052C2511<at>POSTA2K3.mam.gov.tr>
Content-Type: text/plain; charset="us-ascii"


Hi,

When I use a blackboard variable as the first parameter of the cut
command in SAC new version (101.0), it gives error.

I tried this (echo on):
--------------------
SAC> setbb a 600
setbb a 600
SAC> setbb b 700
setbb b 700
SAC> cut %a %b
cut %a %b
==> cut 600
SAC> r
r
WARNING: Start cut greater than stop cut for file 08010_0707_CEV.1
ERROR 1301: No data files read in.
SAC>
-----------------------------
the cut command do not understand the second variable.


This does also not work:
----
SAC> cut %a 700
cut %a 700
==> cut 600
SAC> r
r
WARNING: Start cut greater than stop cut for file 08010_0707_CEV.1
ERROR 1301: No data files read in.
----


But this works
----
SAC> cut 500 %b
cut 500 %b
==> cut 500 700
SAC> r
r
SAC> p
p
----


Do you have any idea?

best regards
Onur Tan




Citerar sac-help-request<at>iris.washington.edu:

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. Two SAC queries - (i) conversion of SAC files to

PITSA-ASCII

and (ii) SAC library issues for multi-channel

cross-correlation

(Fiona Darbyshire)
2. About fft routine in SAC (wwxu)
3. Blackboard Variables and CUT (Onur Tan)


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

Message: 1
Date: Fri, 7 Mar 2008 18:13:09 -0500
From: "Fiona Darbyshire" <f.darbyshire<at>gmail.com>
Subject: [SAC-HELP] Two SAC queries - (i) conversion of SAC files to
PITSA-ASCII and (ii) SAC library issues for multi-channel
cross-correlation
To: sac-help<at>iris.washington.edu
Message-ID:
<256ff3740803071513m5ef1bc36nbfdb2c6317b6d87c<at>mail.gmail.com>
Content-Type: text/plain; charset="iso-8859-1"

Hello,

I have 2 SAC-related queries that are causing me problems on a couple
of fairly urgent projects currently:

i) I inherited a Matlab program for surface wave analysis. The
seismogram files that it reads are in PITSA ASCII format. Is there any

existing software to convert standard SAC to PITSA ASCII? Or will I
have to convert SAC binary to alpha, and try to reformat the resulting

ascii files?

ii) Another inheritance: this time I have been given the source code
(fortran) for John VanDecar's "mccc" multi-channel cross-correlation
program. The program contains a lot of SAC functions, such as "rsach",

"getnhv", "rsac1", "getkhv", "getfhv", "setfhv", "setkhv", "wsac0"...
(and probably a few others I haven't spotted yet). Along with the
Fortran codes I was given a couple of SAC libraries of the type
"libsac.a". The problem is that this set worked OK for a SUN system,
but I am on Linux (CentOS 5). With the SAC library that I was given,
and the SAC library I obtained when downloading SAC executables for my

current system, I was unable to compile the code. I have a really old
SAC library (old enough that I can't remember even which version it

comes from) and that one allowed mccc to compile...

but not to run (segmentation fault). Can anyone advise me on using SAC

libraries correctly and/or making changes to the fortran code that
will allow me to run this program?

I should add that I'm a rather poor programmer - just about literate
in Fortran, okay in simple shell scripts, but not too great with any
other programming languages! Any advice greatly appreciated.

Thanks a lot,

Fiona Darbyshire.
-------------- next part -------------- An HTML attachment was
scrubbed...
URL:

http://www.iris.washington.edu/pipermail/sac-help/attachments/20080307/0
f44b481/attachment-0001.html

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

Message: 2
Date: Sat, 08 Mar 2008 02:59:07 GMT
From: "wwxu" <wwxu<at>mail.iggcas.ac.cn>
Subject: [SAC-HELP] About fft routine in SAC
To: sac-help<at>iris.washington.edu
Message-ID: <20080308025907.3888.eqmail<at>mail.iggcas.ac.cn>
Content-Type: text/plain; format=flowed; charset="gb2312"

Hi ALL,

In SAC, I generate a 5Hz sine wave and do fft transforming as follows:
sac> funcgen sine 5 npts 1024 delta 0.01 fft psp am loglin
and I find that the maximum absolute amplitude is greater than 4 near
frequency 5Hz, but in theory, it should be 0.5 because of amplitude
spectrum symmetry of real series, so could you tell me why?

B.T.W., I also compute the amplitude spectrum using fft program I
wrote and find that the result is not the same with SAC generated, I
checked the program carefully and be sure that it's right and conform
with the theoretical result.

thanks for your attention!


Your sincerely,
Xu W.

************************************************
Dr. WeiWei Xu
Institute of Geology and Geophysics
Chinese Academy of Sciences
P.O.Box 9825
Beijing 100029
China
Tel: 0086-10-62007354(Office)
Fax: 0086-10-62010846
************************************************



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

Message: 3
Date: Sat, 8 Mar 2008 13:29:51 +0200
From: "Onur Tan" <Onur.Tan<at>mam.gov.tr>
Subject: [SAC-HELP] Blackboard Variables and CUT
To: <sac-help<at>iris.washington.edu>
Message-ID:
<110C3E0C4E9ADB4C882229C2A7A1540A052C2511<at>POSTA2K3.mam.gov.tr>
Content-Type: text/plain; charset="us-ascii"


Hi,

When I use a blackboard variable as the first parameter of the cut
command in SAC new version (101.0), it gives error.

I tried this (echo on):
--------------------
SAC> setbb a 600
setbb a 600
SAC> setbb b 700
setbb b 700
SAC> cut %a %b
cut %a %b
==> cut 600
SAC> r
r
WARNING: Start cut greater than stop cut for file 08010_0707_CEV.1
ERROR 1301: No data files read in.
SAC>
-----------------------------
the cut command do not understand the second variable.


This does also not work:
----
SAC> cut %a 700
cut %a 700
==> cut 600
SAC> r
r
WARNING: Start cut greater than stop cut for file 08010_0707_CEV.1
ERROR 1301: No data files read in.
----


But this works
----
SAC> cut 500 %b
cut 500 %b
==> cut 500 700
SAC> r
r
SAC> p
p
----


Do you have any idea?

best regards
Onur Tan




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

_______________________________________________
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 33, Issue 4
***************************************

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