wphase:documentation
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
|
wphase:documentation [2014/03/19 01:03] wphase |
wphase:documentation [2015/01/23 17:16] wphase |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ===== W-phase documentation ===== | + | ====== W-phase documentation ====== |
| - | ==== Installation ==== | + | |
| + | ===== Installation ===== | ||
| + | |||
| + | ==== Getting the code ==== | ||
| + | |||
| + | Currently, the W-phase package is hosted as a github repository. Ask Zacharie Duputel <zacharie.duputel@unistra.fr> or Luis Rivera <luis.rivera@unistra.fr> for access to the repository. | ||
| + | |||
| + | Check out W-phase repository: | ||
| + | <code> | ||
| + | git clone https://github.com/eost/wphase.git wphase_package | ||
| + | </code> | ||
| + | |||
| + | To update your W-phase repository (pull changes) | ||
| + | <code> | ||
| + | cd /to/the/wphase/directory/ | ||
| + | git pull origin master | ||
| + | </code> | ||
| + | |||
| + | |||
| + | ==== Dependencies ==== | ||
| + | |||
| + | The w-phase package have only been tested on Unix and Linux computers. You will need the following: | ||
| + | - csh shell | ||
| + | - gcc and gfortran | ||
| + | - python2.7 | ||
| + | - You have to install pylab and basemap (matplotlib toolkit) to run some python scripts which make figures. | ||
| + | - GMT4 | ||
| + | |||
| + | ==== Building the code ==== | ||
| To install the code, we must first setup a few environment variables. If you use csh or tcsh: | To install the code, we must first setup a few environment variables. If you use csh or tcsh: | ||
| Line 38: | Line 66: | ||
| ----- | ----- | ||
| - | ==== How to run the W-phase code ==== | + | ===== How to run W-phase ===== |
| - | === Preparing directories, i_master and CMTSOLUTION === | + | ==== Preparing directories, i_master and CMTSOLUTION ==== |
| Before each inversion, it is necessary to create a “run” directory containing two input ascii files: | Before each inversion, it is necessary to create a “run” directory containing two input ascii files: | ||
| * CMTSOLUTION: a file containing the PDE, the event centroid location and timing and optionally a “reference” moment tensor solution. | * CMTSOLUTION: a file containing the PDE, the event centroid location and timing and optionally a “reference” moment tensor solution. | ||
| Line 51: | Line 79: | ||
| ----- | ----- | ||
| - | === Extracting data from (mini)SEED file(s) === | + | ==== Extracting data from SEED ==== |
| Once the i_master and CMTSOLUTION files are created, we can extract waveforms and instrument response parameters and perform a rough screening by epicentral distance. This can be done using: | Once the i_master and CMTSOLUTION files are created, we can extract waveforms and instrument response parameters and perform a rough screening by epicentral distance. This can be done using: | ||
| <code>${WPHASE_HOME}/bin/extract.csh</code> | <code>${WPHASE_HOME}/bin/extract.csh</code> | ||
| Line 60: | Line 88: | ||
| ----- | ----- | ||
| - | === Calculating Synthetics, deconvolution and filtering === | + | ==== Calculating Synthetics, deconvolution and filtering ==== |
| The next step is to calculate the kernel functions associated with the 6 elements of the seismic moment tensor for the stations listed in 'scr_dat_fil_list' and to convolve them with the moment rate function (MRF) specified in the CMTSOLUTION file (time shift and half duration). We must then deconvolve the instrument response from the data and band pass filter each waveform in the frequency pass band specified in i_master. This is performed using: | The next step is to calculate the kernel functions associated with the 6 elements of the seismic moment tensor for the stations listed in 'scr_dat_fil_list' and to convolve them with the moment rate function (MRF) specified in the CMTSOLUTION file (time shift and half duration). We must then deconvolve the instrument response from the data and band pass filter each waveform in the frequency pass band specified in i_master. This is performed using: | ||
| <code>${WPHASE_HOME}/bin/prepare_wp.csh</code> | <code>${WPHASE_HOME}/bin/prepare_wp.csh</code> | ||
| Line 70: | Line 98: | ||
| ----- | ----- | ||
| - | === Inversion === | + | ==== Inversion ==== |
| The inversion can then be performed using the Kernels functions in ./GF and data files listed in i_wpinversion. | The inversion can then be performed using the Kernels functions in ./GF and data files listed in i_wpinversion. | ||
| Line 105: | Line 133: | ||
| ----- | ----- | ||
| - | === RUNA scripts === | + | ==== RUNA scripts ==== |
| Data extraction, screening and inversion described above can be performed by running one of the RUNA3*csh scripts (usually in bin). These scripts perform data extraction/screening and calculate a moment tensor solution after median data screening and after rejecting traces associated with large misfit using the threshold 5.0 3.0 0.9 (i.e. -th, see 2. in section 4): | Data extraction, screening and inversion described above can be performed by running one of the RUNA3*csh scripts (usually in bin). These scripts perform data extraction/screening and calculate a moment tensor solution after median data screening and after rejecting traces associated with large misfit using the threshold 5.0 3.0 0.9 (i.e. -th, see 2. in section 4): | ||
| <code> | <code> | ||
| Line 135: | Line 163: | ||
| ----- | ----- | ||
| - | === Grid searches === | + | ==== Grid searches ==== |
| In the grid-search scheme, there is a first global rough exploration which is followed by finer samplings around minimal points. If the optimum is found near a bound, the explored space is extended. | In the grid-search scheme, there is a first global rough exploration which is followed by finer samplings around minimal points. If the optimum is found near a bound, the explored space is extended. | ||
| Line 146: | Line 174: | ||
| ----- | ----- | ||
| - | === Plot routines === | + | ==== Plot routines ==== |
| All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module basemap is also needed for plotting maps but it is optional (even if we recommend to install it). | All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module basemap is also needed for plotting maps but it is optional (even if we recommend to install it). | ||
| Line 159: | Line 187: | ||
| - | ==== File formats ==== | + | ===== Data formats ===== |
| === CMTSOLUTION FILE === | === CMTSOLUTION FILE === | ||
wphase/documentation.txt · Last modified: 2022/01/10 07:39 by wphase
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
