response files

Branden Christensen

unread,

Sep 22, 2012, 11:10:45 AM9/22/12

to Earthworm Community Forum

Hi all:

Does anyone know if there is a program to create RESP files from a list of values? PDCC does not seem to have this option.

Or to generate a RESP file from the output of seisan's resp program?

Thanks for any insights.

Kind Regards,

Branden C. Christensen

Director, OSOP

branden.christensen@osop.com.pa

www.osop.com.pa

Volcan, Panama

mwit...@memphis.edu

unread,

Sep 22, 2012, 2:00:50 PM9/22/12

to Earthworm Community Forum

I used vi the last time I had to do it for something that wasn't in the nrl.

SEED is alot like UNIX; its very user friendly but is picky about who its
friends are.

Mitch

Center for Earthquake Research and Information (CERI)
University of Memphis Ph: 901-678-4940
Memphis, TN 38152 Fax: 901-678-4734

On Sat, 22 Sep 2012, Branden Christensen wrote:

> Hi all:
>
>
> Does anyone know if there is a program to create RESP files from a list of
> values? PDCC does not seem to have this option.
>
> Or to generate a RESP file from the output of seisan's resp program?
>
>
> Thanks for any insights.
>
>
> Kind Regards,
>
>

> *Branden C. Christensen*
> Director, OSOP
> *branden.c...@osop.com.pa*
> www.osop.com.pa
> Volcan, Panama
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to earthwo...@googlegroups.com
>
> To unsubscribe from this group, send an email to
> earthworm_for...@googlegroups.com
>
> For more options, visit this group at
> http://groups.google.com/group/earthworm_forum?hl=en
>

Paul Friberg

unread,

Sep 22, 2012, 2:30:04 PM9/22/12

to earthwo...@googlegroups.com

Look into writing your own RESP files (they are just ASCII) and then rolling dataless with ISTI's SHAPE.

Cheers,

Paul

On Sat, Sep 22, 2012 at 2:00 PM, <mwit...@memphis.edu> wrote:

I used vi the last time I had to do it for something that wasn't in the nrl.

SEED is alot like UNIX; its very user friendly but is picky about who its
friends are.

Mitch

Center for Earthquake Research and Information (CERI)
University of Memphis Ph: 901-678-4940
Memphis, TN 38152 Fax: 901-678-4734

On Sat, 22 Sep 2012, Branden Christensen wrote:

Hi all:

Does anyone know if there is a program to create RESP files from a list of
values? PDCC does not seem to have this option.

Or to generate a RESP file from the output of seisan's resp program?

Thanks for any insights.

Kind Regards,

*Branden C. Christensen*
Director, OSOP
*branden.c...@osop.com.pa*
www.osop.com.pa
Volcan, Panama

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--
===================================
Paul Friberg p.fr...@isti.com
ISTI==Instrumental Software Technologies, Inc.
Phone 845.256.9290 FAX 845.256.9299

Branden Christensen

unread,

Sep 22, 2012, 2:36:46 PM9/22/12

to earthwo...@googlegroups.com

Mitch and Paul:

Ha ha. Thanks for taking the piss. I know they are ascii and easy to edit. I was looking for a less mundane solution.

SHAPE? Never heard of it. I will investigate...

Enjoy the weekend.

B

To post to this group, send an email to earthwo...@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_for...@googlegroups.com

Paul Friberg

unread,

Sep 22, 2012, 4:32:11 PM9/22/12

to earthwo...@googlegroups.com

SHAPE: linked via the ORFEUS seed tool list: http://bigboy.isti.com/~sid/shape/

Aleksey Emanov

unread,

Sep 23, 2012, 2:36:30 AM9/23/12

to earthwo...@googlegroups.com

Hi, Branden,

if you have a datalessSEED already you could use one of OBSPY converters - obspy-dataless2resp .

$ obspy-dataless2resp

Usage: obspy-dataless2resp [options] filename

Options:

-h, --help show this help message and exit

-d, --debug show debugging information

-q, --quiet non verbose mode

-z, --zipped Pack files of one station into a ZIP archive.

$

Regards,

Aleksey Emanov

22.09.2012, в 22:10, Branden Christensen написал(а):

Paul Friberg

unread,

Sep 23, 2012, 8:15:48 AM9/23/12

to earthwo...@googlegroups.com

Thanks for that ObsPy tip Aleksey!

Note that if you have IRIS's rdseed or jrdseed installed you can get RESP files out from a dataless too.

Cheers,

Paul

Mitchell M Withers (mwithers)

unread,

Sep 23, 2012, 11:10:29 AM9/23/12

to earthwo...@googlegroups.com

That's one of my SEED pet-peeves. If you have a chicken you can make an egg...

From: earthwo...@googlegroups.com [earthwo...@googlegroups.com] on behalf of Paul Friberg [p.fr...@gmail.com]
Sent: Sunday, September 23, 2012 7:15 AM
To: earthwo...@googlegroups.com
Subject: Re: [Earthworm Forum] response files

Jean-Marie SAUREL

unread,

Sep 24, 2012, 8:00:17 AM9/24/12

to earthwo...@googlegroups.com

Hi Branden,

What kind of values do you have ?
Is this poles and zeros ?
Or maybe, poles and zeros for the whole response, like the PZ SAC files ?

I don't practice Seisan.

Regards.

Jean-Marie.

Le 22/09/2012 16:10, Branden Christensen a �crit :

> Hi all:
>
>
> Does anyone know if there is a program to create RESP files from a list
> of values? PDCC does not seem to have this option.
>
> Or to generate a RESP file from the output of seisan's resp program?
>
>
> Thanks for any insights.
>
>
> Kind Regards,
>

> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.pa <mailto:branden.c...@osop.com.pa>*
> www.osop.com.pa <http://www.osop.com.pa>

> Volcan, Panama
>
>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to earthwo...@googlegroups.com
>
> To unsubscribe from this group, send an email to
> earthworm_for...@googlegroups.com
>
> For more options, visit this group at
> http://groups.google.com/group/earthworm_forum?hl=en

--

--------------------------------------
Institut de Physique du Globe de Paris
Observatoire Volcanologique et Sismologique de Martinique
Morne des Cadets
97250 Fonds St Denis
0596784141
Martinique - FWI

Sid Hellman

unread,

Sep 24, 2012, 8:34:54 AM9/24/12

to earthwo...@googlegroups.com

with SHAPE you can go from RESP to dataless to RESP (via rdseed) to dataless, ad infinitum. You just need to supply the missing information that the dataless needs in the configuration file (HAT file).

Not necessarily for the faint of heart, but it works.

Sid

   Sidney Hellman - s.he...@isti.com - http://www.isti.com
   (518)602-0001 - office    (518)602-0002 - fax
   Instrumental Software Technologies, Inc.
Seismic Software and Services

Branden Christensen

unread,

Sep 24, 2012, 8:39:06 AM9/24/12

to earthwo...@googlegroups.com

Thank you all for your suggests.

Hi Jean-Marie:

Unfortunately I do not have poles and zeros nor direct access to the instrumentation. If I had some poles and zeros I would be home free as I could simply take the brute force method suggested by others and edit one of my existing RESP files. I am starting with the natural period, generator constant, damping ratio, amplifier gain, recording gain, filter information and trying to arrive at a RESP file.

