[SAC-HELP] Instrument responce files
Milton Plasencia
mpplasencia at gmail.com
Wed Nov 27 00:47:32 PST 2013
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.iris.washington.edu/pipermail/sac-help/attachments/20131127/94c02f78/attachment.html>
More information about the sac-help
mailing list