#RF TOOLBOX MATLAB HELP PDF WINDOWS#
It first links SAC data in the local directory to the earthquake catalog automatically, and then cuts the waveform in proper time windows by calculating the theoretical arrival times with Taup Toolkit (Crotwell et al. 2008) developed under MATLAB environment by French scientist is an open-source GUI toolbox for shear-wave splitting analysis. 2013), the FuncLab toolbox for RF analysis (Eagar and Fouch 2012), and SplitLab toolbox for shear-wave splitting analysis (Wüstefeld et al. Thanks to the rapid development of computer science, more GUI-based seismological softwares bring much convenience and power efficiency in seismic data analysis, such as the AIMBAT for picking arrival times (Lou et al. However, users have to write many small codes to obtain the preprocessed data first, e.g., selecting teleseismic events, removing average and linear trend from the waveforms, filtering, and selecting phases in proper time windows. These codes generally import the preprocessed seismic data to calculate RFs from teleseismic radial and vertical components. Many open-source codes have been released for RF analysis in recent years, such as the CPS 330 system, the RF module under Python, and the process RFmatlab module under MATLAB. Zhu and Kanamori ( 2000) further proposed the H- k stacking method with a grid-search algorithm that determines the v P/ v S ratio and Moho depth ( H) with the obtained P-wave RFs at a station, which is very important in determining one-dimensional structure in the Earth. The routines of RF analysis include the following procedures: (a) selecting the proper events (e.g., epicentral distance, magnitude) from the earthquake catalog and matching them with seismic data in the local directory (b) preprocessing the waveforms, such as removing average shift of amplitude, removing the linear trend, and filtering the waveforms with certain bandpass filters (c) picking P-wave arrivals and (d) calculating P-wave RFs by deconvolving the radial component with the vertical component.
2002), which makes it one of the most popular methods in seismology. The RF method is stable and efficient in determining the crustal and upper mantle structures (e.g., Dueker and Sheehan 1997 Kind et al. The time delay between the Ps converted phase and direct P-wave is sensitive to the depth of the velocity discontinuities and the average S-wave velocity above it (Zhu and Kanamori 2000). P-wave RF deconvolves the radial component with the vertical component and obtains the Ps phase converted at velocity discontinuities (e.g., Moho) in the Earth. Receiver function (RF) analysis, proposed by Langston ( 1979), is one of routine tools to determine the one-dimensional structure of the crust and mantle. It should look like the image below.Seismology is becoming the most popular tool in studying the interior structure of the Earth thanks to the development of digital seismology and the rapidly accumulated seismic data in recent years. Step 5: Transfer the matlab_test_output.pdf to your computer and open it (see Uploading & Downloading Files for instructions on transferring files). Step 4: Upon successful completion of the job, you will find three new files in your directory: matlab_output.txt, matlab_error.txt, and matlab_test_output.pdf. Step 3: Submit the batch job by typing sbatch matlab_batch.sh
#RF TOOLBOX MATLAB HELP PDF CODE#
If your code can run on all 20 or 24 cores of a node, please add "#SBATCH -exclusive" to the script above, to request the entire node, in order to avoid overloading a compute node. If your matlab code can run multi-threaded, please request as many cores as you need with -ntasks=x instead of 1. This assumes a single-threaded matlab process. Step 2: Create a batch submission file called matlab_batch.sh and insert the following lines of code: Step 1: Create a sample Matlab file called matlab_test.m and insert the following lines of code:
0 kommentar(er)