Saludos,

Branden C. Christensen

Director, OSOP

branden.christensen@osop.com.pa

www.osop.com.pa

Volcan, Panama

On Mon, Sep 24, 2012 at 12:00 PM, Jean-Marie SAUREL <sau...@ipgp.fr> wrote:

Hi Branden,

What kind of values do you have ?
Is this poles and zeros ?
Or maybe, poles and zeros for the whole response, like the PZ SAC files ?

I don't practice Seisan.

Regards.

Jean-Marie.

Le 22/09/2012 16:10, Branden Christensen a écrit :

Hi all:

Does anyone know if there is a program to create RESP files from a list
of values? PDCC does not seem to have this option.

Or to generate a RESP file from the output of seisan's resp program?

Thanks for any insights.

Kind Regards,

*Branden C. Christensen*
Director, OSOP

*branden.c...@osop.com.pa <mailto:branden.christensen@osop.com.pa>*
www.osop.com.pa <http://www.osop.com.pa>

Volcan, Panama

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--
--------------------------------------
Institut de Physique du Globe de Paris
Observatoire Volcanologique et Sismologique de Martinique
Morne des Cadets
97250 Fonds St Denis
0596784141
Martinique - FWI

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

mwit...@memphis.edu

unread,

Sep 24, 2012, 8:44:38 AM9/24/12

to earthwo...@googlegroups.com

Oh, so you have the genome and you're trying to make a chicken.

Mitch

Center for Earthquake Research and Information (CERI)
University of Memphis Ph: 901-678-4940
Memphis, TN 38152 Fax: 901-678-4734

On Mon, 24 Sep 2012, Branden Christensen wrote:

> Thank you all for your suggests.
>
> Hi Jean-Marie:
>
>
> Unfortunately I do not have poles and zeros nor direct access to the
> instrumentation. If I had some poles and zeros I would be home free as I
> could simply take the brute force method suggested by others and edit one
> of my existing RESP files. I am starting with the natural period, generator
> constant, damping ratio, amplifier gain, recording gain, filter information
> and trying to arrive at a RESP file.
>
>
> Saludos,
>
>

> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.pa*
> www.osop.com.pa
> Volcan, Panama

>
>
>
>
>
> On Mon, Sep 24, 2012 at 12:00 PM, Jean-Marie SAUREL <sau...@ipgp.fr> wrote:
>
>> Hi Branden,
>>
>> What kind of values do you have ?
>> Is this poles and zeros ?
>> Or maybe, poles and zeros for the whole response, like the PZ SAC files ?
>>
>> I don't practice Seisan.
>>
>> Regards.
>>
>> Jean-Marie.
>>
>> Le 22/09/2012 16:10, Branden Christensen a écrit :
>>
>>> Hi all:
>>>
>>>
>>> Does anyone know if there is a program to create RESP files from a list
>>> of values? PDCC does not seem to have this option.
>>>
>>> Or to generate a RESP file from the output of seisan's resp program?
>>>
>>>
>>> Thanks for any insights.
>>>
>>>
>>> Kind Regards,
>>>
>>> *Branden C. Christensen*
>>> Director, OSOP

>>> *branden.c...@osop.com.**pa <branden.c...@osop.com.pa><mailto:
>>> branden.christensen@**osop.com.pa <branden.c...@osop.com.pa>>*

>>> www.osop.com.pa <http://www.osop.com.pa>
>>>
>>> Volcan, Panama
>>>
>>>
>>>
>>> --
>>> You received this message because you are subscribed to the Google
>>> Groups "Earthworm Community Forum" group.
>>>

>>> To post to this group, send an email to earthworm_forum@googlegroups.**
>>> com <earthwo...@googlegroups.com>

>>>
>>> To unsubscribe from this group, send an email to

>>> earthworm_forum+unsubscribe@**googlegroups.com<earthworm_forum%2Bunsu...@googlegroups.com>

>>>
>>> For more options, visit this group at

>>> http://groups.google.com/**group/earthworm_forum?hl=en<http://groups.google.com/group/earthworm_forum?hl=en>
>>>
>>
>>
>> --
>> ------------------------------**--------

>> Institut de Physique du Globe de Paris
>> Observatoire Volcanologique et Sismologique de Martinique
>> Morne des Cadets
>> 97250 Fonds St Denis
>> 0596784141
>> Martinique - FWI
>>
>>
>> --
>> You received this message because you are subscribed to the Google
>> Groups "Earthworm Community Forum" group.
>>

>> To post to this group, send an email to earthworm_forum@googlegroups.**com<earthwo...@googlegroups.com>

>>
>> To unsubscribe from this group, send an email to

>> earthworm_forum+unsubscribe@**googlegroups.com<earthworm_forum%2Bunsu...@googlegroups.com>

>>
>> For more options, visit this group at

>> http://groups.google.com/**group/earthworm_forum?hl=en<http://groups.google.com/group/earthworm_forum?hl=en>

>>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>

> To post to this group, send an email to earthwo...@googlegroups.com

>
> To unsubscribe from this group, send an email to

> earthworm_for...@googlegroups.com

Jean-Marie SAUREL

unread,

Sep 24, 2012, 9:05:40 AM9/24/12

to earthwo...@googlegroups.com, lemarchand

Hola,

Arnaud Lemarchand, a colleague from Paris, uses a bunch of home-made
matlab script to build a POD directory structure with the corresponding
SEED blockettes and then runs POD to build the dataless.

The input of the whole system is quite similar as what you have :
seismometer type, gain, damping ....

The idea being that with those values you are able to construct
corresponding blockettes, then POD (or other softwares, I think) creates
the dataless.

Sort of synthetic egg generator ...... ;)

Jean-Marie.

Le 24/09/2012 13:39, Branden Christensen a �crit :

> Thank you all for your suggests.
>
> Hi Jean-Marie:
>
>
> Unfortunately I do not have poles and zeros nor direct access to the
> instrumentation. If I had some poles and zeros I would be home free as I
> could simply take the brute force method suggested by others and edit
> one of my existing RESP files. I am starting with the natural period,
> generator constant, damping ratio, amplifier gain, recording gain,
> filter information and trying to arrive at a RESP file.
>
>
> Saludos,
>
>

> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.pa <mailto:branden.c...@osop.com.pa>*

> www.osop.com.pa <http://www.osop.com.pa>
> Volcan, Panama
>
>
>
>
>

> On Mon, Sep 24, 2012 at 12:00 PM, Jean-Marie SAUREL <sau...@ipgp.fr
> <mailto:sau...@ipgp.fr>> wrote:
>
> Hi Branden,
>
> What kind of values do you have ?
> Is this poles and zeros ?
> Or maybe, poles and zeros for the whole response, like the PZ SAC
> files ?
>
> I don't practice Seisan.
>
> Regards.
>
> Jean-Marie.
>

> Le 22/09/2012 16:10, Branden Christensen a �crit :

