QSSP2SAC

Convert synthetic seismograms calculated by QSSP to SAC file format.

Features:

1. Convert outputs of QSSP to SAC-format conveniently with one command or python script;
2. Write basic header information like locations of event and stations, names of event and stations, and channel names;
3. Highlight observable types and directions according to QSSP outputs in the name of SAC files.

Python Script

Prerequisites

Module sactrace of obspy.io.sac is used for SAC writing, so the tool box ObsPy (Beyreuther et al., 2010) is needed. See ObsPy homepage for the instructions of installation.

Usage:

$ python qssp2sac.py qssp.inp qssp.out

Note

• qssp.inp and qssp.out are in the SAME directory of qssp2sac.py, in the other word, they should be in the current working directory; (using absolute or relative path will be added in the furture)
• qssp.inp is the input file with parameters for QSSP;
• qssp.out is one of several output files after running QSSP with the input file qssp.inp;
• the file name of qssp.out should follow the format as prefix_observable_direction.dat, in which prefix is set in section RECEIVER PARAMETERS in qssp.inp, observable is the one you select in section RECEIVER PARAMETERS in qssp.inp and direction is from e/n/z or ee/en/ez/nn/nz/zz upon different observables.
• JUST REMEMBER NOT TO RENAME OUTPUT FILES OF QSSP WHEN USING THIS SCRIPT.

Example:

Assuming you run QSSP with the input data file qssp_exmaple.inp and obtain 3 outputs as example_disp_e.dat example_disp_n.dat, and example_disp_z.dat (these four files could be found in folder ./Example). Then you could run the command as

$ cd ./Example
$ mv ../qssp2sac.py ./
$ python qssp2sac.py qssp_example.inp exmaple_disp_e.dat
QSSP input file is qssp_example.inp
QSSP output file is example_disp_e.dat
Create SAC file: EXAMPLE.STA1.DISP.BXE.SAC
Create SAC file: EXAMPLE.STA2.DISP.BXE.SAC

and get two SAC files as EXAMPLE.STA1.DISP.BXE.SAC and EXAMPLE.STA2.DISP.BXE.SAC. You could check header information and plot waveform using program SAC (Seismic Analysis Code).

Reference

M. Beyreuther, R. Barsch, L. Krischer, T. Megies, Y. Behr and J. Wassermann (2010) ObsPy: A Python Toolbox for Seismology SRL, 81(3), 530-533 DOI: 10.1785/gssrl.81.3.530

Fortran

Compilation

A makefile called makefile can be used to compile your code.

$ cd ./FortranSrc/
$ make qssp2sac

Otherwise, you can compile the code straightforwardly as

$ cd ./FortranSrc
$ ifort *.f90 -o qssp2sac

Usage

After compilation, an executable program is created (e.g. qssp2sac). THen you can run it as any other executable programs using ./:

You need to type the name of input file for QSSP, which contains all parameters to synthesize the seismograms with QSSP, and the name of output file of QSSP, which contains the synthetic seismograms.

$ ./qssp2sac
 The input file for QSSP is:
qssp.inp
 The output seismograms of QSSP is:
qssp.out

Note

• qssp.inp and qssp.out are in the SAME directory as qssp2sac, in the other word, they should be in the current working directory; (using absolute or relative path will be added in the furture)
• qssp.inp is the input file with parameters for QSSP;
• qssp.out is one of several output files after running QSSP with the input file qssp.inp;
• the file name of qssp.out should follow the format as prefix_observable_direction.dat, in which prefix is set in section RECEIVER PARAMETERS in qssp.inp, observable is the one you select in section RECEIVER PARAMETERS in qssp.inp and direction is from e/n/z or ee/en/ez/nn/nz/zz upon different observables.
• JUST REMEMBER NOT TO RENAME OUTPUT FILES OF QSSP WHEN USING THIS SCRIPT.

Example

Assuming you run QSSP with the input data file qssp_exmaple.inp and obtain 3 outputs as example_disp_e.dat example_disp_n.dat, and example_disp_z.dat (these four files could be found in folder ./Example). Then you could run the command as

$ cd ./Example
$ mv ../FortranSrc/qssp2sac ./
$ ./qssp2sac
 The input file for QSSP is:
qssp_example.inp
 The output seismograms of QSSP is:
example_disp_e.dat
 Create SAC file: EXAMPLE.STA1.DISP.BXE.SAC
 Create SAC file: EXAMPLE.STA2.DISP.BXE.SAC

and get two SAC files as EXAMPLE.STA1.DISP.BXE.SAC and EXAMPLE.STA2.DISP.BXE.SAC. You could check header information and plot waveform using program SAC (Seismic Analysis Code).