Thread: evalresp transfer failure

Started: 2021-05-20 17:41:55
Last activity: 2021-05-21 08:07:41
Topics: SAC Help
Sheila Peacock
2021-05-20 17:41:55
Dear all,

I extracted data from a SEED file with rdseed and told it to write the RESP files so I could deconvolve the instrument response.

Many of the channels in the SEED file aren't seismic and their RESP files are not useful, but when I tested the response of one of the channels with an impulse, I got an error message as if the "transfer" subprogram had been reading from the wrong file:

SAC> funcgen impulse npts 16384 delta 0.004
SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
Station (-12345 ), Channel (-12345 )
Extracting evresp response for *, *...
Thu May 20 17:26:26 2021 [WARN] check_units; units found ('% - PERCENTAGE') are not supported
Thu May 20 17:26:26 2021 [ERROR] DG.0385F.01.LZ0: undefined units
No transfer function applied!
ERROR 2118: No transfer function applied.

My test directory contains two RESP files:

RESP.DG.0385F.01.HHZ.good RESP.DG.0385F.01.LZ0.bad

If I delete the "bad" file or if I rename it e.g. to badresp.DG.0385F.01.LZO.bad then the problem goes away:

SAC> funcgen impulse npts 16384 delta 0.004
SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
Station (-12345 ), Channel (-12345 )
Extracting evresp response for *, *...

What worries me is that if SAC is not using the RESP file that I told it to use, then if I had several "good" ones in my directory, how can I be sure it used the one I specified in the "transfer "command?

I've attached the "good" and "bad" files. I'm not interested in why the bad file is bad, just why SAC is reading it when I told it to read the good file, if that makes sense.