>
> Hi all:
>
>
> Does anyone know if there is a program to create RESP files from
> a list
> of values? PDCC does not seem to have this option.
>
> Or to generate a RESP file from the output of seisan's resp program?
>
>
> Thanks for any insights.
>
>
> Kind Regards,
>
> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.__pa
> <mailto:branden.c...@osop.com.pa>
> <mailto:branden.christensen@__osop.com.pa
> <mailto:branden.c...@osop.com.pa>>*
> www.osop.com.pa <http://www.osop.com.pa> <http://www.osop.com.pa>

>
> Volcan, Panama
>
>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to

> earthworm_forum@googlegroups.__com
> <mailto:earthwo...@googlegroups.com>

>
> To unsubscribe from this group, send an email to

> earthworm_forum+unsubscribe@__googlegroups.com
> <mailto:earthworm_forum%2Bunsu...@googlegroups.com>

>
> For more options, visit this group at

> http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>
>
>
>
> --
> ------------------------------__--------

> Institut de Physique du Globe de Paris
> Observatoire Volcanologique et Sismologique de Martinique
> Morne des Cadets
> 97250 Fonds St Denis
> 0596784141
> Martinique - FWI
>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to

> earthworm_forum@googlegroups.__com
> <mailto:earthwo...@googlegroups.com>

>
> To unsubscribe from this group, send an email to

> earthworm_forum+unsubscribe@__googlegroups.com
> <mailto:earthworm_forum%2Bunsu...@googlegroups.com>

>
> For more options, visit this group at

> http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>

>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>

> To post to this group, send an email to earthwo...@googlegroups.com

>
> To unsubscribe from this group, send an email to

> earthworm_for...@googlegroups.com

Paul Friberg

unread,

Sep 24, 2012, 9:11:03 AM9/24/12

to earthwo...@googlegroups.com

Jean-Marie,

Nice. Thanks for that info. Has Arnaud linked his tools anywhere on the web? Too bad MatLab is not really open source.....more licensed code you need to run something :(. It would be nice if this were available as an ObsPy script that could crunch gain/period/damping into a simple P&Z file.....perhaps there is and we are not aware of it?

Cheers,

Paul

On Mon, Sep 24, 2012 at 9:05 AM, Jean-Marie SAUREL <sau...@ipgp.fr> wrote:

Hola,

Arnaud Lemarchand, a colleague from Paris, uses a bunch of home-made matlab script to build a POD directory structure with the corresponding SEED blockettes and then runs POD to build the dataless.

The input of the whole system is quite similar as what you have : seismometer type, gain, damping ....

The idea being that with those values you are able to construct corresponding blockettes, then POD (or other softwares, I think) creates the dataless.

Sort of synthetic egg generator ...... ;)

Jean-Marie.

Le 24/09/2012 13:39, Branden Christensen a écrit :

Thank you all for your suggests.

Hi Jean-Marie:

Unfortunately I do not have poles and zeros nor direct access to the
instrumentation. If I had some poles and zeros I would be home free as I
could simply take the brute force method suggested by others and edit
one of my existing RESP files. I am starting with the natural period,
generator constant, damping ratio, amplifier gain, recording gain,
filter information and trying to arrive at a RESP file.

Saludos,

*Branden C. Christensen*
Director, OSOP

*branden.c...@osop.com.pa <mailto:branden.christensen@osop.com.pa>*

www.osop.com.pa <http://www.osop.com.pa>
Volcan, Panama

On Mon, Sep 24, 2012 at 12:00 PM, Jean-Marie SAUREL <sau...@ipgp.fr

<mailto:sau...@ipgp.fr>> wrote:

Hi Branden,

What kind of values do you have ?
Is this poles and zeros ?
Or maybe, poles and zeros for the whole response, like the PZ SAC
files ?

I don't practice Seisan.

Regards.

Jean-Marie.

Le 22/09/2012 16:10, Branden Christensen a écrit :

Hi all:

Does anyone know if there is a program to create RESP files from
a list
of values? PDCC does not seem to have this option.

Or to generate a RESP file from the output of seisan's resp program?

Thanks for any insights.

Kind Regards,

*Branden C. Christensen*
Director, OSOP

*branden.c...@osop.com.__pa
<mailto:branden.christensen@osop.com.pa>
<mailto:branden.christensen@__osop.com.pa
<mailto:branden.christensen@osop.com.pa>>*

www.osop.com.pa <http://www.osop.com.pa> <http://www.osop.com.pa>

Volcan, Panama

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to
earthworm_forum@googlegroups.__com

<mailto:earthworm_forum@googlegroups.com>

To unsubscribe from this group, send an email to
earthworm_forum+unsubscribe@__googlegroups.com
<mailto:earthworm_forum%2Bunsu...@googlegroups.com>

For more options, visit this group at
http://groups.google.com/__group/earthworm_forum?hl=en
<http://groups.google.com/group/earthworm_forum?hl=en>

--
------------------------------__--------

Institut de Physique du Globe de Paris
Observatoire Volcanologique et Sismologique de Martinique
Morne des Cadets
97250 Fonds St Denis
0596784141
Martinique - FWI

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to
earthworm_forum@googlegroups.__com

<mailto:earthworm_forum@googlegroups.com>

To unsubscribe from this group, send an email to
earthworm_forum+unsubscribe@__googlegroups.com
<mailto:earthworm_forum%2Bunsu...@googlegroups.com>

For more options, visit this group at
http://groups.google.com/__group/earthworm_forum?hl=en

<http://groups.google.com/group/earthworm_forum?hl=en>

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--
--------------------------------------
Institut de Physique du Globe de Paris
Observatoire Volcanologique et Sismologique de Martinique
Morne des Cadets
97250 Fonds St Denis
0596784141
Martinique - FWI

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--

Ruben Soares Luís

unread,

Sep 24, 2012, 9:12:30 AM9/24/12

to earthwo...@googlegroups.com

One can use Octave :)

To post to this group, send an email to earthwo...@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_for...@googlegroups.com

Jean-Marie SAUREL

unread,

Sep 24, 2012, 9:24:44 AM9/24/12

to earthwo...@googlegroups.com, lemarchand

Paul,

I'm not sure how Arnaud's scripts are universal.
I think they are very tighed to the different signal path we have in
IPGP observatories (sometimes straigh analog transmission, sometimes
with a bit of digital transmission in the middle).

The fact is that Matlab (and perhaps Octave, or Scilab) provides some
tools to synthetise some common filters based on damping, cut-off and
then can ouput the result in a wide variety of formats (poles and zeros,
off course, but also amplitude/phase).
I'm not sure at all, but that might be the path Arnaud took.

Maybe this should be feasible mixing some analog filter design Python
libraries (is they exists) with ObsPy.

Regards.

Jean-Marie.

Le 24/09/2012 14:11, Paul Friberg a �crit :

