Sheila Peacock
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.
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.
Attachments
-
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/
-
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/
-
Milton Plasencia
