文件列表:

hk1.3\hk\knt.c (11687, 2010-02-05)
hk1.3\hk\k_stack.o (6688, 2008-08-28)
hk1.3\hk\Complex.o (1856, 2008-08-28)
hk1.3\hk\k_stack.c (7155, 2011-01-12)
hk1.3\hk\example\hk.grd.std (17888, 2008-08-28)
hk1.3\hk\example\pp.1.ti (4232, 2005-01-14)
hk1.3\hk\example\pp.0.ri (4232, 2005-01-14)
hk1.3\hk\example\pp.6.ri (4232, 2005-01-14)
hk1.3\hk\example\KUL.r (5432, 2006-09-07)
hk1.3\hk\example\pp.6.ti (4232, 2005-01-14)
hk1.3\hk\example\pp.0.ti (4232, 2005-01-14)
hk1.3\hk\example\pp.1.ri (4232, 2005-01-14)
hk1.3\hk\example\pp.4.ti (4232, 2005-01-14)
hk1.3\hk\example\hk.grd (17888, 2008-08-28)
hk1.3\hk\example\pp.18.ri (4232, 2005-01-14)
hk1.3\hk\example\pp.18.ti (4232, 2005-01-14)
hk1.3\hk\example\KUL.z (5432, 2006-09-07)
hk1.3\hk\example\KUL.t (5432, 2006-09-07)
hk1.3\hk\example\KUL.ri (4232, 2006-09-07)
hk1.3\hk\example\pp.4.ri (4232, 2005-01-14)
hk1.3\hk\example\KUL.ti (4232, 2006-09-07)
hk1.3\hk\rcvFn.c (3045, 2009-12-22)
hk1.3\hk\sacio.c (13512, 2010-04-01)
hk1.3\hk\decon.c (5064, 2007-04-16)
hk1.3\hk\sacio.o (6156, 2008-08-28)
hk1.3\hk\inversion.h (696, 2010-03-12)
hk1.3\hk\fft.o (9448, 2008-08-28)
hk1.3\hk\Complex.h (2247, 2009-12-22)
hk1.3\hk\k_stack (810776, 2008-08-28)
hk1.3\hk\fft.c (13170, 2009-12-22)
hk1.3\hk\matrix.c (1340, 1999-02-25)
hk1.3\hk\iter_decon.c (8116, 2008-03-06)
hk1.3\hk\Complex.c (1012, 2008-07-27)
hk1.3\hk\knt.h (2427, 2009-12-21)
hk1.3\hk\hk_plt.pl (1214, 2005-01-14)
hk1.3\hk\grid2d.o (2000, 2008-08-28)
hk1.3\hk\grdmin.o (2644, 2008-08-28)
hk1.3\hk\sac.h (14615, 2009-04-11)
hk1.3\hk\Makefile (557, 2008-03-06)
... ...

# Information for using receiver function stacking package # @(#)README 1.0 11/30/2006 Copyright (c) 1996-2006, L. Zhu (lupei@eas.slu.edu) Permission to use, copy, modify, and distribute this package and supporting documentation for any purpose without fee is hereby granted, provided that the above copyright notice appear in all copies, that both that copyright notice and this permission notice appear in supporting documentation. In case that you modify any codes in this package, you have to change the name of the modified code. You are welcome to send me a note about your modification and any suggestion. In case that you redistribute this package to others, please send me the contacting info (email addr. preferred) so that future updates of this package can be sent to them. 1. Content of the Package The package contains following main codes (some might be in a separate package): decon.c time and frequency domain deconvolution. iter_decon.c time-domain iterative deconvolution. rcvFn.c Forward computation of receiver functions. k_stack.c H-K stacking. hk_plt.pl a Perl script to plot results of H-k stacking. grdmin.c find the minimum on a 2D grid. ccpStack3D.c CCP stacking. 3DSlice.c Slice a 2D cross-section in a 3D CCP volume. 2. Installation To compile the codes see Makefile. Note that you need to have GMT installed to be able to compile some programs. Two related environment variables, GMT_INC and GMT_LIBS, are needed. On my computer running Suss 10.1, they are set to: GMT_INC = -I/opt/GMT4.0/include -I/opt/netcdf-3.5.0/Linux/include GMT_LIBS = -L/opt/GMT4.0/Linux/lib -lgmt -lpsl -L/opt/netcdf-3.5.0/Linux/lib -lnetcdf -lm -s 3. Usage run the command with -help to show the usage of the command. In general, your waveform data should be organized by events in the form of event_ID/sta.[r,t,z]. The ray parameter p needs to be stored in SAC header USER0. It is recommended that the first P arrival be picked and stored in SAC header A For the time window, I normally chose 10 sec before the P arrival and 80 sec after, with a sampling rate of 10 sample/sec. 3.1 Deconvolution Subdirectory example/ contains a three-component waveform record (KUL.[r,t,z]). You can use the following command ./iter_decon -F1/3/-5 -N100 -C-2/-10/80 -T0.1 example/KUL.z example/KUL.[r,t] to obtain receiver functions KUL.ri and KUL.ti. 3.2 H-K stacking The subdirectory also has several averaged receiver functions of different incidence angle (pp.*.ri) that you can use ./k_stack -R20/60/1.5/2.0 -I0.5/0.01 -Gexample/hk.grd example/pp.*.ri to produce H-k stack grd file "example/hk.grd" and the variance file "example/hk.grd.var". You can then use the following command to find the stacking maximum: ./grdmin -D example/hk.grd and plot the result by using ./grdmin -D example/hk.grd | ./hk_plt.pl > junk.ps 3.2 CCP stacking The following three commands produce the CCP image in Zhu et al (2006) GJI paper: ccpStack3D -R-50/50/-25/25/0/50 -C28.030/38.161 -A3 -I1/5/0.5 -S1/4/1 -E1 -Vwa.vs -Htop.layer -Gccp.3d *.ri 3DSlice ccp.3d -Gjunk.grd -R-50/50/1 -C28.030/38.161 -A3 grdimage ccp.grd -Cccp.cpt -JX5/2.5 -R-50/50/-50/0 -Ba10f2/a10g30f2:"Depth (km)":WSe -K -Y2 > ccp.ps 4. Reference Zhu and Kanamori (2000), JGR. Zhu (2000), EPSL. Zhu et al (2006), GJI. 5. Revision History Dec 1, 2006 put the first release on line 6 FAQs

近期下载者：

相关文件：

评论：[我要评论] [举报此文件]

收藏者：