> Jean-Marie,
>
> Nice. Thanks for that info. Has Arnaud linked his tools anywhere on the
> web? Too bad MatLab is not really open source.....more licensed code you
> need to run something :(. It would be nice if this were available as an
> ObsPy script that could crunch gain/period/damping into a simple P&Z
> file.....perhaps there is and we are not aware of it?
>
> Cheers,
>
> Paul
>
> On Mon, Sep 24, 2012 at 9:05 AM, Jean-Marie SAUREL <sau...@ipgp.fr
> <mailto:sau...@ipgp.fr>> wrote:
>
> Hola,
>
> Arnaud Lemarchand, a colleague from Paris, uses a bunch of home-made
> matlab script to build a POD directory structure with the
> corresponding SEED blockettes and then runs POD to build the dataless.
>
> The input of the whole system is quite similar as what you have :
> seismometer type, gain, damping ....
>
> The idea being that with those values you are able to construct
> corresponding blockettes, then POD (or other softwares, I think)
> creates the dataless.
>
>
> Sort of synthetic egg generator ...... ;)
>
> Jean-Marie.
>

> Le 24/09/2012 13:39, Branden Christensen a �crit :

>
> Thank you all for your suggests.
>
> Hi Jean-Marie:
>
>
> Unfortunately I do not have poles and zeros nor direct access to the
> instrumentation. If I had some poles and zeros I would be home
> free as I
> could simply take the brute force method suggested by others and
> edit
> one of my existing RESP files. I am starting with the natural
> period,
> generator constant, damping ratio, amplifier gain, recording gain,
> filter information and trying to arrive at a RESP file.
>
>
> Saludos,
>
>
> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.__pa
> <mailto:branden.c...@osop.com.pa>
> <mailto:branden.christensen@__osop.com.pa

> <mailto:branden.c...@osop.com.pa>>*

> www.osop.com.pa <http://www.osop.com.pa> <http://www.osop.com.pa>
> Volcan, Panama
>
>
>
>
>

> On Mon, Sep 24, 2012 at 12:00 PM, Jean-Marie SAUREL
> <sau...@ipgp.fr <mailto:sau...@ipgp.fr>

> <mailto:sau...@ipgp.fr <mailto:sau...@ipgp.fr>>> wrote:
>
> Hi Branden,
>
> What kind of values do you have ?
> Is this poles and zeros ?
> Or maybe, poles and zeros for the whole response, like the
> PZ SAC
> files ?
>
> I don't practice Seisan.
>
> Regards.
>
> Jean-Marie.
>

> Le 22/09/2012 16:10, Branden Christensen a �crit :

>
> Hi all:
>
>
> Does anyone know if there is a program to create RESP
> files from
> a list
> of values? PDCC does not seem to have this option.
>
> Or to generate a RESP file from the output of seisan's
> resp program?
>
>
> Thanks for any insights.
>
>
> Kind Regards,
>
> *Branden C. Christensen*
> Director, OSOP
> *branden.c...@osop.com

> <mailto:branden.c...@osop.com>.____pa
> <mailto:branden.christensen@__osop.com.pa
> <mailto:branden.c...@osop.com.pa>>
> <mailto:branden.christensen@
> <mailto:branden.christensen@>____osop.com.pa <http://osop.com.pa>
> <mailto:branden.christensen@__osop.com.pa
> <mailto:branden.c...@osop.com.pa>>>*

> www.osop.com.pa <http://www.osop.com.pa>
> <http://www.osop.com.pa> <http://www.osop.com.pa>
>
>
> Volcan, Panama
>
>
>
> --
> You received this message because you are subscribed to
> the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to

> earthworm_forum@googlegroups.____com
> <mailto:earthworm_forum@__googlegroups.com
> <mailto:earthwo...@googlegroups.com>>

>
>
> To unsubscribe from this group, send an email to

> earthworm_forum+unsubscribe@____googlegroups.com
> <http://googlegroups.com>
> <mailto:earthworm_forum%__2Buns...@googlegroups.com
> <mailto:earthworm_forum%252Buns...@googlegroups.com>__>

>
>
> For more options, visit this group at

> http://groups.google.com/____group/earthworm_forum?hl=en
> <http://groups.google.com/__group/earthworm_forum?hl=en>

> <http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>>
>
>
>
> --

> ------------------------------____--------

>
> Institut de Physique du Globe de Paris
> Observatoire Volcanologique et Sismologique de Martinique
> Morne des Cadets
> 97250 Fonds St Denis
> 0596784141
> Martinique - FWI
>
>
> --
> You received this message because you are subscribed to the
> Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to

> earthworm_forum@googlegroups.____com
> <mailto:earthworm_forum@__googlegroups.com
> <mailto:earthwo...@googlegroups.com>>

>
>
> To unsubscribe from this group, send an email to

> earthworm_forum+unsubscribe@____googlegroups.com
> <http://googlegroups.com>
> <mailto:earthworm_forum%__2Buns...@googlegroups.com
> <mailto:earthworm_forum%252Buns...@googlegroups.com>__>

>
>
> For more options, visit this group at

> http://groups.google.com/____group/earthworm_forum?hl=en
> <http://groups.google.com/__group/earthworm_forum?hl=en>

>
> <http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>>
>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to

> earthworm_forum@googlegroups.__com
> <mailto:earthwo...@googlegroups.com>

>
> To unsubscribe from this group, send an email to
> earthworm_forum+unsubscribe@__googlegroups.com
> <mailto:earthworm_forum%2Bunsu...@googlegroups.com>
>
> For more options, visit this group at
> http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>
>
>
>
> --
> ------------------------------__--------
> Institut de Physique du Globe de Paris
> Observatoire Volcanologique et Sismologique de Martinique
> Morne des Cadets
> 97250 Fonds St Denis
> 0596784141
> Martinique - FWI
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to
> earthworm_forum@googlegroups.__com

> <mailto:earthwo...@googlegroups.com>

>
> To unsubscribe from this group, send an email to
> earthworm_forum+unsubscribe@__googlegroups.com
> <mailto:earthworm_forum%2Bunsu...@googlegroups.com>
>
> For more options, visit this group at
> http://groups.google.com/__group/earthworm_forum?hl=en
> <http://groups.google.com/group/earthworm_forum?hl=en>
>
>
>
>
> --

> ===================================
> Paul Friberg p.fr...@isti.com <mailto:p.fr...@isti.com>

> ISTI==Instrumental Software Technologies, Inc.
> Phone 845.256.9290 FAX 845.256.9299
>

> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>

> To post to this group, send an email to earthwo...@googlegroups.com

>
> To unsubscribe from this group, send an email to

> earthworm_for...@googlegroups.com

Branden Christensen

unread,

Sep 24, 2012, 1:27:43 PM9/24/12

to earthwo...@googlegroups.com

Hi all:

I have made significant progress in trying to go from the genome to the chicken :)

Here is how I am doing it:

I feed SEISAN's resp program the following values (option GSE FAP, SEISMOMETER):

seismometer period: 0.222

generator constant (V/m/s): 27.5

Damping ratio: 0.647

Amplifier Gain (dB): 0.0 (or 1)

Recording Gain (counts/Volt): 8388.608

Filter Constants, F = 12.5, Poles = 1

Seisan spits out a GSE file with the following poles and zeros (displacement):

-0.18311806E+02 0.21580477E+02

-0.18311806E+02 -0.21580477E+02

-0.78539818E+02 -0.68661689E-05

0.00000000E+00 0.00000000E+00

I then manually insert this information into a RESP file:

# << IRIS SEED Reader, Release 5.2 >>

#

# ======== CHANNEL RESPONSE DATA ========

B050F03 Station: SD04

B050F16 Network: CX

B052F03 Location: ??

B052F04 Channel: SHZ

B052F22 Start date: 2011,001,00:00:00.0000

B052F23 End date: No Ending Time

# =======================================

# + +--------------------------------------------+ +

# + | Response (Poles & Zeros), SD04 ch SHZ | +

# + +--------------------------------------------+ +

#

B053F03 Transfer function type: A [Laplace Transform (Rad/sec)]

B053F04 Stage sequence number: 1

B053F05 Response in units lookup: M/S - Velocity in Meters Per Second

B053F06 Response out units lookup: V - Volts

B053F07 A0 normalization factor: 1.1E-03

B053F08 Normalization frequency: 1

B053F09 Number of zeroes: 2

B053F14 Number of poles: 3

# Complex zeroes:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 2 -0.785398E+02 -0.686617E-05 0.000000E+00 0.000000E+00

#

# + +---------------------------------------+ +

# + | Channel Gain, SD04 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 1

B058F04 Gain: 2.75E+01

B058F05 Frequency of gain: 1.000000E+00 HZ

B058F06 Number of calibrations: 0

#

# + +-------------------------------------------+ +

# + | Response (Coefficients), SD04 ch SHZ | +

# + +-------------------------------------------+ +

#

B054F03 Transfer function type: D

B054F04 Stage sequence number: 2

B054F05 Response in units lookup: V - Volts

B054F06 Response out units lookup: COUNTS - Digital Counts

B054F07 Number of numerators: 0

B054F10 Number of denominators: 0

# Numerator coefficients:

# i, coefficient, error

#

# + +------------------------------+ +

# + | Decimation, SD04 ch SHZ | +

# + +------------------------------+ +

#

B057F03 Stage sequence number: 2

B057F04 Input sample rate: 5.000000E+01

B057F05 Decimation factor: 1

B057F06 Decimation offset: 0

B057F07 Estimated delay (seconds): 0.000000E+00

B057F08 Correction applied (seconds): 0.000000E+00

#

# + +---------------------------------------+ +

# + | Channel Gain, SD04 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 2

B058F04 Gain: 8.388608E+03

B058F05 Frequency of gain: 0.000000E+00 HZ

B058F06 Number of calibrations: 0

#

# + +---------------------------------------+ +

# + | Channel Sensitivity, SD04 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 0

B058F04 Sensitivity: 2.30687E+05

B058F05 Frequency of sensitivity: 1.000000E+00 HZ

B058F06 Number of calibrations: 0

#

Where:
- One of the zeros has been deleted because RESP files required Velocity

- "Sensitivity" is the Generator Constant (V/m/s) * Recording Gain (c/V) so counts/meters/second

Now, the "A0 normalization factor" is also known as the "Poles and zeros constant". I do not know how to calculate it. It controls the location of the response curve in amplitude space. I have used PQLX to play with it and have determined (see attached sequence), based on my experience with PSDs, that the value is ~1.1E+04. That is to say in the +04 magnitude range looks reasonable).

Any ideas how to calculate it directly?

I hope this comes in useful for someone in the future!

- Branden

ao_norm_factor_1.1E+02.png

ao_norm_factor_1.1E+03.png

ao_norm_factor_1.1E+04.png

ao_norm_factor_1.1E+05.png

ao_norm_factor_1.1E+06.png

Jean-Marie SAUREL

unread,

Sep 24, 2012, 1:58:27 PM9/24/12

to earthwo...@googlegroups.com

Branden,

The A0 normalization factor is the factor by which you multiply the
poles&zeros based curve to achieve a gain of 1 at the desired frequency
(typically 1Hz).
So the idea is to plot the pole and zero stage response curve with A0=1
(with JplotResp, for example), then pick the amplitude at the
normalization frequency, and reverse it (1/picked amplitude).
Normaly, I would expect it would do the job.

Alternatively, I think there is a gse2seed utility somewhere on IRIS or
Orfeus website.
You should be able to use it to have directly a SEED dataless. The only
trick is to have some waveform within the GSE, but I don't think this
would be an issue.
A colleague from El Salvador has been using this method, a year or two ago.

Jean-Marie.

Le 24/09/2012 18:27, Branden Christensen a �crit :

> Hi all:
>
>
> I have made significant progress in trying to go from the genome to the
> chicken :)
>
> Here is how I am doing it:
>
>
> I feed SEISAN's resp program the following values (option GSE FAP,
> SEISMOMETER):
>
> seismometer period: 0.222
> generator constant (V/m/s): 27.5
> Damping ratio: 0.647
> Amplifier Gain (dB): 0.0 (or 1)
> Recording Gain (counts/Volt): 8388.608
> Filter Constants, F = 12.5, Poles = 1
>
> Seisan spits out a GSE file with the following poles and zeros
> (displacement):
>
> -0.18311806E+02 0.21580477E+02
> -0.18311806E+02 -0.21580477E+02
> -0.78539818E+02 -0.68661689E-05
> 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.00000000E+00
> 0.00000000E+00 0.00000000E+00
>
>
> I then manually insert this information into a RESP file:
>
> #<< IRIS SEED Reader, Release 5.2 >>
> #

> #======== CHANNEL RESPONSE DATA ========

> B050F03 Station: SD04
> B050F16 Network: CX
> B052F03 Location: ??
> B052F04 Channel: SHZ
> B052F22 Start date: 2011,001,00:00:00.0000
> B052F23 End date: No Ending Time
> #=======================================
> #+ +--------------------------------------------+
> +

> #+ | Response (Poles & Zeros), SD04 ch SHZ |
> +
> #+ +--------------------------------------------+

> +
> #
> B053F03 Transfer function type: A [Laplace Transform
> (Rad/sec)]
> B053F04 Stage sequence number: 1
> B053F05 Response in units lookup: M/S - Velocity in
> Meters Per Second
> B053F06 Response out units lookup: V - Volts
> B053F07 A0 normalization factor: 1.1E-03
> B053F08 Normalization frequency: 1
> B053F09 Number of zeroes: 2
> B053F14 Number of poles: 3

> #Complex zeroes:

> # i real imag real_error imag_error
> B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
> B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

> #Complex poles:

> # i real imag real_error imag_error
> B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00
> B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00
> B053F15-18 2 -0.785398E+02 -0.686617E-05 0.000000E+00 0.000000E+00
> #
> #+ +---------------------------------------+
> +

> #+ | Channel Gain, SD04 ch SHZ |
> +
> #+ +---------------------------------------+

> +
> #
> B058F03 Stage sequence number: 1
> B058F04 Gain: 2.75E+01
> B058F05 Frequency of gain: 1.000000E+00 HZ
> B058F06 Number of calibrations: 0
> #
> #+ +-------------------------------------------+
> +

> #+ | Response (Coefficients), SD04 ch SHZ |
> +
> #+ +-------------------------------------------+

> +
> #
> B054F03 Transfer function type: D
> B054F04 Stage sequence number: 2
> B054F05 Response in units lookup: V - Volts
> B054F06 Response out units lookup: COUNTS - Digital Counts
> B054F07 Number of numerators: 0
> B054F10 Number of denominators: 0

