calculating_Green_functions
所属分类:数学计算
开发工具:Fortran
文件大小:1606KB
下载次数:55
上传日期:2012-10-31 20:01:31
上 传 者:
wuyu19880216
说明: 理论地震学:平行层状介质理论地震图计算源程序
(Theoretical seismology: the parallel layered media seismograms source)
文件列表:
Programs_for_calculating_Green_functions\Bessel.f (5053, 1999-08-07)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\Draw_Seismogram_Apsel10_11.m (831, 2001-12-22)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\fig5_10.m (764, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\fig5_11.m (771, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\Inbasic.dat (1230, 2001-12-22)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\Inbasic.dat.bak (1230, 2001-12-21)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\inmedia.dat (375, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\In_r_r.dat (331, 2001-12-21)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\In_r_r.dat.bak (332, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\In_s_d.dat (365, 2001-12-21)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\In_s_d.dat.bak (367, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\Uf.dat (17920, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.10~5.11\Ur.dat (17920, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\Draw_Seismogram_Apsel12_13.m (831, 2001-12-22)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\fig5_12.m (762, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\fig5_13.m (759, 2000-10-20)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\Inbasic.dat (1230, 2001-12-22)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\inmedia.dat (375, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\In_r_r.dat (331, 2001-12-22)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\In_r_r.dat.bak (332, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\In_s_d.dat (367, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\Uf.dat (17920, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.12~5.13\Ur.dat (17920, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\fig5_14.m (772, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\Inbasic.dat (973, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\inmedia.dat (375, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\In_r_r.dat (332, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\In_s_d.dat (367, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\Uf.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\Ur.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.14\Uz.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\fig5_15.m (773, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\Inbasic.dat (973, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\inmedia.dat (375, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\In_r_r.dat (333, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\In_s_d.dat (367, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\Uf.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\Ur.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.15\Uz.dat (7168, 1999-09-14)
Programs_for_calculating_Green_functions\Comparisons_2\Apsel(1979)\figure5.16\fig5_16.m (773, 1999-09-14)
... ...
INSTRUCTION
This program is to calculate the synthesized seismograms in multilayered medium due to
an arbitary source(double couple, single force and explosive source). For detail, please
refer to Prof. Xiaofei Chen's paper(1999). This program is originally created by Prof.
Chen and modified by H.M.Zhang.
To run this program, you must create a workspace in FORTRAN POWERSTATION or Visual Fortran
including the all the *.f files. The head file 'green_com.inc' will be included automatically
when you insert the *.f files into the workspace.
Let us give an explanation for the input- and output-files:
INPUT FILES:
~~~~~~~~~~~~
(1). inmedia.dat:
In this file you set the model. you should set the total layer number and the depth of
interface, density, velocities and quality factors for S wave and P wave in each layer
respectively.
(2). inbasic.dat:
In this file you set several important control parameters. We'll give an explanation one
by one:
m ------------Power exponent of sample point number, i.e., the total sample point number
in time domain is MT=2^m.
Twin ---------Time window width, the time interval is [0, Twin].
criterion ----Criterion for k integral. It is only valid when the self-adaptive Filon's
integration method is used (r0>1000km). An appropriate value is 0.01.
SourceType ---Source type, 'D' or 'd' for double couple source, 'S' or 's' for single
force source, and 'E' or 'e' for explosive source.
Type_STF -----Type of source time function. You may add special source time function to
function 'sw' in gsubs.f.
rise time ----Rise time for source time function.
WindowSwitch -'Yes' for add spectral window and 'No' for else.
!!! CXF made modification: 'ON' instead of 'Yes', and 'OFF' instead of 'No' !!! (2001.12.21)
WindowType ---The type of the spectral window(filter), and there is only one alternative
iterm - Hamming - at the present time. It's active only when WindowSwitch
are set as 'Yes'.
f1,f2,f3,f4 --Parameters of Hamming window. They're active only when WindowSwitch are set
as 'Yes'.
OutDomain ----Parameter to control the domain of output data, alternative items are:
'time' and 'frequency'.
Case ---------Source-receiver configaration case: 'r-array' for receiver array case,
and 'vsp' for 'VSP' case.
(3). in-r-r.dat:
This data file is active only when parameter 'Case' in inbasic.dat is set to be 'r-array'.
r0 -----------Distance (in km) between epicenter and receiver.
fai0 ---------Azimuth angle (in degree) of the receiver array.
z0 -----------Depth (in km) of the receiver(s).
dz0 ----------Interval of the receivers in r-direction, it's active only when n_z != 1.
n_z ----------Number of receivers.
(4). in-r-vsp.dat:
This data file is active only when parameter 'Case' in inbasic.dat is set to be 'vsp'.
r0 -----------Distance (in km) between epicenter and receiver.
fai0 ---------Azimuth angle (in degree) of the receiver array.
z0 -----------Depth (in km) of the first receiver(s).
dz0 ----------Interval of the receivers in z-direction, it's active only when n_z != 1.
n_z ----------Number of receivers.
(5). in-s-d.dat:
This data file is active only when parameter 'SourceType ' in inbasic.dat is set to be 'D'
or 'd'.
D0s ----------Slip amount (in km).
A0s ----------Size of fault plane (in km^2).
dip ----------Dip angle (in degree).
rake ---------Rake angle (in degree).
strike -------Stirke angle (in degree).
zs -----------Source depth (in km).
(6). in-s-s.dat:
This data file is active only when parameter 'SourceType ' in inbasic.dat is set to be 'S'
or 's'.
s0 -----------Intensity of the single-point source.
nx -----------X-direction component of unit vector in direction of the force.
ny -----------Y-direction component of unit vector in direction of the force.
nz -----------Z-direction component of unit vector in direction of the force.
zs -----------Source depth (in km).
(7). in-s-e.dat:
This data file is active only when parameter 'SourceType ' in inbasic.dat is set to be 'E'
or 'e'.
Me -----------Source moment.
zs -----------Source depth (in km).
OUTPUT FILES:
~~~~~~~~~~~~~
(1). Ur.dat:
Store the radial components of synthesized seismogram (in time domain or freqency domain).
a.When parameter 'OutDomain' in inbasic.dat is set to be 'time', There are mt lines and
(jz+1) columns in Ur.dat, where mt is the sample point number in time domain, and jz is
the number of receivers. The first column is time, and the second to (jz+1)th column are
radial components in time domain recorded at different receivers, respectively.
b.When parameter 'OutDomain' in inbasic.dat is set to be 'frequency', There are nf2 lines
and (1+2*jz) columns in Ur.dat, where nf2=mt/2. The first column is frequency, and the
second to (1+2*jz)th column are radial components in freqency domain recorded at different
receivers, respectively.
(2). Uf.dat:
It is the same as ur.dat, except for transverse component.
(3). Uz.dat:
Also the same as ur.dat, except for vertical component.
(4). time.dat:
Record the run time of this program.
NOTE:
~~~~~
1. This program is still need some modification. At present time, only the case of
DOUBLE COUPLE can be implemented and jz must be assigned to 1. The user may make
corresponding modification as you need.
2. Although many examples have been calculated to validate this program, we are not
sure at the present time that this program will give accurate and reliable result
for all cases you encounter, because this program still need for more validating
and developing works.
3. Please don't hesitate to contact with us and discuss any question you encounter
in use.
April,2000.
近期下载者:
相关文件:
收藏者: