Thread: Plotrecordsection problem

Started: 2020-02-24 16:24:13
Last activity: 2020-02-25 12:52:10
Topics: SAC Help
Juan Rueda
Plotrecordsection problem
2020-02-24 16:24:13
Hello.
I have a problem with the plotrecorsection command. I can't draw the travel-times on the seismograms, even with the simple example code.
I do not understand what's happening.
Thanks

SAC> datagen sub regional *.z
/usr/local/sac/aux/datagen/regional/elk.z ...knb.z ...lac.z ...mnv.z
SAC> sss
Signal Stacking Subprocess.
SAC/SSS> traveltime model iasp91 phase Pn Sn
traveltime: depth: 5.000000 km
traveltime: Pn at 57.169170 s [ t = 57.169170 s ]
traveltime: Sn at 101.646294 s [ t = 101.646294 s ]
traveltime: Pn at 66.056442 s [ t = 66.056442 s ]
traveltime: Sn at 117.625183 s [ t = 117.625183 s ]
traveltime: Pn at 58.111580 s [ t = 58.111580 s ]
traveltime: Sn at 103.340820 s [ t = 103.340820 s ]
traveltime: error finding phase Pn
traveltime: error finding phase Sn
SAC/SSS> prs ttime on orient portrait
  • George Helffrich
    Re: Plotrecordsection problem
    2020-02-24 19:04:28
    Dear JJrueda -

    Provided that you define a time window, this works properly with SAC/BRIS, e.g.

    SAC> datagen sub regional *.z
    elk.z knb.z lac.z mnv.z
    SAC> sss
    Signal Stacking Subprocess.
    SAC/SSS> traveltime model iasp91 phase Pn Sn
    SAC/SSS> timewindow -10 300
    SAC/SSS> prs ttime on orient portrait

    As the messages indicate, there may be problems with the traveltime command in SAC2000, however.

    On 24 Feb 2020, at 16:24, jjrueda<at>fomento.es wrote:

    Hello.
    I have a problem with the plotrecorsection command. I can't draw the travel-times on the seismograms, even with the simple example code.
    I do not understand what's happening.
    Thanks

    SAC> datagen sub regional *.z
    /usr/local/sac/aux/datagen/regional/elk.z ...knb.z ...lac.z ...mnv.z
    SAC> sss
    Signal Stacking Subprocess.
    SAC/SSS> traveltime model iasp91 phase Pn Sn
    traveltime: depth: 5.000000 km
    traveltime: Pn at 57.169170 s [ t = 57.169170 s ]
    traveltime: Sn at 101.646294 s [ t = 101.646294 s ]
    traveltime: Pn at 66.056442 s [ t = 66.056442 s ]
    traveltime: Sn at 117.625183 s [ t = 117.625183 s ]
    traveltime: Pn at 58.111580 s [ t = 58.111580 s ]
    traveltime: Sn at 103.340820 s [ t = 103.340820 s ]
    traveltime: error finding phase Pn
    traveltime: error finding phase Sn
    SAC/SSS> prs ttime on orient portrait

    ----------------------
    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/



    George Helffrich
    george<at>elsi.jp


    • Juan Rueda
      RE: Plotrecordsection problem
      2020-02-24 19:54:20
      Dear George.
      Unfortunately, with SAC/IRIS 101.6a I don't get the same result of figure 8.2 of your book. Only plot the seismograms not the travel-time.
      Thanks
      Juan

      ________________________________
      De: sac-help-bounce<at>iris.washington.edu <sac-help-bounce<at>iris.washington.edu> en nombre de George Helffrich <george<at>elsi.jp>
      Enviado: lunes, 24 de febrero de 2020 20:05
      Para: SAC Help <sac-help<at>lists.ds.iris.edu>
      Asunto: Re: [IRIS][sac-help] Plotrecordsection problem

      Dear JJrueda -

      Provided that you define a time window, this works properly with SAC/BRIS, e.g.

      SAC> datagen sub regional *.z
      elk.z knb.z lac.z mnv.z
      SAC> sss
      Signal Stacking Subprocess.
      SAC/SSS> traveltime model iasp91 phase Pn Sn
      SAC/SSS> timewindow -10 300
      SAC/SSS> prs ttime on orient portrait

      As the messages indicate, there may be problems with the traveltime command in SAC2000, however.

      On 24 Feb 2020, at 16:24, jjrueda<at>fomento.es<jjrueda<at>fomento.es> wrote:

      Hello.
      I have a problem with the plotrecorsection command. I can't draw the travel-times on the seismograms, even with the simple example code.
      I do not understand what's happening.
      Thanks

      SAC> datagen sub regional *.z
      /usr/local/sac/aux/datagen/regional/elk.z ...knb.z ...lac.z ...mnv.z
      SAC> sss
      Signal Stacking Subprocess.
      SAC/SSS> traveltime model iasp91 phase Pn Sn
      traveltime: depth: 5.000000 km
      traveltime: Pn at 57.169170 s [ t = 57.169170 s ]
      traveltime: Sn at 101.646294 s [ t = 101.646294 s ]
      traveltime: Pn at 66.056442 s [ t = 66.056442 s ]
      traveltime: Sn at 117.625183 s [ t = 117.625183 s ]
      traveltime: Pn at 58.111580 s [ t = 58.111580 s ]
      traveltime: Sn at 103.340820 s [ t = 103.340820 s ]
      traveltime: error finding phase Pn
      traveltime: error finding phase Sn
      SAC/SSS> prs ttime on orient portrait

      ----------------------
      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/



      George Helffrich
      george<at>elsi.jp<george<at>elsi.jp>


      • Geo .
        Re: Plotrecordsection problem
        2020-02-25 12:52:10
        Hi Juan,
        Try the following and Pn and Sn should show up except for the seismogram
        closest to the source:

        SAC> datagen sub regional *.z

        /usr/local/sac/aux/datagen/regional/elk.z ...knb.z ...lac.z ...mnv.z

        SAC> sss

        Signal Stacking Subprocess.

        SAC/SSS> traveltime model iasp91 picks 0 phase Pn Sn

        traveltime: depth: 5.000000 km

        SAC/SSS> prs ttime on

        However, note that you are writing Pn and Sn time headers to Tn starting
        from T0 now.

        Januka


        Januka Attanayake
        Research Fellow | Earthquake Seismology
        Homepage: *https://sites.google.com/site/janukaattanayake
        https://sites.google.com/site/janukaattanayake*

        School of Earth Sciences | McCoy Bldg. 200
        University of Melbourne | Parkville 3010 VIC
        Australia




        On Tue, 25 Feb 2020 at 06:56, Juan Rueda <jjrueda<at>fomento.es> wrote:

        Dear George.
        Unfortunately, with SAC/IRIS 101.6a I don't get the same result of figure
        8.2 of your book. Only plot the seismograms not the travel-time.
        Thanks
        Juan

        ------------------------------
        *De:* sac-help-bounce<at>iris.washington.edu <
        sac-help-bounce<at>iris.washington.edu> en nombre de George Helffrich <
        george<at>elsi.jp>
        *Enviado:* lunes, 24 de febrero de 2020 20:05
        *Para:* SAC Help <sac-help<at>lists.ds.iris.edu>
        *Asunto:* Re: [IRIS][sac-help] Plotrecordsection problem

        Dear JJrueda -

        Provided that you define a time window, this works properly with
        SAC/BRIS, e.g.

        SAC> datagen sub regional *.z
        elk.z knb.z lac.z mnv.z
        SAC> sss
        Signal Stacking Subprocess.
        SAC/SSS> traveltime model iasp91 phase Pn Sn
        SAC/SSS> timewindow -10 300
        SAC/SSS> prs ttime on orient portrait

        As the messages indicate, there may be problems with the traveltime
        command in SAC2000, however.

        On 24 Feb 2020, at 16:24, jjrueda<at>fomento.es wrote:

        Hello.
        I have a problem with the plotrecorsection command. I can't draw the
        travel-times on the seismograms, even with the simple example code.
        I do not understand what's happening.
        Thanks

        SAC> datagen sub regional *.z
        /usr/local/sac/aux/datagen/regional/elk.z ...knb.z ...lac.z ...mnv.z
        SAC> sss
        Signal Stacking Subprocess.
        SAC/SSS> traveltime model iasp91 phase Pn Sn
        traveltime: depth: 5.000000 km
        traveltime: Pn at 57.169170 s [ t = 57.169170 s ]
        traveltime: Sn at 101.646294 s [ t = 101.646294 s ]
        traveltime: Pn at 66.056442 s [ t = 66.056442 s ]
        traveltime: Sn at 117.625183 s [ t = 117.625183 s ]
        traveltime: Pn at 58.111580 s [ t = 58.111580 s ]
        traveltime: Sn at 103.340820 s [ t = 103.340820 s ]
        traveltime: error finding phase Pn
        traveltime: error finding phase Sn
        SAC/SSS> prs ttime on orient portrait

        ----------------------
        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/




        George Helffrich
        george<at>elsi.jp


        ----------------------
        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/


NSF SAGE Facility has been certified as a trustworthy data repository by the CoreTrustSeal Standards and Certification Board
08:49:54 v.b3198453