> #Numerator coefficients:

> # i, coefficient, error
> #
> #+ +------------------------------+
> +

> #+ | Decimation, SD04 ch SHZ |
> +
> #+ +------------------------------+

> +
> #
> B057F03 Stage sequence number: 2
> B057F04 Input sample rate: 5.000000E+01
> B057F05 Decimation factor: 1
> B057F06 Decimation offset: 0
> B057F07 Estimated delay (seconds): 0.000000E+00
> B057F08 Correction applied (seconds): 0.000000E+00
> #
> #+ +---------------------------------------+
> +

> #+ | Channel Gain, SD04 ch SHZ |
> +
> #+ +---------------------------------------+

> +
> #
> B058F03 Stage sequence number: 2
> B058F04 Gain: 8.388608E+03
> B058F05 Frequency of gain: 0.000000E+00 HZ
> B058F06 Number of calibrations: 0
> #
> #+ +---------------------------------------+
> +

> #+ | Channel Sensitivity, SD04 ch SHZ |
> +
> #+ +---------------------------------------+

> +
> #
> B058F03 Stage sequence number: 0
> B058F04 Sensitivity: 2.30687E+05
> B058F05 Frequency of sensitivity: 1.000000E+00 HZ
> B058F06 Number of calibrations: 0
> #
>
> Where:
> - One of the zeros has been deleted because RESP files required Velocity
> - "Sensitivity" is the Generator Constant (V/m/s) * Recording Gain (c/V)
> so counts/meters/second
>
> Now, the "A0 normalization factor" is also known as the "Poles and zeros
> constant". I do not know how to calculate it. It controls the location
> of the response curve in amplitude space. I have used PQLX to play with
> it and have determined (see attached sequence), based on my experience
> with PSDs, that the value is ~1.1E+04. That is to say in the +04
> magnitude range looks reasonable).
>
> Any ideas how to calculate it directly?
>
>
> I hope this comes in useful for someone in the future!
>
>
> - Branden
>

Branden Christensen

unread,

Sep 24, 2012, 2:55:41 PM9/24/12

to earthwo...@googlegroups.com

Jean-Marie:

The set your AO = 1 in the RESP file and plot with JPlotResp option worked great! It yields an amplitude of 800 at 1 Hz. So the AO factor = 1/800 or 1.25E-03. As you might suspect from the PQLX plots I shared in my last email, this places the PSDs WAY above the high noise model, suggesting that my generator constant, damping ratio, amplifier gain, recording gain and/ or filters are wrong. So time to investigate why my instrumentation specs are spurious.

I tried the gse2seed option for a few hours and couldn't get it to output anything other than a blank file.

Saludos,

B

On Mon, Sep 24, 2012 at 5:58 PM, Jean-Marie SAUREL <sau...@ipgp.fr> wrote:

Branden,

The A0 normalization factor is the factor by which you multiply the poles&zeros based curve to achieve a gain of 1 at the desired frequency (typically 1Hz).
So the idea is to plot the pole and zero stage response curve with A0=1 (with JplotResp, for example), then pick the amplitude at the normalization frequency, and reverse it (1/picked amplitude).
Normaly, I would expect it would do the job.

Alternatively, I think there is a gse2seed utility somewhere on IRIS or Orfeus website.
You should be able to use it to have directly a SEED dataless. The only trick is to have some waveform within the GSE, but I don't think this would be an issue.
A colleague from El Salvador has been using this method, a year or two ago.

Jean-Marie.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

--
--------------------------------------
Institut de Physique du Globe de Paris
Observatoire Volcanologique et Sismologique de Martinique
Morne des Cadets
97250 Fonds St Denis
0596784141
Martinique - FWI

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthworm_forum@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_forum+unsubscribe@googlegroups.com

Philip Crotwell

unread,

Sep 24, 2012, 4:15:18 PM9/24/12

to earthwo...@googlegroups.com

The overall gain of the seismometer (at the given frequency) should be the product of the A0 and the stage 1 gain value. So, if you fix A0 to give you 1 at 1 Hz, then you need to set the stage 1 gain to compensate.

Gotta love how easy and transparent all this is now. I remember back in the old days when we had to type poles and zeros in by hand... oh wait... never mind.

Philip

To post to this group, send an email to earthwo...@googlegroups.com

To unsubscribe from this group, send an email to

earthworm_for...@googlegroups.com

Sid Hellman

unread,

Sep 24, 2012, 4:33:02 PM9/24/12

to earthwo...@googlegroups.com

HI Folks,

We (ISTI) just posted SHAPE directly to our web site. http://www.isti.com/products/shape, or just go to the Products tab to get it.

Enjoy,

Sid

Aleksey Emanov

unread,

Sep 24, 2012, 7:19:12 PM9/24/12

to earthwo...@googlegroups.com

Hi, Branden,

The table part of response in your RESP file is not flat. You can see it with JplotResp plotting it as a velocity.

I guess you should have two zeros and two poles in the final RESP file.

Regards,

Aleksey Emanov

25.09.2012, в 0:27, Branden Christensen написал(а):

--
You received this message because you are subscribed to the Google
Groups "Earthworm Community Forum" group.

To post to this group, send an email to earthwo...@googlegroups.com

To unsubscribe from this group, send an email to
earthworm_for...@googlegroups.com

For more options, visit this group at
http://groups.google.com/group/earthworm_forum?hl=en

<ao_norm_factor_1.1E+02.png><ao_norm_factor_1.1E+03.png><ao_norm_factor_1.1E+04.png><ao_norm_factor_1.1E+05.png><ao_norm_factor_1.1E+06.png>

Branden Christensen

unread,

Sep 24, 2012, 7:26:37 PM9/24/12

to earthwo...@googlegroups.com

Aleksey:

Thanks for looking into this in such detail!

I also noticed this. I think something is amiss in the instrument values I have. So perhaps the poles and number of poles is wrong.

Saludos,

Branden C. Christensen

Director, OSOP

branden.christensen@osop.com.pa

www.osop.com.pa

Volcan, Panama

Aleksey Emanov

unread,

Sep 24, 2012, 7:38:37 PM9/24/12

to earthwo...@googlegroups.com

Hi, Branden,

We calculate poles and zeros in Excel and I've Just using your period and dumping. So I got

two zeros in (0.0) and two poles in

-18,31180583	21,58047732
-18,31180583	-21,58047732

Also poles should be complex conjugate numbers so it is well.

Regards,

Aleksey Emanov

P.S. If you wish I could send the example excel file to you, but sorry, all of the comments and remarks are in russian, so

it not so useful.

25.09.2012, в 6:26, Branden Christensen написал(а):

Jean-Marie SAUREL

unread,

Sep 24, 2012, 8:03:32 PM9/24/12

to earthwo...@googlegroups.com

Branden,

I just figured that your Seisan outputs values in displacement, while you
are building a RESP file in velocity....
That may explain the slope Aleksey noted, and probably also the strange
curves you saw in PQLX.

Check also the units, as for example SAC uses nanometers and no meters in
all their response (velocity, displacement, acceleration). Maybe Seisan
also uses something else than meter (and thus the gain problem).

