The original code was written by Brian Kennett. I have slightly changed the code so that it can be compiled successfully. I also introduce some Python functions and scripts to run the ERZSOL3 more flexible and easier.
Weijia Sun
Email: swj (at) mail.iggcas.ac.cn
Updated on 20 August 2017
- Change line 136 of qbessel.f from
>>> COMPLEX FUNCTION BESHS0*16(X, IFAIL)to
>>> COMPLEX FUNCTION BESHS0(X, IFAIL), so that the compilation could be successful.
- Change the output seismogram from unformatted to text format.
- Add some python scripts convert the wavefroms into SAC format.
Updated on 24 Feburary 2018
- Enlarge array size from 600 to 3600.
sac.user0 :
sac.user1 :
*.cmd : input parameter file
*.mod : velocity model
*.dst : range and azimuth file
*.wav : wavelet file
- XS : float, array_like
- XS -- Station position, ranges in km
- AZ : float, array_like
- Azimuth from North (X-axis) in degree
- CZC : character*4
- : "Z"
- CRC : character*4
- : "R"
- CTC : character*4
- : "T"
- DELT: float
- Time sampling interval in second
- NT : int
- Sampling number, must be FFT length, e.g., 1024, 2048, 4096
- PR : float, unit might be s/km
- Reduced slowness
- STMIN : float
- starting reduced time
- TCAL : float
- Reduced time, STMIN+PR*XS(JX)