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 tohttp://www.iris.washington.edu/pipermail/sac-help/attachments/20080307/0f44b481/attachment-0001.html
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:
------------------------------
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
sac> fft
sac> 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
***************************************
-
Onur Tan and Tuna Eken
This is an issue with SAC version 101.0. The blackboard variable
section
of the code base has a bug regarding the concatenation of blackboard
variables.
Please upgrade to 101.1 to resolve this issue.
Cheers
Brian Savage
On Mar 8, 2008, at 5:03 PM , Tuna Eken wrote:
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
http://www.iris.washington.edu/pipermail/sac-help/attachments/
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:
20080307/0f44b481/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
sac> fft
sac> 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