Buenas noches.

Jean-Marie.

> Hi, Branden,
>
> We calculate poles and zeros in Excel and I've Just using your period and
> dumping. So I got
> two zeros in (0.0) and two poles in
>
> -18,31180583 21,58047732
> -18,31180583 -21,58047732
>
> Also poles should be complex conjugate numbers so it is well.
>
>
> Regards,
> Aleksey Emanov
>
> P.S. If you wish I could send the example excel file to you, but sorry,
> all of the comments and remarks are in russian, so
> it not so useful.
>
> 25.09.2012, в 6:26, Branden Christensen написал(а):
>
>> Aleksey:
>>
>>
>> Thanks for looking into this in such detail!
>>
>> I also noticed this. I think something is amiss in the instrument values
>> I have. So perhaps the poles and number of poles is wrong.
>>
>>
>> Saludos,
>>
>> Branden C. Christensen
>> Director, OSOP

>> branden.c...@osop.com.pa
>> www.osop.com.pa
>> Volcan, Panama

Aleksey Emanov

unread,

Sep 24, 2012, 8:14:30 PM9/24/12

to earthwo...@googlegroups.com

Hi,

Also for such poles you should have A0 close to 1. And choose some different than 1 Hz normalization

frequency since you have 4.5 Hz eigen frequency of the sensor. Normalization frequency should be on the

flat passband, but not on the sloped transition band. So 8 or even 10 Hz should be fine.

Regards,

Aleksey Emanov

25.09.2012, в 6:26, Branden Christensen написал(а):

Branden Christensen

unread,

Sep 24, 2012, 10:24:50 PM9/24/12

to earthwo...@googlegroups.com

All:

Great insights guys.

This whole exercise has really pushed my understanding of poles and zeros to the edge. I see from your comments and suggestions that I am not the only one challenged here. Time to read Scherbaum AGAIN! I plan to make a comprehensive wiki page on methods/ software available to get an egg from a chicken and a chicken from a genome as Mitch put it.

Jean-Marie, please note that the only things I use from the Seisan output are poles and zeros which are not affected by displacement v. velocity, an extra zero aside in the case of displacement. (Or this is how I understand it). Though, as Aleksey pointed out and I discuss below, I had an extra spurious pole in the mix.

And it is worth noting that poles only come in pairs where there is an imaginary component. Theory requires complex conjugate pairs. Here is an example for an STS-2/ Q330 combination:

# Complex zeroes:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 2 -1.515000E+01 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 3 -1.766000E+02 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 4 -4.631000E+02 -4.305000E+02 0.000000E+00 0.000000E+00

B053F10-13 5 -4.631000E+02 4.305000E+02 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -3.700000E-02 -3.700000E-02 0.000000E+00 0.000000E+00

B053F15-18 1 -3.700000E-02 3.700000E-02 0.000000E+00 0.000000E+00

B053F15-18 2 -1.564000E+01 0.000000E+00 0.000000E+00 0.000000E+00

B053F15-18 3 -9.734000E+01 -4.007000E+02 0.000000E+00 0.000000E+00

B053F15-18 4 -9.734000E+01 4.007000E+02 0.000000E+00 0.000000E+00

B053F15-18 5 -3.748000E+02 0.000000E+00 0.000000E+00 0.000000E+00

B053F15-18 6 -5.203000E+02 0.000000E+00 0.000000E+00 0.000000E+00

B053F15-18 7 -1.053000E+04 -1.005000E+04 0.000000E+00 0.000000E+00

B053F15-18 8 -1.053000E+04 1.005000E+04 0.000000E+00 0.000000E+00

B053F15-18 9 -1.330000E+04 0.000000E+00 0.000000E+00 0.000000E+00

B053F15-18 10 -2.550970E+02 0.000000E+00 0.000000E+00 0.000000E+00

Notice, for instance, that there is no pair specified for the pole:

B053F15-18 9 -1.330000E+04 0.000000E+00 0.000000E+00 0.000000E+00

No imaginary component, no need to invoke the rule for complex conjugate pairs.

This does explain why my curve was not flat in velocity as Aleksey so astutely pointed out (though not why it was shifted in dB as seen in PQLX...this is a gain issue I am still hunting down).

I was using this:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 2 -0.785398E+02 -0.686617E-05 0.000000E+00 0.000000E+00

When the only valid possibilities are:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

and

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 2 -0.785398E+02 -0.686617E-05 0.000000E+00 0.000000E+00

B053F15-18 2 -0.785398E+02 0.686617E-05 0.000000E+00 0.000000E+00

I chose the former route of eliminating what looks like a spurious pole and the velocity curve looks great now. This means the SEISAN FAP and GSE FAP outputs from resp were wrong. Programs are still only as smart as their users :D

Unfortunately, my work is not done. To get the flat part of the curve to have a gain of 1 I still need an A0 factor to the -06. So one of my gains must be off by a light year.

B

Aleksey Emanov

unread,

Sep 24, 2012, 11:34:32 PM9/24/12

to earthwo...@googlegroups.com

Dear Branden,

Not exactly understand which the possibility you use of two mentioned.
The first possibility is Ok, but the second is is even more sloped than the case with 3 poles.
Every zero increase the angle of response curve counter clockwise from horizontal at the point of zero coordinates.
And than every pole decrease clockwise at the point of pole coordinates. So to get flat portion of the response
you have to use the same number of poles like zeros number.

Also do not forget to change number of poles in the header.

B053F09 Number of zeroes: 2

B053F14 Number of poles: 2

I've checked first case with JplotResp and with A0=1.0 I got Amp=1.00265 an 8.111 Hz.

Regards,
Aleksey Emanov

25.09.2012, в 9:24, Branden Christensen написал(а):

Branden Christensen

unread,

Sep 24, 2012, 11:52:01 PM9/24/12

to earthwo...@googlegroups.com

Aleksey:

The first.

That is funny, when I set A0 = 1 and use JPlotResp I get an amplitude of 2.313E+06 at 8.11 Hz. Here is the updated resp file:

# << IRIS SEED Reader, Release 5.2 >>

#

# ======== CHANNEL RESPONSE DATA ========

B050F03 Station: SD09

B050F16 Network: CX

B052F03 Location: ??

B052F04 Channel: SHZ

B052F22 Start date: 2011,001,00:00:00.0000

B052F23 End date: No Ending Time

# =======================================

# + +--------------------------------------------+ +

# + | Response (Poles & Zeros), SD09 ch SHZ | +

# + +--------------------------------------------+ +

#

B053F03 Transfer function type: A [Laplace Transform (Rad/sec)]

B053F04 Stage sequence number: 1

B053F05 Response in units lookup: M/S - Velocity in Meters Per Second

B053F06 Response out units lookup: V - Volts

B053F07 A0 normalization factor: 1.0E+01

B053F08 Normalization frequency: 1.0E+01

B053F09 Number of zeroes: 2

B053F14 Number of poles: 2

# Complex zeroes:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

#

# + +---------------------------------------+ +

# + | Channel Gain, SD09 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 1

B058F04 Gain: 2.75E+01

