Thread: question about reading/writing sac files with fortran

Started: 2021-12-08 16:31:13
Last activity: 2021-12-10 14:44:56
Topics: SAC Help
Hello,

I recently downloaded the source code for SAC and I am interested in
reading and writing spectral amplitude (and phase) information from
seismogram time series with Fortran 90.

In the included examples, the program rsac1f.f makes reference to the
subroutine 'rsac1', but I encounter compilation errors because the
subroutine rsac1 can't be located. Do you know where I can find the f90
source code for this function? I hope I am posting my question to the right
place. Thanks so much!

Marlon



--

Marlon D. Ramos, Ph.D.

he/him/his

National Research Council (NRC) Postdoctoral Fellow

Air Force Research Lab, Kirtland AF Base

Albuquerque, NM

  • In SAC, "rsac1" is a SAC command, not a subroutine.   File ~/examples/README 
    shows how one can write programs in Fortran and C to call SAC commands, like rsac1.

    On 12/9/21 1:22 PM, Marlon Ramos (via IRIS) wrote:
    Hello,

    I recently downloaded the source code for SAC and I am interested in reading
    and writing spectral amplitude (and phase) information from seismogram time
    series with Fortran 90.

    In the included examples, the program rsac1f.f makes reference to the
    subroutine 'rsac1', but I encounter compilation errors because the subroutine
    rsac1 can't be located. Do you know where I can find the f90 source code for
    this function? I hope I am posting my question to the right place. Thanks so
    much!

    Marlon



    --

    Marlon D. Ramos, Ph.D.

    he/him/his

    National Research Council (NRC) Postdoctoral Fellow

    Air Force Research Lab, Kirtland AF Base

    Albuquerque, NM



  • George Helffrich
    2021-12-10 10:27:54
    Dear Marlon -

    Please see the online resources provided with the SAC book (The Seismic Analysis Code, ISBN 9781107613195, Cambridge Univ. Press). Among them are Fortran 90 routines to read and write SAC files.

    On 9 Dec 2021, at 18:22, Marlon Ramos (via IRIS) <sac-help-bounce<at>lists.ds.iris.edu> wrote:

    Hello,

    I recently downloaded the source code for SAC and I am interested in reading and writing spectral amplitude (and phase) information from seismogram time series with Fortran 90.

    In the included examples, the program rsac1f.f makes reference to the subroutine 'rsac1', but I encounter compilation errors because the subroutine rsac1 can't be located. Do you know where I can find the f90 source code for this function? I hope I am posting my question to the right place. Thanks so much!

    Marlon



    --
    Marlon D. Ramos, Ph.D.
    he/him/his
    National Research Council (NRC) Postdoctoral Fellow
    Air Force Research Lab, Kirtland AF Base
    Albuquerque, NM


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


16:14:28 v.01697673