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.

近期下载者

相关文件


收藏者