Thread: deconvolution

Started: 2016-10-13 19:11:24
Last activity: 2016-10-26 12:15:46
Topics: SAC Help
Roberto Cabieces Diaz
2016-10-13 19:11:24
Hi users;

I would like to know if someone could help me in this topic,

I have to deconvolve a sort of files (6 moths of files, file per day) and
I would like to use a macro to read , for example every Z component files
and deconvolve to Acceleration.

(using the response of the z channel that I got previously with rdseed)

There are a simple answer for that or I have to do it one by one,
This is because i have to study the noise of a station in this channel.

thank you for any response, Roberto

  • Milton Plasencia
    2016-10-26 12:15:46
    Hi Roberto,

    I do a simple macro, i tried it with a short sac file.

    in the macro use a do wild loop and process all files in directory “.” that end with “.SAC”,
    adapt it to your needs, for example put your Z files in a separate dir and adapt freqlimits to
    the band of interest.

    So, save the macro to a file and put it in the working dir,
    run sac and then
    sac> m macro.file

    -----------
    DO FILE WILD DIR . *.SAC
    r $FILE
    rmean
    transfer from polezero s YOUR_SACPZ_FILE to acc freqlimits 0.01 0.02 4 8
    w append .corr
    ENDDO
    message " all is done “
    ——————

    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 156 (Udine)
    Tel: +39 040 2140 256 (Trieste)
    Cel.: +39 331 6481 935

    E-mail: mplasencia<at>inogs.it
    GPG key: 22FCFFA8

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

    On Oct 26, 2016, at 02:19, Roberto Cabieces Diaz <rcabdia<at>gmail.com> wrote:

    Hi users;

    I would like to know if someone could help me in this topic,

    I have to deconvolve a sort of files (6 moths of files, file per day) and I would like to use a macro to read , for example every Z component files and deconvolve to Acceleration.

    (using the response of the z channel that I got previously with rdseed)

    There are a simple answer for that or I have to do it one by one,
    This is because i have to study the noise of a station in this channel.

    thank you for any response, Roberto

    ----------------------
    SAC Help (http://ds.iris.edu/message-center/topic/sac-help/)

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


14:02:27 v.28f29ea9