B058F05 Frequency of gain: 1.0E+01 HZ

B058F06 Number of calibrations: 0

#

# + +-------------------------------------------+ +

# + | Response (Coefficients), SD09 ch SHZ | +

# + +-------------------------------------------+ +

#

B054F03 Transfer function type: D

B054F04 Stage sequence number: 2

B054F05 Response in units lookup: V - Volts

B054F06 Response out units lookup: COUNTS - Digital Counts

B054F07 Number of numerators: 0

B054F10 Number of denominators: 0

# Numerator coefficients:

# i, coefficient, error

#

# + +------------------------------+ +

# + | Decimation, SD09 ch SHZ | +

# + +------------------------------+ +

#

B057F03 Stage sequence number: 2

B057F04 Input sample rate: 5.0E+01

B057F05 Decimation factor: 1

B057F06 Decimation offset: 0

B057F07 Estimated delay (seconds): 0.0E+00

B057F08 Correction applied (seconds): 0.0E+00

#

# + +---------------------------------------+ +

# + | Channel Gain, SD09 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 2

B058F04 Gain: 8.388608E+03

B058F05 Frequency of gain: 0.0E+00 HZ

B058F06 Number of calibrations: 0

#

# + +---------------------------------------+ +

# + | Channel Sensitivity, SD09 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 0

B058F04 Sensitivity: 2.30687E+05

B058F05 Frequency of sensitivity: 1.0E+01 HZ

Aleksey Emanov

unread,

Sep 24, 2012, 11:59:05 PM9/24/12

to earthwo...@googlegroups.com

Dear Branden,

1. Just plot only 1st stage in JplotResp, skipping others.

2. You have to void (set to 1.0) sensor gain for the fist stage to get a value of A0

Regards,

Aleksey

25.09.2012, в 10:52, Branden Christensen написал(а):

Aleksey Emanov

unread,

Sep 25, 2012, 12:00:18 AM9/25/12

to earthwo...@googlegroups.com

Dear Branden,

I've checked first case with JplotResp and with A0=1.0 I got Amp=1.00265 an 8.111 Hz.

Sorry that is correct for norm frequency 1.0.

The correct representation for your sensor with sensivity 1.0 @10 Hz should be like that:

# ======== CHANNEL RESPONSE DATA ========

B050F03 Station: SD04

B050F16 Network: CX

B052F03 Location: ??

B052F04 Channel: SHZ

B052F22 Start date: 2011,001,00:00:00.0000

B052F23 End date: No Ending Time

# =======================================

# + +--------------------------------------------+ +

# + | Response (Poles & Zeros), SD04 ch SHZ | +

# + +--------------------------------------------+ +

#

B053F03 Transfer function type: A [Laplace Transform (Rad/sec)]

B053F04 Stage sequence number: 1

B053F05 Response in units lookup: M/S - Velocity in Meters Per Second

B053F06 Response out units lookup: V - Volts

B053F07 A0 normalization factor: 2.040752E+01

B053F08 Normalization frequency: 10.0

B053F09 Number of zeroes: 2

B053F14 Number of poles: 2

# Complex zeroes:

# i real imag real_error imag_error

B053F10-13 0 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

B053F10-13 1 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

# Complex poles:

# i real imag real_error imag_error

B053F15-18 0 -0.183118E+02 0.215805E+02 0.000000E+00 0.000000E+00

B053F15-18 1 -0.183118E+02 -0.215805E+02 0.000000E+00 0.000000E+00

#

# + +---------------------------------------+ +

# + | Channel Gain, SD04 ch SHZ | +

# + +---------------------------------------+ +

#

B058F03 Stage sequence number: 1

B058F04 Gain: 1.

B058F05 Frequency of gain: 1.000000E+01 HZ

B058F06 Number of calibrations: 0

Regards,

Aleksey Emanov

25.09.2012, в 10:52, Branden Christensen написал(а):

Branden Christensen

unread,

Sep 25, 2012, 11:11:56 AM9/25/12

to earthwo...@googlegroups.com

Aleksey, Jean-Marie, Mitch, Paul, Philip, Ruben, Sid et al.:

Thank you for your insights. This has been a lively discussion!

In Summary, To go from a genome to a Chicken:

1. Use Seisan's resp program with options GSE FAP and SEISMOMETER to generate poles and zeros from your instrument parameters (seismometer period, sample rate (sample per second), generator constant (V/m/s), damping ratio, amplifier gain (dB), recording gain (c/V), filter constants). Be careful to only use real poles and complex conjugate pairs!

2. Edit an existing RESP file. If you do not have one, use either 1) PDCC to generate a generic unity RESP file or 2) rdseed -R -f <dataless> to generate a RESP file from an existing dataless one.

Most of the fields in a RESP file are self explanatory. An example with explanations can be found at: http://www.iris.edu/KB/questions/126/What+are+the+fields+in+a+RESP+file%3F.

Not so obvious fields include:

a. The stage 1 gain is = to the generator constant (V/m/s)

b. The channel gain is = to the recording gain (c/V)

c. The sensitivity is = to the stage 1 gain * the channel gain

3. Now calculate the A0 normalization factor.

i. Create a copy of the RESP file.

ii. Set A0 to 1

iii. At stage one, set gain to 1

iv. Delete everything after stage 1, thus simplifying the RESP file and removing most of the contents

v. Plot the new edited RESP file in JPlotResp (choose velocity output)

vi. Read the amplitude at the normalization frequency (or any frequency where the response is flat)

vii. 1/amplitude = A0

viii. In the original RESP file, set A0 to 1/amplitude and set the normalization frequency to the frequency you calculated this at.

To confirm the response file is correct, you can use PQLX. Your PSDs should fall within the bounds of the high and low noise model.

Kind Regards,

Jean-Marie SAUREL

unread,

Sep 25, 2012, 1:51:40 PM9/25/12

to earthwo...@googlegroups.com

Hello Branden,

>
> To confirm the response file is correct, you can use PQLX. Your PSDs
> should fall within the bounds of the high and low noise model.
>

In addition, another easy verification is to remove the response file
from a big teleseismic recording (with SAC, for example "transfer from
evalresp to velocity freq 0.001 0.002 20 50") and compare it with the
signal recorded by a near well known broadband station (again with
response removed).

You will see how outstanding the similarities are between the two, even
if you compare an analog transmitted 2Hz geophone with a STS2 around
between 10s and 100s of signal period.
If you have similar waveforms and an amplitude within 30%, I would
consider the results as good and validated.

Regards.

Jean-Marie SAUREL.

>
> Kind Regards,
>
> *Branden C. Christensen*
> Director, OSOP

> *branden.c...@osop.com.pa <mailto:branden.c...@osop.com.pa>*
> www.osop.com.pa <http://www.osop.com.pa>
> Volcan, Panama

>
>
> --
> You received this message because you are subscribed to the Google
> Groups "Earthworm Community Forum" group.
>
> To post to this group, send an email to earthwo...@googlegroups.com
>
> To unsubscribe from this group, send an email to
> earthworm_for...@googlegroups.com
>
> For more options, visit this group at
> http://groups.google.com/group/earthworm_forum?hl=en

Reply all

Reply to author

Forward