Thread: Instrument responce files

Started: 2013-11-27 14:03:22
Last activity: 2014-02-18 20:03:10
Topics: SAC Help
Shekar
2013-11-27 14:03:22
Hai all

Any body has the instrument response file(PZ file) for the following combination

1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

Thank you very much.

shekar


On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



Message: 1

Date: Tue, 26 Nov 2013 06:45:34 -0300

From: Milton Plasencia

Subject: Re: [SAC-HELP] Insert hedder

To: shekar

Cc: sac-help<at>iris.washington.edu

Message-ID:

Content-Type: text/plain; charset="us-ascii"



Hi,

Use

ch O GMT 2009 318 13 03 34 300



where 318 is the Julian day corresponding to month 11 and day 14.



HTH,



Milton P. Plasencia Linares

mpplasencia<at>gmail.com







On Nov 26, 2013, at 6:34 AM, shekar wrote:



Dear all,



I have a SAC file converted from SEISAN.



How to insert the origin time(2009-11-14-13:03:34.300) of the event in sac header?



Thanks in advance



Shekar





On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



2. Remove Response Instrument (amik Amik)



3. Re: Remove Response Instrument (Arthur Snoke)



4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







Message: 1



Date: Wed, 13 Nov 2013 09:03:30 +0000



From: George Helffrich



Subject: Re: [SAC-HELP] sort and plot commands don't play



To: Frederik Tilmann



Cc: sac-help<at>iris.washington.edu



Message-ID:



Content-Type: text/plain; charset=us-ascii







Dear All -







Plotting the traces in a record section in SSS displays traces in increasing distance order.







On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







I have a problem in that the plot commands (p1,ppk) do not show the files in the same order as they are listed in memory after a sort.







What I am trying to do is to display a north- to south record section of noise correlation functions with respect to a target station. The coordinates of both stations are place in evla and stla, but it is not well defined in which variable a given station ends up in (evla, or stla).







So to plot all the functions relative to J10S I use







r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



sort evla ascend stla ascend



lh evla stla







Everything looks good, stations are in the right order in the lh list.







But if I then do



p1



or ppk



they come out in the wrong order, plotting the J03-J10 after J10-J11 (correct would have been J09-J10 as in the lh listing). I attach a screenshot demonstrating this effect.



I tried



sort commit evla ascend stla ascend



but this this did not make a difference.







This is on 101.5c



SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











Any ideas would be appreciated







Regards



Frederik



















Output of lh command







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac - 1



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







evla = 5.253600e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac - 2



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







evla = 5.270155e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac - 3



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







evla = 5.287174e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac - 4



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







evla = 5.305374e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac - 5



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







evla = 5.320426e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac - 6



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







evla = 5.339244e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac - 7



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







evla = 5.355617e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac - 8



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







evla = 5.366659e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac - 9



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







evla = 5.390106e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac - 10



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







evla = 5.411258e+01



Waiting







stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac - 11



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







evla = 5.427260e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac - 12



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







evla = 5.444247e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac - 13



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







evla = 5.460217e+01



stla = 5.481430e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac - 14



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







evla = 5.481430e+01



stla = 5.640510e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac - 15



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







evla = 5.481430e+01



stla = 5.617126e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac - 16



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







evla = 5.481430e+01



stla = 5.603592e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac - 17



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







evla = 5.481430e+01



stla = 5.577234e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac - 18



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







evla = 5.481430e+01



stla = 5.530829e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac - 19



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







evla = 5.481430e+01



stla = 5.512811e+01







FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac - 20



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







evla = 5.481430e+01



Waiting







stla = 5.497167e+01



SAC>



_______________________________________________



sac-help mailing list



sac-help<at>iris.washington.edu



http://www.iris.washington.edu/mailman/listinfo/sac-help















George Helffrich



george.helffrich<at>bris.ac.uk



















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







Message: 2



Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



From: amik Amik



Subject: [SAC-HELP] Remove Response Instrument



To: "sac-help<at>iris.washington.edu"



Message-ID:







Content-Type: text/plain; charset="iso-8859-1"







I am working on a local magnitude research. I want to remove response instrument



from the data and simulate the response of Wood Anderson seismograph. What should I do?



I have a few questions :



1. Which one the true procedure and command in SAC :



?? a) SAC> trans from evalresp to none



?? ? ?and then



?? ? ?SAC> trans from polezero subtype "file name.pz" to none







?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







?? c) SAC> trans to wa







?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







?? f) SAC> trans from evalresp fname RESP.SATION to wa







?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq 0.005 0.01 1.0 2.0







2. What the different utilizing of Response File and Polezero File?



?? I think polezero already included in responfile, is it right ?







May you help me? Thank You very much











Sincerely Yours



Ami



-------------- next part --------------



An HTML attachment was scrubbed...



URL:







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







Message: 3



Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



From: Arthur Snoke



Subject: Re: [SAC-HELP] Remove Response Instrument



To: amik Amik



Cc: "sac-help<at>iris.washington.edu"



Message-ID:



Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







I think the best way for you to understand what the various options do is



to try them andsave the results. The TRANSFER help file at



has answers



to your explicit quesitons. If anything is not clear, let me know.







On Sun, 17 Nov 2013, amik Amik wrote:







I am working on a local magnitude research. I want to remove response



instrument



from the data and simulate the response of Wood Anderson seismograph. What



should I do?



I have a few questions :



1. Which one the true procedure and command in SAC :



?? a) SAC> trans from evalresp to none



?? ? ?and then



?? ? ?SAC> trans from polezero subtype "file name.pz" to none







?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







?? c) SAC> trans to wa







?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







?? f) SAC> trans from evalresp fname RESP.SATION to wa







?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq



0.005 0.01 1.0 2.0







2. What the different utilizing of Response File and Polezero File?



?? I think polezero already included in responfile, is it right ?







May you help me? Thank You very much











Sincerely Yours



Ami















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







Message: 4



Date: Sun, 17 Nov 2013 14:29:15 +0100



From: "Milton P. Plasencia Linares"



Subject: Re: [SAC-HELP] Remove Response Instrument



To: "amik Amik"



Cc: "sac-help<at>iris.washington.edu"



Message-ID:



Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



format="flowed"







look below,







**********************************



Milton P. PLASENCIA LINARES







Dipartimento Centro di Ricerche Sismologiche



Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







Borgo Grotta Gigante 42/C



(34010) Sgonico - TRIESTE - ITALIA



Tel: +39-040-2140 256/141



Fax: +39-040-327307







E-mail: mplasencia<at>inogs.it







ASAIN (Antarctic Seismographic Argentinean Italian Network)



*********************************











Quoting "amik Amik" :







I am working on a local magnitude research. I want to remove



response instrument



from the data and simulate the response of Wood Anderson



seismograph. What should I do?



I have a few questions :



1. Which one the true procedure and command in SAC :



?? a) SAC> trans from evalresp to none



?? ? ?and then



?? ? ?SAC> trans from polezero subtype "file name.pz" to none







?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







?? c) SAC> trans to wa







?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



Option d) This is fine if you use the polezero station file.







?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







?? f) SAC> trans from evalresp fname RESP.SATION to wa







?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



freq 0.005 0.01 1.0 2.0







2. What the different utilizing of Response File and Polezero File?



Response file include the filter (FIR) coefficients







?? I think polezero already included in responfile, is it right ?



Yes.







May you help me? Thank You very much







My suggest ius use:







transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















Sincerely Yours



Ami







Cheers,







Milton















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



This message was sent using IMP, the Internet Messaging Program.



















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







Message: 5



Date: Mon, 18 Nov 2013 08:01:33 -0500



From: Fiona Darbyshire



Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



To: sac-help<at>iris.washington.edu



Message-ID:







Content-Type: text/plain; charset=ISO-8859-1







Hello,







Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







Up to now, e.g. when binning receiver functions, I have used BAZ and



GCARC to group them together, and that's always seemed to work, so



I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



you a unique position from the point of interest (the station).







This time, I have some data output that is expressed in AZ and DIST. I



assumed (perhaps naively?) that this would also allow successful



spatial grouping, but it turns out not to be the case; one can have



the same forward azimuth for several different locations.







Here's the example where I found the issue:



Station lat, lon: 82.5033, -62.35 (ALE)







Event1: lat 80.205, lon -1.091



Event 2: lat 86.876, lon 54.373



The two events are almost 1000km apart.







Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







I'd like to know more about the calculation, in particular because (i)



I'd like to reassure myself that BAZ/distance does indeed give a



unique location, and (ii) I don't understand how AZ is calculated. One



website I looked at (Matlab) talks about rhumb lines vs great-circles;



is that the issue? Will an AZ/distance pair always be non-unique?



(Also most websites define back-azimuth as simply 180 opposite to



azimuth, but I guess they aren't talking about spherical geometry



then...)







Thanks for any insight.







Fiona Darbyshire.







Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







Message: 6



Date: Mon, 18 Nov 2013 11:38:44 -0500



From: Brian Savage



Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



To: SAC HELP



Message-ID:



Content-Type: text/plain; charset="windows-1252"







Fiona,







