Note: We found a bug in "grfault.f" and "grpoint.f" on Sep.22, 1995.
Errors occurred when you put several sources in different layers.
If you have the old codes before Sep.22, 1995, please change them.
Note: We made a newer Green function code for seismic source, which is
in /pub/hisada/grflt3. That code runs much faster than "grfault.f"
in this directory. In addition, you can compute normal mode solutions.
Please try it.
Yoshiaki Hisada (Sep., 1995)
************************************************
*** README (Jan. 13 1995) by Yoshiaki Hisada ***
************************************************
* "grpoint.f" is a FORTAN code to compute Green's function due to point
sources. "grpoint.in" is an example data for "grpoint.f".
Please check "grpoint.manual" for the details.
* "grfault.f" is a FORTAN code to compute Green's function due to dislocation
sources. "grfault.in" is an example data for "grfault.f".
Please check "grfault.manual" for the details.
* "grfft.f" is a FORTRAN code for FFT, which transforms the frequency domain
data by "grpoint.f" or "grfault.f" into the time domain waveforms.
The output data will be stored in "wave_dat". Note that the data in the
first line in "wave_dat" are "the number of the data", "the time
increment", "the minimum amplitude", and "the maximum amplitude".
The waveform data are stored from the second line.
* "wave.m" is a "m" file for Matlab to plot the waveforms from "wave_dat".
A postscript file named "wave.ps" will be created.
* For comments or questions, please e-mail to "hisada@cc.kogakuin.ac.jp",
or mail to
Yoshiaki Hisada
Dept. of Architecture, Kogakuin Univ.,
Nishi-Shinjuku 1-24-2, Tokyo 163-91, Japan
(FAX: 3-3340-0140).