m3csdirest
Multiple 3-Component Station DIRection ESTimation (m3csdirest)
Science Score: 67.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
✓DOI references
Found 7 DOI reference(s) in README -
✓Academic publication links
Links to: zenodo.org -
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (10.9%) to scientific vocabulary
Repository
Multiple 3-Component Station DIRection ESTimation (m3csdirest)
Basic Info
- Host: GitHub
- Owner: stevenjgibbons
- License: gpl-3.0
- Language: Fortran
- Default Branch: main
- Size: 1.07 MB
Statistics
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 1
Metadata Files
README.md
m3csdirest
Multiple 3-Component Station DIRection ESTimation (m3csdirest)
Release v1.0.0 is permanently stored on Zenodo with DOI 10.5281/zenodo.8270726
https://doi.org/10.5281/zenodo.8270726
m3csdirest is a FORTRAN77 program that tries to estimate the direction of incoming P-waves on 3-component stations. It reads in SAC-format files and so you need to have SAC available in order to compile it in its current form. If you are not able to obtain SAC for whatever reason, you will need to replace the lines which call the subroutine RSAC1 with some other routine which fills the variables RZVALS, NPTS, BEG, and DT. It also requires the LAPACK and BLAS libraries for linear algebra routines that interpolate to find the optimal backazimuth.
The program should compile by typing
make m3csdirest
in the directory created once you have edited the makefile.
This program writes out ASCII output and is not intended for operational use in any routine system. It is designed as a proof of concept program.
Example
(1) Go into the directory KEVNZ20101011 and run the python script getKEV.py
This should generate 3 SAC-files: IUKEVBHE10.sac IUKEVBHN10.sac IUKEVBHZ_10.sac
(2) Assuming that you have managed to compile the program m3csdirest and have placed it in your path, run the script runKEVexample.sh
```
!/bin/sh
Note that the KEV data collected have 40 Hz sampling
so 4 seconds is 160 samples.
This script should calculate the output for the files
generated in the directory KEV_NZ20101011
inputfile=KEVm3csdirest.input cat << EOF > ${inputfile} KEVNZ20101011/IUKEVBHZ10.sac KEV KEVNZ20101011/IUKEVBHE10.sac 90.0 KEVNZ20101011/IUKEVBHN_10.sac 0.0 EOF NSTAT=1 NCCLEN=160 NSKIP=40 NAZIB2=36 ICFLAG=1 OUTFLAG=1111 m3csdirest $NSTAT $NCCLEN $NSKIP $NAZIB2 $ICFLAG $OUTFLAG < $inputfile > KEV.out ```
Auxiliary program - aziwsa
This program just tries to make a weighted average of azimuith values.
An example is done by running
sh run_aziwsa.sh
python plot_aziwsa_out.py
Relevant paper
If you use the code, I would like you to refer to the following publication:
Gibbons, S.J. (2023),
Direction Estimates for Short‐Period P‐Waves on Three‐Component Stations and Arrays,
The Seismic Record. 3(4), pp. 299–310
https://doi.org/10.1785/0320230036
!
Owner
- Name: Steven J. Gibbons
- Login: stevenjgibbons
- Kind: user
- Location: Oslo, Norway
- Company: NGI
- Website: https://scholar.google.no/citations?user=dO6xNCQAAAAJ
- Repositories: 2
- Profile: https://github.com/stevenjgibbons
I am an applied geoscientist with expertise in seismology, acoustics, computational fluid dynamics, signal processing, inverse theory, and natural hazards.
Citation (CITATION.cff)
cff-version: 1.2.0 message: "If you use this software, please cite it as below." authors: - family-names: Gibbons given-names: Steven J. orcid: "https://orcid.org/0000-0002-7822-0244" title: "m3csdirest" version: v1.0.0 doi: https://doi.org/10.5281/zenodo.8270726 date-released: 2023-08-21 url: "https://github.com/stevenjgibbons/m3csdirest"