The implementation of the computation of Great Circle Distance (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's formula found in Geodesy by R.Bomford (1980). It assumes a Radius of 6378.160 km and a flattening of 1/298.2466081. I can send you the appropriate pages from R. Bomford (1980) if you would like.







The azimuth and back-azimuth will generally not be offset by180 degrees, and some pairs of point will have "odd-looking" (not event close to 180 degrees) azimuth/back azimuth pairs.







Note: There are more robust methods of calculating these quantities and the current implementation in SAC as well as many other implementations, up until recently [1], have trouble with particular pairs of points. For over ~ 99% of point pairs SAC's implementation produces the correct result, but that last 1% can be incorrect (sometimes badly).







I hope this helps.







Brian Savage







[1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z (open access). Addenda.











On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







Hello,







Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







Up to now, e.g. when binning receiver functions, I have used BAZ and



GCARC to group them together, and that's always seemed to work, so



I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



you a unique position from the point of interest (the station).







This time, I have some data output that is expressed in AZ and DIST. I



assumed (perhaps naively?) that this would also allow successful



spatial grouping, but it turns out not to be the case; one can have



the same forward azimuth for several different locations.







Here's the example where I found the issue:



Station lat, lon: 82.5033, -62.35 (ALE)







Event1: lat 80.205, lon -1.091



Event 2: lat 86.876, lon 54.373



The two events are almost 1000km apart.







Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







I'd like to know more about the calculation, in particular because (i)



I'd like to reassure myself that BAZ/distance does indeed give a



unique location, and (ii) I don't understand how AZ is calculated. One



website I looked at (Matlab) talks about rhumb lines vs great-circles;



is that the issue? Will an AZ/distance pair always be non-unique?



(Also most websites define back-azimuth as simply 180 opposite to



azimuth, but I guess they aren't talking about spherical geometry



then...)







Thanks for any insight.







Fiona Darbyshire.







Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







_______________________________________________



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:







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







_______________________________________________



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 97, Issue 2



***************************************







Get your own FREE website, FREE domain & FREE mobile app with Company email.

Know More >

_______________________________________________

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:



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



_______________________________________________

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 97, Issue 7

***************************************


  • Milton Plasencia
    2013-11-27 13:47:32
    Hi,

    You can consider using the PDCC software by IRIS.
    http://www.iris.edu/dms/nodes/dmc/software/downloads/pdcc/

    For Guralp sensor you must have the TTL value (Tap TableLookup)
    for apply the correct FIR filter combination). look here
    http://www.guralp.com/determining-fir-filter-information-from-gcf-data/

    HTH,

    m.

    * * * * * * * * * * * * * * * * * * * * * * *
    Milton P. PLASENCIA LINARES

    Centro di Ricerche Sismologiche (CRS)
    OGS - Istituto Nazionale di Oceanografia
    e di Geofisica Sperimentale

    Borgo Grotta Gigante 42/C
    (34010) Sgonico - Trieste - Italia
    Tel: +39 040 2140 141 (Udine)
    Tel: +39 040 2140 256 (Trieste)
    Cel.: +39 331 6481 935
    Fax: +39-040-327307

    E-mail: mplasencia<at>inogs.it

    ASAIN (Antarctic Seismographic Argentinean Italian Network)
    *********************************










    On Nov 27, 2013, at 3:03 AM, shekar <m_shekar2k2<at>rediffmail.com> wrote:

    Hai all

    Any body has the instrument response file(PZ file) for the following combination

    1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
    2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

    Thank you very much.

    shekar


    On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



    Message: 1

    Date: Tue, 26 Nov 2013 06:45:34 -0300

    From: Milton Plasencia

    Subject: Re: [SAC-HELP] Insert hedder

    To: shekar

    Cc: sac-help<at>iris.washington.edu

    Message-ID:

    Content-Type: text/plain; charset="us-ascii"



    Hi,

    Use

    ch O GMT 2009 318 13 03 34 300



    where 318 is the Julian day corresponding to month 11 and day 14.



    HTH,



    Milton P. Plasencia Linares

    mpplasencia<at>gmail.com







    On Nov 26, 2013, at 6:34 AM, shekar wrote:



    Dear all,



    I have a SAC file converted from SEISAN.



    How to insert the origin time(2009-11-14-13:03:34.300) of the event in sac header?



    Thanks in advance



    Shekar





    On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



    2. Remove Response Instrument (amik Amik)



    3. Re: Remove Response Instrument (Arthur Snoke)



    4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



    5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



    6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







    Message: 1



    Date: Wed, 13 Nov 2013 09:03:30 +0000



    From: George Helffrich



    Subject: Re: [SAC-HELP] sort and plot commands don't play



    To: Frederik Tilmann



    Cc: sac-help<at>iris.washington.edu



    Message-ID:



    Content-Type: text/plain; charset=us-ascii







    Dear All -







    Plotting the traces in a record section in SSS displays traces in increasing distance order.







    On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







    I have a problem in that the plot commands (p1,ppk) do not show the files in the same order as they are listed in memory after a sort.







    What I am trying to do is to display a north- to south record section of noise correlation functions with respect to a target station. The coordinates of both stations are place in evla and stla, but it is not well defined in which variable a given station ends up in (evla, or stla).







    So to plot all the functions relative to J10S I use







    r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



    sort evla ascend stla ascend



    lh evla stla







    Everything looks good, stations are in the right order in the lh list.







    But if I then do



    p1



    or ppk



    they come out in the wrong order, plotting the J03-J10 after J10-J11 (correct would have been J09-J10 as in the lh listing). I attach a screenshot demonstrating this effect.



    I tried



    sort commit evla ascend stla ascend



    but this this did not make a difference.







    This is on 101.5c



    SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











    Any ideas would be appreciated







    Regards



    Frederik



















    Output of lh command







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac - 1



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







    evla = 5.253600e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac - 2



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







    evla = 5.270155e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac - 3



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







    evla = 5.287174e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac - 4



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







    evla = 5.305374e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac - 5



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







    evla = 5.320426e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac - 6



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







    evla = 5.339244e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac - 7



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







    evla = 5.355617e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac - 8



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







    evla = 5.366659e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac - 9



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







    evla = 5.390106e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac - 10



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







    evla = 5.411258e+01



    Waiting







    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac - 11



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







    evla = 5.427260e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac - 12



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







    evla = 5.444247e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac - 13



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







    evla = 5.460217e+01



    stla = 5.481430e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac - 14



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







    evla = 5.481430e+01



    stla = 5.640510e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac - 15



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







    evla = 5.481430e+01



    stla = 5.617126e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac - 16



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







    evla = 5.481430e+01



    stla = 5.603592e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac - 17



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







    evla = 5.481430e+01



    stla = 5.577234e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac - 18



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







    evla = 5.481430e+01



    stla = 5.530829e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac - 19



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







    evla = 5.481430e+01



    stla = 5.512811e+01







    FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac - 20



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







    evla = 5.481430e+01



    Waiting







    stla = 5.497167e+01



    SAC>



    _______________________________________________



    sac-help mailing list



    sac-help<at>iris.washington.edu



    http://www.iris.washington.edu/mailman/listinfo/sac-help















    George Helffrich



    george.helffrich<at>bris.ac.uk



















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







    Message: 2



    Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



    From: amik Amik



    Subject: [SAC-HELP] Remove Response Instrument



    To: "sac-help<at>iris.washington.edu"



    Message-ID:



    <1384675653.16429.YahooMailNeo<at>web190801.mail.sg3.yahoo.com>



    Content-Type: text/plain; charset="iso-8859-1"







    I am working on a local magnitude research. I want to remove response instrument



    from the data and simulate the response of Wood Anderson seismograph. What should I do?



    I have a few questions :



    1. Which one the true procedure and command in SAC :



    ?? a) SAC> trans from evalresp to none



    ?? ? ?and then



    ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







    ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







    ?? c) SAC> trans to wa







    ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







    ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







    ?? f) SAC> trans from evalresp fname RESP.SATION to wa







    ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq 0.005 0.01 1.0 2.0







    2. What the different utilizing of Response File and Polezero File?



    ?? I think polezero already included in responfile, is it right ?







    May you help me? Thank You very much











    Sincerely Yours



    Ami



    -------------- next part --------------



    An HTML attachment was scrubbed...



    URL:







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







    Message: 3



    Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



    From: Arthur Snoke



    Subject: Re: [SAC-HELP] Remove Response Instrument



    To: amik Amik



    Cc: "sac-help<at>iris.washington.edu"



    Message-ID:



    Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







    I think the best way for you to understand what the various options do is



    to try them andsave the results. The TRANSFER help file at



    has answers



    to your explicit quesitons. If anything is not clear, let me know.







    On Sun, 17 Nov 2013, amik Amik wrote:







    I am working on a local magnitude research. I want to remove response



    instrument



    from the data and simulate the response of Wood Anderson seismograph. What



    should I do?



    I have a few questions :



    1. Which one the true procedure and command in SAC :



    ?? a) SAC> trans from evalresp to none



    ?? ? ?and then



    ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







    ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







    ?? c) SAC> trans to wa







    ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







    ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







    ?? f) SAC> trans from evalresp fname RESP.SATION to wa







    ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq



    0.005 0.01 1.0 2.0







    2. What the different utilizing of Response File and Polezero File?



    ?? I think polezero already included in responfile, is it right ?







    May you help me? Thank You very much











    Sincerely Yours



    Ami















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







    Message: 4



    Date: Sun, 17 Nov 2013 14:29:15 +0100



    From: "Milton P. Plasencia Linares"



    Subject: Re: [SAC-HELP] Remove Response Instrument



    To: "amik Amik"



    Cc: "sac-help<at>iris.washington.edu"



    Message-ID: <20131117142915.910179rx5iy7g6nf<at>webmail.ogs.trieste.it>



    Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



    format="flowed"







    look below,







    **********************************



    Milton P. PLASENCIA LINARES







    Dipartimento Centro di Ricerche Sismologiche



    Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







    Borgo Grotta Gigante 42/C



    (34010) Sgonico - TRIESTE - ITALIA



    Tel: +39-040-2140 256/141



    Fax: +39-040-327307







    E-mail: mplasencia<at>inogs.it







    ASAIN (Antarctic Seismographic Argentinean Italian Network)



    *********************************











    Quoting "amik Amik" :







    I am working on a local magnitude research. I want to remove



    response instrument



    from the data and simulate the response of Wood Anderson



    seismograph. What should I do?



    I have a few questions :



    1. Which one the true procedure and command in SAC :



    ?? a) SAC> trans from evalresp to none



    ?? ? ?and then



    ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







    ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







    ?? c) SAC> trans to wa







    ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



    Option d) This is fine if you use the polezero station file.







    ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







    ?? f) SAC> trans from evalresp fname RESP.SATION to wa







    ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



    freq 0.005 0.01 1.0 2.0







    2. What the different utilizing of Response File and Polezero File?



    Response file include the filter (FIR) coefficients







    ?? I think polezero already included in responfile, is it right ?



    Yes.







    May you help me? Thank You very much







    My suggest ius use:







    transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















    Sincerely Yours



    Ami







    Cheers,







    Milton















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



    This message was sent using IMP, the Internet Messaging Program.



















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







    Message: 5



    Date: Mon, 18 Nov 2013 08:01:33 -0500



    From: Fiona Darbyshire



    Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



    To: sac-help<at>iris.washington.edu



    Message-ID:







    Content-Type: text/plain; charset=ISO-8859-1







    Hello,







    Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







    Up to now, e.g. when binning receiver functions, I have used BAZ and



    GCARC to group them together, and that's always seemed to work, so



    I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



    you a unique position from the point of interest (the station).







    This time, I have some data output that is expressed in AZ and DIST. I



    assumed (perhaps naively?) that this would also allow successful



    spatial grouping, but it turns out not to be the case; one can have



    the same forward azimuth for several different locations.







    Here's the example where I found the issue:



    Station lat, lon: 82.5033, -62.35 (ALE)







    Event1: lat 80.205, lon -1.091



    Event 2: lat 86.876, lon 54.373



    The two events are almost 1000km apart.







    Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



    Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







    I'd like to know more about the calculation, in particular because (i)



    I'd like to reassure myself that BAZ/distance does indeed give a



    unique location, and (ii) I don't understand how AZ is calculated. One



    website I looked at (Matlab) talks about rhumb lines vs great-circles;



    is that the issue? Will an AZ/distance pair always be non-unique?



    (Also most websites define back-azimuth as simply 180 opposite to



    azimuth, but I guess they aren't talking about spherical geometry



    then...)







    Thanks for any insight.







    Fiona Darbyshire.







    Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







    Message: 6



    Date: Mon, 18 Nov 2013 11:38:44 -0500



    From: Brian Savage



    Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



    To: SAC HELP



    Message-ID: <8D3AE955-FCD5-48D4-8643-4E4E25115DEF<at>uri.edu>



    Content-Type: text/plain; charset="windows-1252"







    Fiona,







    The implementation of the computation of Great Circle Distance (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's formula found in Geodesy by R.Bomford (1980). It assumes a Radius of 6378.160 km and a flattening of 1/298.2466081. I can send you the appropriate pages from R. Bomford (1980) if you would like.







    The azimuth and back-azimuth will generally not be offset by180 degrees, and some pairs of point will have "odd-looking" (not event close to 180 degrees) azimuth/back azimuth pairs.







    Note: There are more robust methods of calculating these quantities and the current implementation in SAC as well as many other implementations, up until recently [1], have trouble with particular pairs of points. For over ~ 99% of point pairs SAC's implementation produces the correct result, but that last 1% can be incorrect (sometimes badly).







    I hope this helps.







    Brian Savage







    [1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z (open access). Addenda.











    On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







    Hello,







    Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







    Up to now, e.g. when binning receiver functions, I have used BAZ and



    GCARC to group them together, and that's always seemed to work, so



    I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



    you a unique position from the point of interest (the station).







    This time, I have some data output that is expressed in AZ and DIST. I



    assumed (perhaps naively?) that this would also allow successful



    spatial grouping, but it turns out not to be the case; one can have



    the same forward azimuth for several different locations.







    Here's the example where I found the issue:



    Station lat, lon: 82.5033, -62.35 (ALE)







    Event1: lat 80.205, lon -1.091



    Event 2: lat 86.876, lon 54.373



    The two events are almost 1000km apart.







    Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



    Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







    I'd like to know more about the calculation, in particular because (i)



    I'd like to reassure myself that BAZ/distance does indeed give a



    unique location, and (ii) I don't understand how AZ is calculated. One



    website I looked at (Matlab) talks about rhumb lines vs great-circles;



    is that the issue? Will an AZ/distance pair always be non-unique?



    (Also most websites define back-azimuth as simply 180 opposite to



    azimuth, but I guess they aren't talking about spherical geometry



    then...)







    Thanks for any insight.







    Fiona Darbyshire.







    Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







    _______________________________________________



    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:







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







    _______________________________________________



    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 97, Issue 2



    ***************************************







    Get your own FREE website, FREE domain & FREE mobile app with Company email.

    Know More >

    <2009-11-14-1304-50S.KGD___003_KGD___HH_Z_SAC>_______________________________________________

    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:



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



    _______________________________________________

    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 97, Issue 7

    ***************************************



    Get your own FREE website, FREE domain & FREE mobile app with Company email.
    Know More >
    _______________________________________________
    sac-help mailing list
    sac-help<at>iris.washington.edu
    http://www.iris.washington.edu/mailman/listinfo/sac-help


    • raja babu
      2014-02-18 17:09:37
      Hi all,

      normalization constant in radian = (normalization constant in Hz)2π^
      (numberof poles−number of zeroes)

      for displacement normalization constant in above formula, number of zeroes
      will include one additional zero than velocity normalization
      constant.?????????

      I have confusion. can anyone clear it

      any help will be appreciated..

      and thanks for help...


      On Wed, Nov 27, 2013 at 2:17 PM, Milton Plasencia <mpplasencia<at>gmail.com>wrote:

      Hi,

      You can consider using the PDCC software by IRIS.
      http://www.iris.edu/dms/nodes/dmc/software/downloads/pdcc/

      For Guralp sensor you must have the TTL value (Tap TableLookup)
      for apply the correct FIR filter combination). look here
      http://www.guralp.com/determining-fir-filter-information-from-gcf-data/

      HTH,

      m.

      * * * * * * * * * * * * * * * * * * * * * * *
      Milton P. PLASENCIA LINARES

      Centro di Ricerche Sismologiche (CRS)
      OGS - Istituto Nazionale di Oceanografia
      e di Geofisica Sperimentale

      Borgo Grotta Gigante 42/C
      (34010) Sgonico - Trieste - Italia
      Tel: +39 040 2140 141 (Udine)
      Tel: +39 040 2140 256 (Trieste)
      Cel.: +39 331 6481 935
      Fax: +39-040-327307

      E-mail: mplasencia<at>inogs.it

      ASAIN (Antarctic Seismographic Argentinean Italian Network)
      *********************************










      On Nov 27, 2013, at 3:03 AM, shekar <m_shekar2k2<at>rediffmail.com> wrote:

      Hai all

      Any body has the instrument response file(PZ file) for the following
      combination

      1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
      2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

      Thank you very much.

      shekar


      On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



      Message: 1

      Date: Tue, 26 Nov 2013 06:45:34 -0300

      From: Milton Plasencia

      Subject: Re: [SAC-HELP] Insert hedder

      To: shekar

      Cc: sac-help<at>iris.washington.edu

      Message-ID:

      Content-Type: text/plain; charset="us-ascii"



      Hi,

      Use

      ch O GMT 2009 318 13 03 34 300



      where 318 is the Julian day corresponding to month 11 and day 14.



      HTH,



      Milton P. Plasencia Linares

      mpplasencia<at>gmail.com







      On Nov 26, 2013, at 6:34 AM, shekar wrote:



      Dear all,



      I have a SAC file converted from SEISAN.



      How to insert the origin time(2009-11-14-13:03:34.300) of the event in
      sac header?



      Thanks in advance



      Shekar





      On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



      2. Remove Response Instrument (amik Amik)



      3. Re: Remove Response Instrument (Arthur Snoke)



      4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



      5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



      6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







      Message: 1



      Date: Wed, 13 Nov 2013 09:03:30 +0000



      From: George Helffrich



      Subject: Re: [SAC-HELP] sort and plot commands don't play



      To: Frederik Tilmann



      Cc: sac-help<at>iris.washington.edu



      Message-ID:



      Content-Type: text/plain; charset=us-ascii







      Dear All -







      Plotting the traces in a record section in SSS displays traces in
      increasing distance order.







      On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







      I have a problem in that the plot commands (p1,ppk) do not show the
      files in the same order as they are listed in memory after a sort.







      What I am trying to do is to display a north- to south record section
      of noise correlation functions with respect to a target station. The
      coordinates of both stations are place in evla and stla, but it is not well
      defined in which variable a given station ends up in (evla, or stla).







      So to plot all the functions relative to J10S I use







      r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac
      sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



      sort evla ascend stla ascend



      lh evla stla







      Everything looks good, stations are in the right order in the lh list.







      But if I then do



      p1



      or ppk



      they come out in the wrong order, plotting the J03-J10 after J10-J11
      (correct would have been J09-J10 as in the lh listing). I attach a
      screenshot demonstrating this effect.



      I tried



      sort commit evla ascend stla ascend



      but this this did not make a difference.







      This is on 101.5c



      SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











      Any ideas would be appreciated







      Regards



      Frederik



















      Output of lh command







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac - 1



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







      evla = 5.253600e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac - 2



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







      evla = 5.270155e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac - 3



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







      evla = 5.287174e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac - 4



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







      evla = 5.305374e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac - 5



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







      evla = 5.320426e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac - 6



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







      evla = 5.339244e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac - 7



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







      evla = 5.355617e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac - 8



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







      evla = 5.366659e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac - 9



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







      evla = 5.390106e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac -
      10



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







      evla = 5.411258e+01



      Waiting







      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac -
      11



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







      evla = 5.427260e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac -
      12



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







      evla = 5.444247e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac -
      13



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







      evla = 5.460217e+01



      stla = 5.481430e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac -
      14



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







      evla = 5.481430e+01



      stla = 5.640510e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac -
      15



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







      evla = 5.481430e+01



      stla = 5.617126e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac -
      16



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







      evla = 5.481430e+01



      stla = 5.603592e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac -
      17



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







      evla = 5.481430e+01



      stla = 5.577234e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac -
      18



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







      evla = 5.481430e+01



      stla = 5.530829e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac -
      19



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







      evla = 5.481430e+01



      stla = 5.512811e+01







      FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac -
      20



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







      evla = 5.481430e+01



      Waiting







      stla = 5.497167e+01



      SAC>



      _______________________________________________



      sac-help mailing list



      sac-help<at>iris.washington.edu



      http://www.iris.washington.edu/mailman/listinfo/sac-help















      George Helffrich



      george.helffrich<at>bris.ac.uk



















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







      Message: 2



      Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



      From: amik Amik



      Subject: [SAC-HELP] Remove Response Instrument



      To: "sac-help<at>iris.washington.edu"



      Message-ID:



      <1384675653.16429.YahooMailNeo<at>web190801.mail.sg3.yahoo.com>



      Content-Type: text/plain; charset="iso-8859-1"







      I am working on a local magnitude research. I want to remove response
      instrument



      from the data and simulate the response of Wood Anderson seismograph.
      What should I do?



      I have a few questions :



      1. Which one the true procedure and command in SAC :



      ?? a) SAC> trans from evalresp to none



      ?? ? ?and then



      ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







      ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







      ?? c) SAC> trans to wa







      ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







      ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to
      none







      ?? f) SAC> trans from evalresp fname RESP.SATION to wa







      ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq
      0.005 0.01 1.0 2.0







      2. What the different utilizing of Response File and Polezero File?



      ?? I think polezero already included in responfile, is it right ?







      May you help me? Thank You very much











      Sincerely Yours



      Ami



      -------------- next part --------------



      An HTML attachment was scrubbed...



      URL:







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







      Message: 3



      Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



      From: Arthur Snoke



      Subject: Re: [SAC-HELP] Remove Response Instrument



      To: amik Amik



      Cc: "sac-help<at>iris.washington.edu"



      Message-ID:



      Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







      I think the best way for you to understand what the various options do
      is



      to try them andsave the results. The TRANSFER help file at



      has answers



      to your explicit quesitons. If anything is not clear, let me know.







      On Sun, 17 Nov 2013, amik Amik wrote:







      I am working on a local magnitude research. I want to remove response



      instrument



      from the data and simulate the response of Wood Anderson seismograph.
      What



      should I do?



      I have a few questions :



      1. Which one the true procedure and command in SAC :



      ?? a) SAC> trans from evalresp to none



      ?? ? ?and then



      ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







      ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







      ?? c) SAC> trans to wa







      ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







      ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa
      to none







      ?? f) SAC> trans from evalresp fname RESP.SATION to wa







      ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"
      freq



      0.005 0.01 1.0 2.0







      2. What the different utilizing of Response File and Polezero File?



      ?? I think polezero already included in responfile, is it right ?







      May you help me? Thank You very much











      Sincerely Yours



      Ami















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







      Message: 4



      Date: Sun, 17 Nov 2013 14:29:15 +0100



      From: "Milton P. Plasencia Linares"



      Subject: Re: [SAC-HELP] Remove Response Instrument



      To: "amik Amik"



      Cc: "sac-help<at>iris.washington.edu"



      Message-ID: <20131117142915.910179rx5iy7g6nf<at>webmail.ogs.trieste.it>



      Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



      format="flowed"







      look below,







      **********************************



      Milton P. PLASENCIA LINARES







      Dipartimento Centro di Ricerche Sismologiche



      Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







      Borgo Grotta Gigante 42/C



      (34010) Sgonico - TRIESTE - ITALIA



      Tel: +39-040-2140 256/141



      Fax: +39-040-327307







      E-mail: mplasencia<at>inogs.it







      ASAIN (Antarctic Seismographic Argentinean Italian Network)



      *********************************











      Quoting "amik Amik" :







      I am working on a local magnitude research. I want to remove



      response instrument



      from the data and simulate the response of Wood Anderson



      seismograph. What should I do?



      I have a few questions :



      1. Which one the true procedure and command in SAC :



      ?? a) SAC> trans from evalresp to none



      ?? ? ?and then



      ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







      ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







      ?? c) SAC> trans to wa







      ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



      Option d) This is fine if you use the polezero station file.







      ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa
      to none







      ?? f) SAC> trans from evalresp fname RESP.SATION to wa







      ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



      freq 0.005 0.01 1.0 2.0







      2. What the different utilizing of Response File and Polezero File?



      Response file include the filter (FIR) coefficients







      ?? I think polezero already included in responfile, is it right ?



      Yes.







      May you help me? Thank You very much







      My suggest ius use:







      transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















      Sincerely Yours



      Ami







      Cheers,







      Milton















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



      This message was sent using IMP, the Internet Messaging Program.



















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







      Message: 5



      Date: Mon, 18 Nov 2013 08:01:33 -0500



      From: Fiona Darbyshire



      Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



      To: sac-help<at>iris.washington.edu



      Message-ID:







      Content-Type: text/plain; charset=ISO-8859-1







      Hello,







      Could someone tell me a bit more about how AZ and BAZ are calculated in
      SAC?







      Up to now, e.g. when binning receiver functions, I have used BAZ and



      GCARC to group them together, and that's always seemed to work, so



      I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



      you a unique position from the point of interest (the station).







      This time, I have some data output that is expressed in AZ and DIST. I



      assumed (perhaps naively?) that this would also allow successful



      spatial grouping, but it turns out not to be the case; one can have



      the same forward azimuth for several different locations.







      Here's the example where I found the issue:



      Station lat, lon: 82.5033, -62.35 (ALE)







      Event1: lat 80.205, lon -1.091



      Event 2: lat 86.876, lon 54.373



      The two events are almost 1000km apart.







      Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



      Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







      I'd like to know more about the calculation, in particular because (i)



      I'd like to reassure myself that BAZ/distance does indeed give a



      unique location, and (ii) I don't understand how AZ is calculated. One



      website I looked at (Matlab) talks about rhumb lines vs great-circles;



      is that the issue? Will an AZ/distance pair always be non-unique?



      (Also most websites define back-azimuth as simply 180 opposite to



      azimuth, but I guess they aren't talking about spherical geometry



      then...)







      Thanks for any insight.







      Fiona Darbyshire.







      Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







      Message: 6



      Date: Mon, 18 Nov 2013 11:38:44 -0500



      From: Brian Savage



      Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



      To: SAC HELP



      Message-ID: <8D3AE955-FCD5-48D4-8643-4E4E25115DEF<at>uri.edu>



      Content-Type: text/plain; charset="windows-1252"







      Fiona,







      The implementation of the computation of Great Circle Distance
      (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's
      formula found in Geodesy by R.Bomford (1980). It assumes a Radius of
      6378.160 km and a flattening of 1/298.2466081. I can send you the
      appropriate pages from R. Bomford (1980) if you would like.







      The azimuth and back-azimuth will generally not be offset by180 degrees,
      and some pairs of point will have "odd-looking" (not event close to 180
      degrees) azimuth/back azimuth pairs.







      Note: There are more robust methods of calculating these quantities and
      the current implementation in SAC as well as many other implementations, up
      until recently [1], have trouble with particular pairs of points. For over
      ~ 99% of point pairs SAC's implementation produces the correct result, but
      that last 1% can be incorrect (sometimes badly).







      I hope this helps.







      Brian Savage







      [1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of
      Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K.
      doi:10.1007/s00190-012-0578-z (open access). Addenda.











      On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







      Hello,







      Could someone tell me a bit more about how AZ and BAZ are calculated
      in SAC?







      Up to now, e.g. when binning receiver functions, I have used BAZ and



      GCARC to group them together, and that's always seemed to work, so



      I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



      you a unique position from the point of interest (the station).







      This time, I have some data output that is expressed in AZ and DIST. I



      assumed (perhaps naively?) that this would also allow successful



      spatial grouping, but it turns out not to be the case; one can have



      the same forward azimuth for several different locations.







      Here's the example where I found the issue:



      Station lat, lon: 82.5033, -62.35 (ALE)







      Event1: lat 80.205, lon -1.091



      Event 2: lat 86.876, lon 54.373



      The two events are almost 1000km apart.







      Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



      Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







      I'd like to know more about the calculation, in particular because (i)



      I'd like to reassure myself that BAZ/distance does indeed give a



      unique location, and (ii) I don't understand how AZ is calculated. One



      website I looked at (Matlab) talks about rhumb lines vs great-circles;



      is that the issue? Will an AZ/distance pair always be non-unique?



      (Also most websites define back-azimuth as simply 180 opposite to



      azimuth, but I guess they aren't talking about spherical geometry



      then...)







      Thanks for any insight.







      Fiona Darbyshire.







      Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







      _______________________________________________



      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:







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







      _______________________________________________



      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 97, Issue 2



      ***************************************







      Get your own FREE website, FREE domain & FREE mobile app with Company
      email.

      Know More >


      <2009-11-14-1304-50S.KGD___003_KGD___HH_Z_SAC>_______________________________________________


      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:



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



      _______________________________________________

      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 97, Issue 7

      ***************************************



      http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?>
      Get your own *FREE* website, *FREE* domain & *FREE* mobile app with
      Company email.
      *Know More >*http://track.rediff.com/click?url=___http://businessemail.rediff.com/company-email-hosting-services?sc_cid=sign-1-10-13___&cmp=host&lnk=sign-1-10-13&nsrv1=host
      _______________________________________________
      sac-help mailing list
      sac-help<at>iris.washington.edu
      http://www.iris.washington.edu/mailman/listinfo/sac-help



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



      • Milton Plasencia
        2014-02-18 16:20:04
        Hi Raja,

        See below

        Milton P. Plasencia Linares
        mpplasencia<at>gmail.com



        On Feb 18, 2014, at 4:39 AM, Raja Babu <rjbb013<at>gmail.com> wrote:

        Hi all,

        normalization constant in radian = (normalization constant in Hz)2π^ (numberof poles−number of zeroes)
        Correct, and it is a constant (not change)

        for displacement normalization constant in above formula, number of zeroes will include one additional zero than velocity normalization constant.?????????
        No.

        You must add 1 zero to the P&Z file to convert velocity to displacement, so you multiply by omega(w).

        I have confusion. can anyone clear it
        For a bit of clarity and details of units, see:
        http://vps.isti.com/trac/ew/wiki/localmag

        any help will be appreciated..

        and thanks for help...
        Cheers,
        Milton



        On Wed, Nov 27, 2013 at 2:17 PM, Milton Plasencia <mpplasencia<at>gmail.com> wrote:
        Hi,

        You can consider using the PDCC software by IRIS.
        http://www.iris.edu/dms/nodes/dmc/software/downloads/pdcc/

        For Guralp sensor you must have the TTL value (Tap TableLookup)
        for apply the correct FIR filter combination). look here
        http://www.guralp.com/determining-fir-filter-information-from-gcf-data/

        HTH,

        m.

        * * * * * * * * * * * * * * * * * * * * * * *
        Milton P. PLASENCIA LINARES

        Centro di Ricerche Sismologiche (CRS)
        OGS - Istituto Nazionale di Oceanografia
        e di Geofisica Sperimentale

        Borgo Grotta Gigante 42/C
        (34010) Sgonico - Trieste - Italia
        Tel: +39 040 2140 141 (Udine)
        Tel: +39 040 2140 256 (Trieste)
        Cel.: +39 331 6481 935
        Fax: +39-040-327307

        E-mail: mplasencia<at>inogs.it

        ASAIN (Antarctic Seismographic Argentinean Italian Network)
        *********************************










        On Nov 27, 2013, at 3:03 AM, shekar <m_shekar2k2<at>rediffmail.com> wrote:

        Hai all

        Any body has the instrument response file(PZ file) for the following combination

        1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
        2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

        Thank you very much.

        shekar


        On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



        Message: 1

        Date: Tue, 26 Nov 2013 06:45:34 -0300

        From: Milton Plasencia

        Subject: Re: [SAC-HELP] Insert hedder

        To: shekar

        Cc: sac-help<at>iris.washington.edu

        Message-ID:

        Content-Type: text/plain; charset="us-ascii"



        Hi,

        Use

        ch O GMT 2009 318 13 03 34 300



        where 318 is the Julian day corresponding to month 11 and day 14.



        HTH,



        Milton P. Plasencia Linares

        mpplasencia<at>gmail.com







        On Nov 26, 2013, at 6:34 AM, shekar wrote:



        Dear all,



        I have a SAC file converted from SEISAN.



        How to insert the origin time(2009-11-14-13:03:34.300) of the event in sac header?



        Thanks in advance



        Shekar





        On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



        2. Remove Response Instrument (amik Amik)



        3. Re: Remove Response Instrument (Arthur Snoke)



        4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



        5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



        6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







        Message: 1



        Date: Wed, 13 Nov 2013 09:03:30 +0000



        From: George Helffrich



        Subject: Re: [SAC-HELP] sort and plot commands don't play



        To: Frederik Tilmann



        Cc: sac-help<at>iris.washington.edu



        Message-ID:



        Content-Type: text/plain; charset=us-ascii







        Dear All -







        Plotting the traces in a record section in SSS displays traces in increasing distance order.







        On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







        I have a problem in that the plot commands (p1,ppk) do not show the files in the same order as they are listed in memory after a sort.







        What I am trying to do is to display a north- to south record section of noise correlation functions with respect to a target station. The coordinates of both stations are place in evla and stla, but it is not well defined in which variable a given station ends up in (evla, or stla).







        So to plot all the functions relative to J10S I use







        r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



        sort evla ascend stla ascend



        lh evla stla







        Everything looks good, stations are in the right order in the lh list.







        But if I then do



        p1



        or ppk



        they come out in the wrong order, plotting the J03-J10 after J10-J11 (correct would have been J09-J10 as in the lh listing). I attach a screenshot demonstrating this effect.



        I tried



        sort commit evla ascend stla ascend



        but this this did not make a difference.







        This is on 101.5c



        SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











        Any ideas would be appreciated







        Regards



        Frederik



















        Output of lh command







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac - 1



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







        evla = 5.253600e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac - 2



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







        evla = 5.270155e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac - 3



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







        evla = 5.287174e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac - 4



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







        evla = 5.305374e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac - 5



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







        evla = 5.320426e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac - 6



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







        evla = 5.339244e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac - 7



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







        evla = 5.355617e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac - 8



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







        evla = 5.366659e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac - 9



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







        evla = 5.390106e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac - 10



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







        evla = 5.411258e+01



        Waiting







        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac - 11



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







        evla = 5.427260e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac - 12



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







        evla = 5.444247e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac - 13



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







        evla = 5.460217e+01



        stla = 5.481430e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac - 14



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







        evla = 5.481430e+01



        stla = 5.640510e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac - 15



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







        evla = 5.481430e+01



        stla = 5.617126e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac - 16



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







        evla = 5.481430e+01



        stla = 5.603592e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac - 17



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







        evla = 5.481430e+01



        stla = 5.577234e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac - 18



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







        evla = 5.481430e+01



        stla = 5.530829e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac - 19



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







        evla = 5.481430e+01



        stla = 5.512811e+01







        FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac - 20



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







        evla = 5.481430e+01



        Waiting







        stla = 5.497167e+01



        SAC>



        _______________________________________________



        sac-help mailing list



        sac-help<at>iris.washington.edu



        http://www.iris.washington.edu/mailman/listinfo/sac-help















        George Helffrich



        george.helffrich<at>bris.ac.uk



















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







        Message: 2



        Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



        From: amik Amik



        Subject: [SAC-HELP] Remove Response Instrument



        To: "sac-help<at>iris.washington.edu"



        Message-ID:



        <1384675653.16429.YahooMailNeo<at>web190801.mail.sg3.yahoo.com>



        Content-Type: text/plain; charset="iso-8859-1"







        I am working on a local magnitude research. I want to remove response instrument



        from the data and simulate the response of Wood Anderson seismograph. What should I do?



        I have a few questions :



        1. Which one the true procedure and command in SAC :



        ?? a) SAC> trans from evalresp to none



        ?? ? ?and then



        ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







        ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







        ?? c) SAC> trans to wa







        ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







        ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







        ?? f) SAC> trans from evalresp fname RESP.SATION to wa







        ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq 0.005 0.01 1.0 2.0







        2. What the different utilizing of Response File and Polezero File?



        ?? I think polezero already included in responfile, is it right ?







        May you help me? Thank You very much











        Sincerely Yours



        Ami



        -------------- next part --------------



        An HTML attachment was scrubbed...



        URL:







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







        Message: 3



        Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



        From: Arthur Snoke



        Subject: Re: [SAC-HELP] Remove Response Instrument



        To: amik Amik



        Cc: "sac-help<at>iris.washington.edu"



        Message-ID:



        Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







        I think the best way for you to understand what the various options do is



        to try them andsave the results. The TRANSFER help file at



        has answers



        to your explicit quesitons. If anything is not clear, let me know.







        On Sun, 17 Nov 2013, amik Amik wrote:







        I am working on a local magnitude research. I want to remove response



        instrument



        from the data and simulate the response of Wood Anderson seismograph. What



        should I do?



        I have a few questions :



        1. Which one the true procedure and command in SAC :



        ?? a) SAC> trans from evalresp to none



        ?? ? ?and then



        ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







        ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







        ?? c) SAC> trans to wa







        ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







        ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







        ?? f) SAC> trans from evalresp fname RESP.SATION to wa







        ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq



        0.005 0.01 1.0 2.0







        2. What the different utilizing of Response File and Polezero File?



        ?? I think polezero already included in responfile, is it right ?







        May you help me? Thank You very much











        Sincerely Yours



        Ami















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







        Message: 4



        Date: Sun, 17 Nov 2013 14:29:15 +0100



        From: "Milton P. Plasencia Linares"



        Subject: Re: [SAC-HELP] Remove Response Instrument



        To: "amik Amik"



        Cc: "sac-help<at>iris.washington.edu"



        Message-ID: <20131117142915.910179rx5iy7g6nf<at>webmail.ogs.trieste.it>



        Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



        format="flowed"







        look below,







        **********************************



        Milton P. PLASENCIA LINARES







        Dipartimento Centro di Ricerche Sismologiche



        Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







        Borgo Grotta Gigante 42/C



        (34010) Sgonico - TRIESTE - ITALIA



        Tel: +39-040-2140 256/141



        Fax: +39-040-327307







        E-mail: mplasencia<at>inogs.it







        ASAIN (Antarctic Seismographic Argentinean Italian Network)



        *********************************











        Quoting "amik Amik" :







        I am working on a local magnitude research. I want to remove



        response instrument



        from the data and simulate the response of Wood Anderson



        seismograph. What should I do?



        I have a few questions :



        1. Which one the true procedure and command in SAC :



        ?? a) SAC> trans from evalresp to none



        ?? ? ?and then



        ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







        ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







        ?? c) SAC> trans to wa







        ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



        Option d) This is fine if you use the polezero station file.







        ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







        ?? f) SAC> trans from evalresp fname RESP.SATION to wa







        ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



        freq 0.005 0.01 1.0 2.0







        2. What the different utilizing of Response File and Polezero File?



        Response file include the filter (FIR) coefficients







        ?? I think polezero already included in responfile, is it right ?



        Yes.







        May you help me? Thank You very much







        My suggest ius use:







        transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















        Sincerely Yours



        Ami







        Cheers,







        Milton















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



        This message was sent using IMP, the Internet Messaging Program.



















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







        Message: 5



        Date: Mon, 18 Nov 2013 08:01:33 -0500



        From: Fiona Darbyshire



        Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



        To: sac-help<at>iris.washington.edu



        Message-ID:







        Content-Type: text/plain; charset=ISO-8859-1







        Hello,







        Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







        Up to now, e.g. when binning receiver functions, I have used BAZ and



        GCARC to group them together, and that's always seemed to work, so



        I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



        you a unique position from the point of interest (the station).







        This time, I have some data output that is expressed in AZ and DIST. I



        assumed (perhaps naively?) that this would also allow successful



        spatial grouping, but it turns out not to be the case; one can have



        the same forward azimuth for several different locations.







        Here's the example where I found the issue:



        Station lat, lon: 82.5033, -62.35 (ALE)







        Event1: lat 80.205, lon -1.091



        Event 2: lat 86.876, lon 54.373



        The two events are almost 1000km apart.







        Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



        Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







        I'd like to know more about the calculation, in particular because (i)



        I'd like to reassure myself that BAZ/distance does indeed give a



        unique location, and (ii) I don't understand how AZ is calculated. One



        website I looked at (Matlab) talks about rhumb lines vs great-circles;



        is that the issue? Will an AZ/distance pair always be non-unique?



        (Also most websites define back-azimuth as simply 180 opposite to



        azimuth, but I guess they aren't talking about spherical geometry



        then...)







        Thanks for any insight.







        Fiona Darbyshire.







        Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







        Message: 6



        Date: Mon, 18 Nov 2013 11:38:44 -0500



        From: Brian Savage



        Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



        To: SAC HELP



        Message-ID: <8D3AE955-FCD5-48D4-8643-4E4E25115DEF<at>uri.edu>



        Content-Type: text/plain; charset="windows-1252"







        Fiona,







        The implementation of the computation of Great Circle Distance (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's formula found in Geodesy by R.Bomford (1980). It assumes a Radius of 6378.160 km and a flattening of 1/298.2466081. I can send you the appropriate pages from R. Bomford (1980) if you would like.







        The azimuth and back-azimuth will generally not be offset by180 degrees, and some pairs of point will have "odd-looking" (not event close to 180 degrees) azimuth/back azimuth pairs.







        Note: There are more robust methods of calculating these quantities and the current implementation in SAC as well as many other implementations, up until recently [1], have trouble with particular pairs of points. For over ~ 99% of point pairs SAC's implementation produces the correct result, but that last 1% can be incorrect (sometimes badly).







        I hope this helps.







        Brian Savage







        [1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z (open access). Addenda.











        On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







        Hello,







        Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







        Up to now, e.g. when binning receiver functions, I have used BAZ and



        GCARC to group them together, and that's always seemed to work, so



        I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



        you a unique position from the point of interest (the station).







        This time, I have some data output that is expressed in AZ and DIST. I



        assumed (perhaps naively?) that this would also allow successful



        spatial grouping, but it turns out not to be the case; one can have



        the same forward azimuth for several different locations.







        Here's the example where I found the issue:



        Station lat, lon: 82.5033, -62.35 (ALE)







        Event1: lat 80.205, lon -1.091



        Event 2: lat 86.876, lon 54.373



        The two events are almost 1000km apart.







        Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



        Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







        I'd like to know more about the calculation, in particular because (i)



        I'd like to reassure myself that BAZ/distance does indeed give a



        unique location, and (ii) I don't understand how AZ is calculated. One



        website I looked at (Matlab) talks about rhumb lines vs great-circles;



        is that the issue? Will an AZ/distance pair always be non-unique?



        (Also most websites define back-azimuth as simply 180 opposite to



        azimuth, but I guess they aren't talking about spherical geometry



        then...)







        Thanks for any insight.







        Fiona Darbyshire.







        Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







        _______________________________________________



        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:







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







        _______________________________________________



        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 97, Issue 2



        ***************************************







        Get your own FREE website, FREE domain & FREE mobile app with Company email.

        Know More >

        <2009-11-14-1304-50S.KGD___003_KGD___HH_Z_SAC>_______________________________________________


        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:



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



        _______________________________________________

        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 97, Issue 7

        ***************************************



        Get your own FREE website, FREE domain & FREE mobile app with Company email.
        Know More >
        _______________________________________________
        sac-help mailing list
        sac-help<at>iris.washington.edu
        http://www.iris.washington.edu/mailman/listinfo/sac-help


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


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


        • Satish_iitb
          2014-02-18 18:49:03
          Hi Everybody,

          some discussion is going on pole zeros and response

          I have doubt on using of "SAC transfer",
          *Using rdseed:*
          In general pole zeros in displacement (we don't have to add extra zeros
          because rdseed take-care of this) and RESP in velocity
          *On SAC: we do*
          r sacfilename
          rmean
          taper
          transfer from polezero subtype PZfile to *vel* freqlimits 0.001 0.002 1 10
          (velocity pole zeros, It's okay because seismogram in vel) I think we have
          to use Velocity pole zeros...
          *OR*
          transfer from polezero subtype PZfile to *none* freqlimits 0.001 0.002 1
          10 (using displacement pole zeros not sure why?)

          *Since, we have seismogram in velocity so our pole zeros in velocity why
          displacement?*

          *Thanks for looking the question and giving Important time...*




          On Tue, Feb 18, 2014 at 8:20 AM, Milton Plasencia <mpplasencia<at>gmail.com>wrote:

          Hi Raja,

          See below

          Milton P. Plasencia Linares
          mpplasencia<at>gmail.com



          On Feb 18, 2014, at 4:39 AM, Raja Babu <rjbb013<at>gmail.com> wrote:

          Hi all,

          normalization constant in radian = (normalization constant in Hz)2π^
          (numberof poles−number of zeroes)

          Correct, and it is a constant (not change)


          for displacement normalization constant in above formula, number of zeroes
          will include one additional zero than velocity normalization
          constant.?????????

          No.

          You must add 1 zero to the P&Z file to convert velocity to displacement,
          so you multiply by omega(w).


          I have confusion. can anyone clear it

          For a bit of clarity and details of units, see:
          http://vps.isti.com/trac/ew/wiki/localmag


          any help will be appreciated..

          and thanks for help...

          Cheers,
          Milton



          On Wed, Nov 27, 2013 at 2:17 PM, Milton Plasencia <mpplasencia<at>gmail.com>wrote:

          Hi,

          You can consider using the PDCC software by IRIS.
          http://www.iris.edu/dms/nodes/dmc/software/downloads/pdcc/

          For Guralp sensor you must have the TTL value (Tap TableLookup)
          for apply the correct FIR filter combination). look here
          http://www.guralp.com/determining-fir-filter-information-from-gcf-data/

          HTH,

          m.

          * * * * * * * * * * * * * * * * * * * * * * *
          Milton P. PLASENCIA LINARES

          Centro di Ricerche Sismologiche (CRS)
          OGS - Istituto Nazionale di Oceanografia
          e di Geofisica Sperimentale

          Borgo Grotta Gigante 42/C
          (34010) Sgonico - Trieste - Italia
          Tel: +39 040 2140 141 (Udine)
          Tel: +39 040 2140 256 (Trieste)
          Cel.: +39 331 6481 935
          Fax: +39-040-327307

          E-mail: mplasencia<at>inogs.it

          ASAIN (Antarctic Seismographic Argentinean Italian Network)
          *********************************










          On Nov 27, 2013, at 3:03 AM, shekar <m_shekar2k2<at>rediffmail.com> wrote:

          Hai all

          Any body has the instrument response file(PZ file) for the following
          combination

          1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
          2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

          Thank you very much.

          shekar


          On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



          Message: 1

          Date: Tue, 26 Nov 2013 06:45:34 -0300

          From: Milton Plasencia

          Subject: Re: [SAC-HELP] Insert hedder

          To: shekar

          Cc: sac-help<at>iris.washington.edu

          Message-ID:

          Content-Type: text/plain; charset="us-ascii"



          Hi,

          Use

          ch O GMT 2009 318 13 03 34 300



          where 318 is the Julian day corresponding to month 11 and day 14.



          HTH,



          Milton P. Plasencia Linares

          mpplasencia<at>gmail.com







          On Nov 26, 2013, at 6:34 AM, shekar wrote:



          Dear all,



          I have a SAC file converted from SEISAN.



          How to insert the origin time(2009-11-14-13:03:34.300) of the event in
          sac header?



          Thanks in advance



          Shekar





          On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



          2. Remove Response Instrument (amik Amik)



          3. Re: Remove Response Instrument (Arthur Snoke)



          4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



          5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



          6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







          Message: 1



          Date: Wed, 13 Nov 2013 09:03:30 +0000



          From: George Helffrich



          Subject: Re: [SAC-HELP] sort and plot commands don't play



          To: Frederik Tilmann



          Cc: sac-help<at>iris.washington.edu



          Message-ID:



          Content-Type: text/plain; charset=us-ascii







          Dear All -







          Plotting the traces in a record section in SSS displays traces in
          increasing distance order.







          On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







          I have a problem in that the plot commands (p1,ppk) do not show the
          files in the same order as they are listed in memory after a sort.







          What I am trying to do is to display a north- to south record section
          of noise correlation functions with respect to a target station. The
          coordinates of both stations are place in evla and stla, but it is not well
          defined in which variable a given station ends up in (evla, or stla).







          So to plot all the functions relative to J10S I use







          r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac
          sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



          sort evla ascend stla ascend



          lh evla stla







          Everything looks good, stations are in the right order in the lh list.







          But if I then do



          p1



          or ppk



          they come out in the wrong order, plotting the J03-J10 after J10-J11
          (correct would have been J09-J10 as in the lh listing). I attach a
          screenshot demonstrating this effect.



          I tried



          sort commit evla ascend stla ascend



          but this this did not make a difference.







          This is on 101.5c



          SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











          Any ideas would be appreciated







          Regards



          Frederik



















          Output of lh command







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac -
          1



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







          evla = 5.253600e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac -
          2



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







          evla = 5.270155e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac -
          3



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







          evla = 5.287174e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac -
          4



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







          evla = 5.305374e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac -
          5



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







          evla = 5.320426e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac -
          6



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







          evla = 5.339244e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac -
          7



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







          evla = 5.355617e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac -
          8



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







          evla = 5.366659e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac -
          9



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







          evla = 5.390106e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac -
          10



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







          evla = 5.411258e+01



          Waiting







          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac -
          11



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







          evla = 5.427260e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac -
          12



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







          evla = 5.444247e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac -
          13



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







          evla = 5.460217e+01



          stla = 5.481430e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac -
          14



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







          evla = 5.481430e+01



          stla = 5.640510e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac -
          15



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







          evla = 5.481430e+01



          stla = 5.617126e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac -
          16



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







          evla = 5.481430e+01



          stla = 5.603592e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac -
          17



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







          evla = 5.481430e+01



          stla = 5.577234e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac -
          18



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







          evla = 5.481430e+01



          stla = 5.530829e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac -
          19



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







          evla = 5.481430e+01



          stla = 5.512811e+01







          FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac -
          20



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







          evla = 5.481430e+01



          Waiting







          stla = 5.497167e+01



          SAC>



          _______________________________________________



          sac-help mailing list



          sac-help<at>iris.washington.edu



          http://www.iris.washington.edu/mailman/listinfo/sac-help















          George Helffrich



          george.helffrich<at>bris.ac.uk



















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







          Message: 2



          Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



          From: amik Amik



          Subject: [SAC-HELP] Remove Response Instrument



          To: "sac-help<at>iris.washington.edu"



          Message-ID:



          <1384675653.16429.YahooMailNeo<at>web190801.mail.sg3.yahoo.com>



          Content-Type: text/plain; charset="iso-8859-1"







          I am working on a local magnitude research. I want to remove response
          instrument



          from the data and simulate the response of Wood Anderson seismograph.
          What should I do?



          I have a few questions :



          1. Which one the true procedure and command in SAC :



          ?? a) SAC> trans from evalresp to none



          ?? ? ?and then



          ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







          ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







          ?? c) SAC> trans to wa







          ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







          ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to
          none







          ?? f) SAC> trans from evalresp fname RESP.SATION to wa







          ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"
          freq 0.005 0.01 1.0 2.0







          2. What the different utilizing of Response File and Polezero File?



          ?? I think polezero already included in responfile, is it right ?







          May you help me? Thank You very much











          Sincerely Yours



          Ami



          -------------- next part --------------



          An HTML attachment was scrubbed...



          URL:







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







          Message: 3



          Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



          From: Arthur Snoke



          Subject: Re: [SAC-HELP] Remove Response Instrument



          To: amik Amik



          Cc: "sac-help<at>iris.washington.edu"



          Message-ID:



          Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







          I think the best way for you to understand what the various options do
          is



          to try them andsave the results. The TRANSFER help file at



          has answers



          to your explicit quesitons. If anything is not clear, let me know.







          On Sun, 17 Nov 2013, amik Amik wrote:







          I am working on a local magnitude research. I want to remove response



          instrument



          from the data and simulate the response of Wood Anderson seismograph.
          What



          should I do?



          I have a few questions :



          1. Which one the true procedure and command in SAC :



          ?? a) SAC> trans from evalresp to none



          ?? ? ?and then



          ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







          ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







          ?? c) SAC> trans to wa







          ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







          ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa
          to none







          ?? f) SAC> trans from evalresp fname RESP.SATION to wa







          ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"
          freq



          0.005 0.01 1.0 2.0







          2. What the different utilizing of Response File and Polezero File?



          ?? I think polezero already included in responfile, is it right ?







          May you help me? Thank You very much











          Sincerely Yours



          Ami















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







          Message: 4



          Date: Sun, 17 Nov 2013 14:29:15 +0100



          From: "Milton P. Plasencia Linares"



          Subject: Re: [SAC-HELP] Remove Response Instrument



          To: "amik Amik"



          Cc: "sac-help<at>iris.washington.edu"



          Message-ID: <20131117142915.910179rx5iy7g6nf<at>webmail.ogs.trieste.it>



          Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



          format="flowed"







          look below,







          **********************************



          Milton P. PLASENCIA LINARES







          Dipartimento Centro di Ricerche Sismologiche



          Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







          Borgo Grotta Gigante 42/C



          (34010) Sgonico - TRIESTE - ITALIA



          Tel: +39-040-2140 256/141



          Fax: +39-040-327307







          E-mail: mplasencia<at>inogs.it







          ASAIN (Antarctic Seismographic Argentinean Italian Network)



          *********************************











          Quoting "amik Amik" :







          I am working on a local magnitude research. I want to remove



          response instrument



          from the data and simulate the response of Wood Anderson



          seismograph. What should I do?



          I have a few questions :



          1. Which one the true procedure and command in SAC :



          ?? a) SAC> trans from evalresp to none



          ?? ? ?and then



          ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







          ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







          ?? c) SAC> trans to wa







          ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



          Option d) This is fine if you use the polezero station file.







          ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa
          to none







          ?? f) SAC> trans from evalresp fname RESP.SATION to wa







          ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



          freq 0.005 0.01 1.0 2.0







          2. What the different utilizing of Response File and Polezero File?



          Response file include the filter (FIR) coefficients







          ?? I think polezero already included in responfile, is it right ?



          Yes.







          May you help me? Thank You very much







          My suggest ius use:







          transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















          Sincerely Yours



          Ami







          Cheers,







          Milton















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



          This message was sent using IMP, the Internet Messaging Program.



















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







          Message: 5



          Date: Mon, 18 Nov 2013 08:01:33 -0500



          From: Fiona Darbyshire



          Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



          To: sac-help<at>iris.washington.edu



          Message-ID:







          Content-Type: text/plain; charset=ISO-8859-1







          Hello,







          Could someone tell me a bit more about how AZ and BAZ are calculated in
          SAC?







          Up to now, e.g. when binning receiver functions, I have used BAZ and



          GCARC to group them together, and that's always seemed to work, so



          I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



          you a unique position from the point of interest (the station).







          This time, I have some data output that is expressed in AZ and DIST. I



          assumed (perhaps naively?) that this would also allow successful



          spatial grouping, but it turns out not to be the case; one can have



          the same forward azimuth for several different locations.







          Here's the example where I found the issue:



          Station lat, lon: 82.5033, -62.35 (ALE)







          Event1: lat 80.205, lon -1.091



          Event 2: lat 86.876, lon 54.373



          The two events are almost 1000km apart.







          Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



          Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







          I'd like to know more about the calculation, in particular because (i)



          I'd like to reassure myself that BAZ/distance does indeed give a



          unique location, and (ii) I don't understand how AZ is calculated. One



          website I looked at (Matlab) talks about rhumb lines vs great-circles;



          is that the issue? Will an AZ/distance pair always be non-unique?



          (Also most websites define back-azimuth as simply 180 opposite to



          azimuth, but I guess they aren't talking about spherical geometry



          then...)







          Thanks for any insight.







          Fiona Darbyshire.







          Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







          Message: 6



          Date: Mon, 18 Nov 2013 11:38:44 -0500



          From: Brian Savage



          Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



          To: SAC HELP



          Message-ID: <8D3AE955-FCD5-48D4-8643-4E4E25115DEF<at>uri.edu>



          Content-Type: text/plain; charset="windows-1252"







          Fiona,







          The implementation of the computation of Great Circle Distance
          (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's
          formula found in Geodesy by R.Bomford (1980). It assumes a Radius of
          6378.160 km and a flattening of 1/298.2466081. I can send you the
          appropriate pages from R. Bomford (1980) if you would like.







          The azimuth and back-azimuth will generally not be offset by180
          degrees, and some pairs of point will have "odd-looking" (not event close
          to 180 degrees) azimuth/back azimuth pairs.







          Note: There are more robust methods of calculating these quantities and
          the current implementation in SAC as well as many other implementations, up
          until recently [1], have trouble with particular pairs of points. For over
          ~ 99% of point pairs SAC's implementation produces the correct result, but
          that last 1% can be incorrect (sometimes badly).







          I hope this helps.







          Brian Savage







          [1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of
          Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K.
          doi:10.1007/s00190-012-0578-z (open access). Addenda.











          On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







          Hello,







          Could someone tell me a bit more about how AZ and BAZ are calculated
          in SAC?







          Up to now, e.g. when binning receiver functions, I have used BAZ and



          GCARC to group them together, and that's always seemed to work, so



          I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



          you a unique position from the point of interest (the station).







          This time, I have some data output that is expressed in AZ and DIST. I



          assumed (perhaps naively?) that this would also allow successful



          spatial grouping, but it turns out not to be the case; one can have



          the same forward azimuth for several different locations.







          Here's the example where I found the issue:



          Station lat, lon: 82.5033, -62.35 (ALE)







          Event1: lat 80.205, lon -1.091



          Event 2: lat 86.876, lon 54.373



          The two events are almost 1000km apart.







          Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



          Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







          I'd like to know more about the calculation, in particular because (i)



          I'd like to reassure myself that BAZ/distance does indeed give a



          unique location, and (ii) I don't understand how AZ is calculated. One



          website I looked at (Matlab) talks about rhumb lines vs great-circles;



          is that the issue? Will an AZ/distance pair always be non-unique?



          (Also most websites define back-azimuth as simply 180 opposite to



          azimuth, but I guess they aren't talking about spherical geometry



          then...)







          Thanks for any insight.







          Fiona Darbyshire.







          Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







          _______________________________________________



          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:







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







          _______________________________________________



          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 97, Issue 2



          ***************************************







          Get your own FREE website, FREE domain & FREE mobile app with Company
          email.

          Know More >


          <2009-11-14-1304-50S.KGD___003_KGD___HH_Z_SAC>_______________________________________________


          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:



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



          _______________________________________________

          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 97, Issue 7

          ***************************************



          http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?>
          Get your own *FREE* website, *FREE* domain & *FREE* mobile app with
          Company email.
          *Know More >*http://track.rediff.com/click?url=___http://businessemail.rediff.com/company-email-hosting-services?sc_cid=sign-1-10-13___&cmp=host&lnk=sign-1-10-13&nsrv1=host
          _______________________________________________
          sac-help mailing list
          sac-help<at>iris.washington.edu
          http://www.iris.washington.edu/mailman/listinfo/sac-help



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


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



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




          --
          ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          *Satish Maurya*


          * Ph.D. studentDept. of Seismology*
          *IPGP, France*
          ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
          ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

          • Milton Plasencia
            2014-02-18 20:03:10
            See below,

            Milton P. Plasencia Linares
            mpplasencia<at>gmail.com



            On Feb 18, 2014, at 10:49 AM, satish maurya <maurya2278satish<at>gmail.com> wrote:

            Hi Everybody,

            some discussion is going on pole zeros and response

            I have doubt on using of "SAC transfer",
            Using rdseed:
            In general pole zeros in displacement (we don't have to add extra zeros because rdseed take-care of this)
            Correct, if you download seed volumes and convert to sac with P&Z output using rdseed, it add the necessary Zeros and you
            P&Z files give you displacement ( add 2 Zeros if you have acceleration and only one if you have velocity),
            (after applied transfer command in sac), displacement is used to calculate the magnitude.
            rdseed use meters instead sac use nm.

            and RESP in velocity
            On SAC: we do
            r sacfilename
            rmean
            taper
            transfer from polezero subtype PZfile to vel freqlimits 0.001 0.002 1 10 (velocity pole zeros, It's okay because seismogram in vel) I think we have to use Velocity pole zeros...
            Ok, but it depend of your application
            OR

            This is correct to remove the influence instrument of the signal.
            transfer from polezero subtype PZfile to none freqlimits 0.001 0.002 1 10 (using displacement pole zeros not sure why?)

            Since, we have seismogram in velocity so our pole zeros in velocity why displacement?
            For calculate magnitude.

            Thanks for looking the question and giving Important time...
            HTH,

            Milton




            On Tue, Feb 18, 2014 at 8:20 AM, Milton Plasencia <mpplasencia<at>gmail.com> wrote:
            Hi Raja,

            See below

            Milton P. Plasencia Linares
            mpplasencia<at>gmail.com



            On Feb 18, 2014, at 4:39 AM, Raja Babu <rjbb013<at>gmail.com> wrote:

            Hi all,

            normalization constant in radian = (normalization constant in Hz)2π^ (numberof poles−number of zeroes)
            Correct, and it is a constant (not change)


            for displacement normalization constant in above formula, number of zeroes will include one additional zero than velocity normalization constant.?????????
            No.

            You must add 1 zero to the P&Z file to convert velocity to displacement, so you multiply by omega(w).


            I have confusion. can anyone clear it
            For a bit of clarity and details of units, see:
            http://vps.isti.com/trac/ew/wiki/localmag


            any help will be appreciated..

            and thanks for help...
            Cheers,
            Milton



            On Wed, Nov 27, 2013 at 2:17 PM, Milton Plasencia <mpplasencia<at>gmail.com> wrote:
            Hi,

            You can consider using the PDCC software by IRIS.
            http://www.iris.edu/dms/nodes/dmc/software/downloads/pdcc/

            For Guralp sensor you must have the TTL value (Tap TableLookup)
            for apply the correct FIR filter combination). look here
            http://www.guralp.com/determining-fir-filter-information-from-gcf-data/

            HTH,

            m.

            * * * * * * * * * * * * * * * * * * * * * * *
            Milton P. PLASENCIA LINARES

            Centro di Ricerche Sismologiche (CRS)
            OGS - Istituto Nazionale di Oceanografia
            e di Geofisica Sperimentale

            Borgo Grotta Gigante 42/C
            (34010) Sgonico - Trieste - Italia
            Tel: +39 040 2140 141 (Udine)
            Tel: +39 040 2140 256 (Trieste)
            Cel.: +39 331 6481 935
            Fax: +39-040-327307

            E-mail: mplasencia<at>inogs.it

            ASAIN (Antarctic Seismographic Argentinean Italian Network)
            *********************************










            On Nov 27, 2013, at 3:03 AM, shekar <m_shekar2k2<at>rediffmail.com> wrote:

            Hai all

            Any body has the instrument response file(PZ file) for the following combination

            1. Taurus digitizer+ Trillium 240 seismometer(240 sec, 1500v/m/s)
            2. Reftek 72A+ CMG 40T(30 sec, 800 v/m/s)

            Thank you very much.

            shekar


            On Tue, 26 Nov 2013 15:21:35 +0530 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: Insert hedder (Milton Plasencia)





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



            Message: 1

            Date: Tue, 26 Nov 2013 06:45:34 -0300

            From: Milton Plasencia

            Subject: Re: [SAC-HELP] Insert hedder

            To: shekar

            Cc: sac-help<at>iris.washington.edu

            Message-ID:

            Content-Type: text/plain; charset="us-ascii"



            Hi,

            Use

            ch O GMT 2009 318 13 03 34 300



            where 318 is the Julian day corresponding to month 11 and day 14.



            HTH,



            Milton P. Plasencia Linares

            mpplasencia<at>gmail.com







            On Nov 26, 2013, at 6:34 AM, shekar wrote:



            Dear all,



            I have a SAC file converted from SEISAN.



            How to insert the origin time(2009-11-14-13:03:34.300) of the event in sac header?



            Thanks in advance



            Shekar





            On Mon, 18 Nov 2013 22:14:49 +0530 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: sort and plot commands don't play (George Helffrich)



            2. Remove Response Instrument (amik Amik)



            3. Re: Remove Response Instrument (Arthur Snoke)



            4. Re: Remove Response Instrument (Milton P. Plasencia Linares)



            5. Calculation of AZ and BAZ in SAC (Fiona Darbyshire)



            6. Re: Calculation of AZ and BAZ in SAC (Brian Savage)











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







            Message: 1



            Date: Wed, 13 Nov 2013 09:03:30 +0000



            From: George Helffrich



            Subject: Re: [SAC-HELP] sort and plot commands don't play



            To: Frederik Tilmann



            Cc: sac-help<at>iris.washington.edu



            Message-ID:



            Content-Type: text/plain; charset=us-ascii







            Dear All -







            Plotting the traces in a record section in SSS displays traces in increasing distance order.







            On 13 Nov 2013, at 08:49, Frederik Tilmann wrote:







            I have a problem in that the plot commands (p1,ppk) do not show the files in the same order as they are listed in memory after a sort.







            What I am trying to do is to display a north- to south record section of noise correlation functions with respect to a target station. The coordinates of both stations are place in evla and stla, but it is not well defined in which variable a given station ends up in (evla, or stla).







            So to plot all the functions relative to J10S I use







            r sumtraces-nonorm-0.03-0.10/*ZWZW*J10?J[0-9][0-9]*_sum.sac sumtraces-nonorm-0.03-0.10/*ZWZW*J[0-9][0-9]?J10?*_sum.sac



            sort evla ascend stla ascend



            lh evla stla







            Everything looks good, stations are in the right order in the lh list.







            But if I then do



            p1



            or ppk



            they come out in the wrong order, plotting the J03-J10 after J10-J11 (correct would have been J09-J10 as in the lh listing). I attach a screenshot demonstrating this effect.



            I tried



            sort commit evla ascend stla ascend



            but this this did not make a difference.







            This is on 101.5c



            SEISMIC ANALYSIS CODE [02/01/2012 (Version 101.5c)]











            Any ideas would be appreciated







            Regards



            Frederik



















            Output of lh command







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ23B..HHZHHZ_sum.sac - 1



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







            evla = 5.253600e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ22B..HHZHHZ_sum.sac - 2



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







            evla = 5.270155e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ21B..HHZHHZ_sum.sac - 3



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







            evla = 5.287174e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ20B..HHZHHZ_sum.sac - 4



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







            evla = 5.305374e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ19B..HHZHHZ_sum.sac - 5



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







            evla = 5.320426e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ18B..HHZHHZ_sum.sac - 6



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







            evla = 5.339244e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ17B..HHZHHZ_sum.sac - 7



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







            evla = 5.355617e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ25B..HHZHHZ_sum.sac - 8



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







            evla = 5.366659e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ15B..HHZHHZ_sum.sac - 9



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







            evla = 5.390106e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ14S..HHZHHZ_sum.sac - 10



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







            evla = 5.411258e+01



            Waiting







            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ13B..HHZHHZ_sum.sac - 11



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







            evla = 5.427260e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ12B..HHZHHZ_sum.sac - 12



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







            evla = 5.444247e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J10BJ11S..HHZHHZ_sum.sac - 13



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







            evla = 5.460217e+01



            stla = 5.481430e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J09BJ10B..HHZHHZ_sum.sac - 14



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







            evla = 5.481430e+01



            stla = 5.640510e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J08SJ10B..HHZHHZ_sum.sac - 15



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







            evla = 5.481430e+01



            stla = 5.617126e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J07BJ10B..HHZHHZ_sum.sac - 16



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







            evla = 5.481430e+01



            stla = 5.603592e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J04BJ10B..HHZHHZ_sum.sac - 17



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







            evla = 5.481430e+01



            stla = 5.577234e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J03EJ10B..HHZHHZ_sum.sac - 18



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







            evla = 5.481430e+01



            stla = 5.530829e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J02SJ10B..HHZHHZ_sum.sac - 19



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







            evla = 5.481430e+01



            stla = 5.512811e+01







            FILE: sumtraces-nonorm-0.03-0.10/mat__ZWZW.J01EJ10B..HHZHHZ_sum.sac - 20



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







            evla = 5.481430e+01



            Waiting







            stla = 5.497167e+01



            SAC>



            _______________________________________________



            sac-help mailing list



            sac-help<at>iris.washington.edu



            http://www.iris.washington.edu/mailman/listinfo/sac-help















            George Helffrich



            george.helffrich<at>bris.ac.uk



















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







            Message: 2



            Date: Sun, 17 Nov 2013 16:07:33 +0800 (SGT)



            From: amik Amik



            Subject: [SAC-HELP] Remove Response Instrument



            To: "sac-help<at>iris.washington.edu"



            Message-ID:



            <1384675653.16429.YahooMailNeo<at>web190801.mail.sg3.yahoo.com>



            Content-Type: text/plain; charset="iso-8859-1"







            I am working on a local magnitude research. I want to remove response instrument



            from the data and simulate the response of Wood Anderson seismograph. What should I do?



            I have a few questions :



            1. Which one the true procedure and command in SAC :



            ?? a) SAC> trans from evalresp to none



            ?? ? ?and then



            ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







            ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







            ?? c) SAC> trans to wa







            ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







            ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







            ?? f) SAC> trans from evalresp fname RESP.SATION to wa







            ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq 0.005 0.01 1.0 2.0







            2. What the different utilizing of Response File and Polezero File?



            ?? I think polezero already included in responfile, is it right ?







            May you help me? Thank You very much











            Sincerely Yours



            Ami



            -------------- next part --------------



            An HTML attachment was scrubbed...



            URL:







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







            Message: 3



            Date: Sun, 17 Nov 2013 07:43:55 -0500 (EST)



            From: Arthur Snoke



            Subject: Re: [SAC-HELP] Remove Response Instrument



            To: amik Amik



            Cc: "sac-help<at>iris.washington.edu"



            Message-ID:



            Content-Type: text/plain; charset="iso-8859-1"; Format="flowed"







            I think the best way for you to understand what the various options do is



            to try them andsave the results. The TRANSFER help file at



            has answers



            to your explicit quesitons. If anything is not clear, let me know.







            On Sun, 17 Nov 2013, amik Amik wrote:







            I am working on a local magnitude research. I want to remove response



            instrument



            from the data and simulate the response of Wood Anderson seismograph. What



            should I do?



            I have a few questions :



            1. Which one the true procedure and command in SAC :



            ?? a) SAC> trans from evalresp to none



            ?? ? ?and then



            ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







            ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







            ?? c) SAC> trans to wa







            ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa







            ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







            ?? f) SAC> trans from evalresp fname RESP.SATION to wa







            ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz" freq



            0.005 0.01 1.0 2.0







            2. What the different utilizing of Response File and Polezero File?



            ?? I think polezero already included in responfile, is it right ?







            May you help me? Thank You very much











            Sincerely Yours



            Ami















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







            Message: 4



            Date: Sun, 17 Nov 2013 14:29:15 +0100



            From: "Milton P. Plasencia Linares"



            Subject: Re: [SAC-HELP] Remove Response Instrument



            To: "amik Amik"



            Cc: "sac-help<at>iris.washington.edu"



            Message-ID: <20131117142915.910179rx5iy7g6nf<at>webmail.ogs.trieste.it>



            Content-Type: text/plain;charset=ISO-8859-1;DelSp="Yes";



            format="flowed"







            look below,







            **********************************



            Milton P. PLASENCIA LINARES







            Dipartimento Centro di Ricerche Sismologiche



            Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS







            Borgo Grotta Gigante 42/C



            (34010) Sgonico - TRIESTE - ITALIA



            Tel: +39-040-2140 256/141



            Fax: +39-040-327307







            E-mail: mplasencia<at>inogs.it







            ASAIN (Antarctic Seismographic Argentinean Italian Network)



            *********************************











            Quoting "amik Amik" :







            I am working on a local magnitude research. I want to remove



            response instrument



            from the data and simulate the response of Wood Anderson



            seismograph. What should I do?



            I have a few questions :



            1. Which one the true procedure and command in SAC :



            ?? a) SAC> trans from evalresp to none



            ?? ? ?and then



            ?? ? ?SAC> trans from polezero subtype "file name.pz" to none







            ?? b) SAC> trans to polezero subtype "wood anderson file name.pz"







            ?? c) SAC> trans to wa







            ?? d) SAC> trans from polezero subtype "sation's file name.pz" to wa



            Option d) This is fine if you use the polezero station file.







            ?? e) SAC> trans from polezero subtype "sation's file name.pz" to wa to none







            ?? f) SAC> trans from evalresp fname RESP.SATION to wa







            ?? g) SAC> SAC> trans from polezero subtype "sation's file name.pz"



            freq 0.005 0.01 1.0 2.0







            2. What the different utilizing of Response File and Polezero File?



            Response file include the filter (FIR) coefficients







            ?? I think polezero already included in responfile, is it right ?



            Yes.







            May you help me? Thank You very much







            My suggest ius use:







            transfer from polezero subtype "PZ file" to wa freq 0.005 0.01 1.0 2.0















            Sincerely Yours



            Ami







            Cheers,







            Milton















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



            This message was sent using IMP, the Internet Messaging Program.



















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







            Message: 5



            Date: Mon, 18 Nov 2013 08:01:33 -0500



            From: Fiona Darbyshire



            Subject: [SAC-HELP] Calculation of AZ and BAZ in SAC



            To: sac-help<at>iris.washington.edu



            Message-ID:







            Content-Type: text/plain; charset=ISO-8859-1







            Hello,







            Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







            Up to now, e.g. when binning receiver functions, I have used BAZ and



            GCARC to group them together, and that's always seemed to work, so



            I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



            you a unique position from the point of interest (the station).







            This time, I have some data output that is expressed in AZ and DIST. I



            assumed (perhaps naively?) that this would also allow successful



            spatial grouping, but it turns out not to be the case; one can have



            the same forward azimuth for several different locations.







            Here's the example where I found the issue:



            Station lat, lon: 82.5033, -62.35 (ALE)







            Event1: lat 80.205, lon -1.091



            Event 2: lat 86.876, lon 54.373



            The two events are almost 1000km apart.







            Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



            Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







            I'd like to know more about the calculation, in particular because (i)



            I'd like to reassure myself that BAZ/distance does indeed give a



            unique location, and (ii) I don't understand how AZ is calculated. One



            website I looked at (Matlab) talks about rhumb lines vs great-circles;



            is that the issue? Will an AZ/distance pair always be non-unique?



            (Also most websites define back-azimuth as simply 180 opposite to



            azimuth, but I guess they aren't talking about spherical geometry



            then...)







            Thanks for any insight.







            Fiona Darbyshire.







            Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al















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







            Message: 6



            Date: Mon, 18 Nov 2013 11:38:44 -0500



            From: Brian Savage



            Subject: Re: [SAC-HELP] Calculation of AZ and BAZ in SAC



            To: SAC HELP



            Message-ID: <8D3AE955-FCD5-48D4-8643-4E4E25115DEF<at>uri.edu>



            Content-Type: text/plain; charset="windows-1252"







            Fiona,







            The implementation of the computation of Great Circle Distance (degrees), Azimuth, Back Azimuth, and Distance (km) is based off of Rudoe's formula found in Geodesy by R.Bomford (1980). It assumes a Radius of 6378.160 km and a flattening of 1/298.2466081. I can send you the appropriate pages from R. Bomford (1980) if you would like.







            The azimuth and back-azimuth will generally not be offset by180 degrees, and some pairs of point will have "odd-looking" (not event close to 180 degrees) azimuth/back azimuth pairs.







            Note: There are more robust methods of calculating these quantities and the current implementation in SAC as well as many other implementations, up until recently [1], have trouble with particular pairs of points. For over ~ 99% of point pairs SAC's implementation produces the correct result, but that last 1% can be incorrect (sometimes badly).







            I hope this helps.







            Brian Savage







            [1] Karney, C. F. F. (2013). "Algorithms for geodesics". Journal of Geodesy 87 (1): 43?42. arXiv:1109.4448. Bibcode:2013JGeod..87...43K. doi:10.1007/s00190-012-0578-z (open access). Addenda.











            On Nov 18, 2013, at 8:01 AM, Fiona Darbyshire wrote:







            Hello,







            Could someone tell me a bit more about how AZ and BAZ are calculated in SAC?







            Up to now, e.g. when binning receiver functions, I have used BAZ and



            GCARC to group them together, and that's always seemed to work, so



            I've always assumed that the BAZ/GCARC (or DIST, I imagine) pair gives



            you a unique position from the point of interest (the station).







            This time, I have some data output that is expressed in AZ and DIST. I



            assumed (perhaps naively?) that this would also allow successful



            spatial grouping, but it turns out not to be the case; one can have



            the same forward azimuth for several different locations.







            Here's the example where I found the issue:



            Station lat, lon: 82.5033, -62.35 (ALE)







            Event1: lat 80.205, lon -1.091



            Event 2: lat 86.876, lon 54.373



            The two events are almost 1000km apart.







            Event 1: AZ 313, BAZ 72, GCARC 9, DIST 1006



            Event 2: AZ 314, BAZ 17, GCARC 9, DIST 1042







            I'd like to know more about the calculation, in particular because (i)



            I'd like to reassure myself that BAZ/distance does indeed give a



            unique location, and (ii) I don't understand how AZ is calculated. One



            website I looked at (Matlab) talks about rhumb lines vs great-circles;



            is that the issue? Will an AZ/distance pair always be non-unique?



            (Also most websites define back-azimuth as simply 180 opposite to



            azimuth, but I guess they aren't talking about spherical geometry



            then...)







            Thanks for any insight.







            Fiona Darbyshire.







            Centre de recherche GEOTOP, Universit? du Qu?bec ? Montr?al







            _______________________________________________



            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:







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







            _______________________________________________



            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 97, Issue 2



            ***************************************







            Get your own FREE website, FREE domain & FREE mobile app with Company email.

            Know More >

            <2009-11-14-1304-50S.KGD___003_KGD___HH_Z_SAC>_______________________________________________


            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:



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



            _______________________________________________

            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 97, Issue 7

            ***************************************



            Get your own FREE website, FREE domain & FREE mobile app with Company email.
            Know More >
            _______________________________________________
            sac-help mailing list
            sac-help<at>iris.washington.edu
            http://www.iris.washington.edu/mailman/listinfo/sac-help


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


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


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




            --
            ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
            ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
            Satish Maurya
            Ph.D. student
            Dept. of Seismology
            IPGP, France
            ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
            ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
            _______________________________________________
            sac-help mailing list
            sac-help<at>iris.washington.edu
            http://www.iris.washington.edu/mailman/listinfo/sac-help


  • Sheila Peacock
    2013-11-27 16:10:26
    Dear Shekar,

    Guralp seismometers (CMG-40T) all have different response
    files. If you know the serial number of your seismometer,
    email it in the subject line to caldoc<at>guralp.com
    and you will get a data sheet for it by return email.
    Beware, though, that the poles and zeros listed therein
    are in velocity in Hz, not displacement in radians.
    You have to multiply them by 2pi and put in an additional
    zero at (0.0,0.0) to get them into displacement units in
    a form that SAC expects.

    You will then have to calculate the "constant" (or normalising
    factor) using the revised poles-and-zeros, rather than using the
    value on the Guralp sheet.

    Some Guralp seismometers have poles and zeros listed in
    acceleration units - check the units in which the
    "velocity output" is given on the first page of the
    calibration sheet. If so then you have to add TWO
    zeros at (0.0,0.0). I think a GMG-40T calibration is likely to
    be in velocity units, though.

    Sheila Peacock,
    Blacknest.


    PS does anyone on the list know of a good self-contained
    program for calculating the constant (normalising factor)
    of a set of poles and zeros? Or even a good website that
    describes how it's done? Encapsulating the explanation
    into an email without recourse to mathematical symbols
    is something of a challenge.

15:59:47 v.e73c6799