Regards,
Sheila Peacock,
Blacknest.




  • Milton Plasencia
    2021-05-20 20:45:02
    Hi Sheila,

    My two cents,

    I read the help for the transfer command..i recommended you read it.

    To use evalresp and specify a particular file you must use "fname" to call
    a specific resp file

    SAC> transfer to evalresp fname RESP.DG.0385F.01.HHZ

    Station (-12345 ), Channel (-12345 )

    Using response from user-specified file: (RESP.DG.0385F.01.HHZ). "Here the
    use resp file"

    Extracting evresp response for *, *...

    if you do not use fname, evalresp "scanned" the dir and perhaps use another
    resp file...

    I don't know why evalresp works with "subtype" . It is works for polezero
    option, I need to read one more time for the transfer help...

    One of the resp files has the % symbol , and it is not recognized by the
    command transfer, it needs a physical unit..

    Cheers,

    Milton




    ************************************
    Milton P. PLASENCIA LINARES
    Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
    Borgo Grotta Gigante 42/C
    (34010) Sgonico - Trieste - Italia
    Tel: +39 040 2140 156 (Udine)
    Tel: +39 040 2140 256 (Trieste)
    Cell.: +39 331 6481 935

    E-mail: mplasencia<at>inogs.it

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


    On Thu, 20 May 2021 at 18:49, Sheila Peacock (via IRIS) <
    sac-help-bounce<at>lists.ds.iris.edu> wrote:

    Dear all,

    I extracted data from a SEED file with rdseed and told it to write the
    RESP files so I could deconvolve the instrument response.

    Many of the channels in the SEED file aren't seismic and their RESP files
    are not useful, but when I tested the response of one of the channels with
    an impulse, I got an error message as if the "transfer" subprogram had been
    reading from the wrong file:

    SAC> funcgen impulse npts 16384 delta 0.004
    SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
    Station (-12345 ), Channel (-12345 )
    Extracting evresp response for *, *...
    Thu May 20 17:26:26 2021 [WARN] check_units; units found ('% -
    PERCENTAGE') are not supported
    Thu May 20 17:26:26 2021 [ERROR] DG.0385F.01.LZ0: undefined units
    No transfer function applied!
    ERROR 2118: No transfer function applied.

    My test directory contains two RESP files:

    RESP.DG.0385F.01.HHZ.good RESP.DG.0385F.01.LZ0.bad

    If I delete the "bad" file or if I rename it e.g. to
    badresp.DG.0385F.01.LZO.bad then the problem goes away:

    SAC> funcgen impulse npts 16384 delta 0.004
    SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
    Station (-12345 ), Channel (-12345 )
    Extracting evresp response for *, *...

    What worries me is that if SAC is not using the RESP file that I told it
    to use, then if I had several "good" ones in my directory, how can I be
    sure it used the one I specified in the "transfer "command?

    I've attached the "good" and "bad" files. I'm not interested in why the
    bad file is bad, just why SAC is reading it when I told it to read the good
    file, if that makes sense.

    Regards,
    Sheila Peacock,
    Blacknest.



    ----------------------
    SAC Help
    Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
    Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu

    Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
    Update subscription preferences at http://ds.iris.edu/account/profile/


    • Sheila Peacock
      2021-05-21 08:07:41
      Milton - thanks - my mistake assuming that the TRANSFER type EVALRESP used the same command option as POLEZERO.

      Sheila.


      On 20 May 2021, at 19:45, Milton Plasencia <mplasencia<at>inogs.it> wrote:

      Hi Sheila,

      My two cents,

      I read the help for the transfer command..i recommended you read it.

      To use evalresp and specify a particular file you must use "fname" to call a specific resp file

      SAC> transfer to evalresp fname RESP.DG.0385F.01.HHZ
      Station (-12345 ), Channel (-12345 )
      Using response from user-specified file: (RESP.DG.0385F.01.HHZ). "Here the use resp file"
      Extracting evresp response for *, *...

      if you do not use fname, evalresp "scanned" the dir and perhaps use another resp file...

      I don't know why evalresp works with "subtype" . It is works for polezero option, I need to read one more time for the transfer help...

      One of the resp files has the % symbol , and it is not recognized by the command transfer, it needs a physical unit..

      Cheers,

      Milton




      ************************************
      Milton P. PLASENCIA LINARES
      Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
      Borgo Grotta Gigante 42/C
      (34010) Sgonico - Trieste - Italia
      Tel: +39 040 2140 156 (Udine)
      Tel: +39 040 2140 256 (Trieste)
      Cell.: +39 331 6481 935

      E-mail: mplasencia<at>inogs.it <mplasencia<at>inogs.it>

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


      On Thu, 20 May 2021 at 18:49, Sheila Peacock (via IRIS) <sac-help-bounce<at>lists.ds.iris.edu <sac-help-bounce<at>lists.ds.iris.edu>> wrote:
      Dear all,

      I extracted data from a SEED file with rdseed and told it to write the RESP files so I could deconvolve the instrument response.

      Many of the channels in the SEED file aren't seismic and their RESP files are not useful, but when I tested the response of one of the channels with an impulse, I got an error message as if the "transfer" subprogram had been reading from the wrong file:

      SAC> funcgen impulse npts 16384 delta 0.004
      SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
      Station (-12345 ), Channel (-12345 )
      Extracting evresp response for *, *...
      Thu May 20 17:26:26 2021 [WARN] check_units; units found ('% - PERCENTAGE') are not supported
      Thu May 20 17:26:26 2021 [ERROR] DG.0385F.01.LZ0: undefined units
      No transfer function applied!
      ERROR 2118: No transfer function applied.

      My test directory contains two RESP files:

      RESP.DG.0385F.01.HHZ.good RESP.DG.0385F.01.LZ0.bad

      If I delete the "bad" file or if I rename it e.g. to badresp.DG.0385F.01.LZO.bad then the problem goes away:

      SAC> funcgen impulse npts 16384 delta 0.004
      SAC> transfer to evalresp subtype RESP.DG.0385F.01.HHZ.good
      Station (-12345 ), Channel (-12345 )
      Extracting evresp response for *, *...

      What worries me is that if SAC is not using the RESP file that I told it to use, then if I had several "good" ones in my directory, how can I be sure it used the one I specified in the "transfer "command?

      I've attached the "good" and "bad" files. I'm not interested in why the bad file is bad, just why SAC is reading it when I told it to read the good file, if that makes sense.

      Regards,
      Sheila Peacock,
      Blacknest.



      ----------------------
      SAC Help
      Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>

      Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
      Update subscription preferences at http://ds.iris.edu/account/profile/


15:53:50 v.01697673