This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
wphase:tutorial [2015/01/23 19:51] wphase |
wphase:tutorial [2021/10/22 18:52] (current) wphase |
||
---|---|---|---|
Line 1: | Line 1: | ||
======= W-phase tutorial ======= | ======= W-phase tutorial ======= | ||
- | You may also refer to [[wphase:documentation|the W-phase documentation page]], where scripts and executables are described in [[wphase:documentation#How to run W-phase|this section]], file-formats are described [[wphase:documentation#Data formats|here]]. | + | You may also refer to [[wphase:documentation|the W-phase documentation page]], where scripts and executables are described in [[wphase:documentation#How to run W-phase|this section]], input and output file-formats are described [[wphase:documentation#Data formats|here]]. |
This tutorial page only provides a typical example of W phase inversion run. | This tutorial page only provides a typical example of W phase inversion run. | ||
- | ===== Example of run for the Mw8.8 2010 Maule earthquake===== | + | ===== Example of run for the 2010 Maule earthquake===== |
- | ==== Input files ==== | + | We consider here an example of run for the Mw8.8 2010 Maule earthquake. |
- | Let’s consider the Mw8.8 2010 Maule earthquake. | + | |
+ | ==== Input files ==== | ||
The first thing to do is to create i_master and CMTSOLUTION files. We can use as input the following CMTSOLUTION file | The first thing to do is to create i_master and CMTSOLUTION files. We can use as input the following CMTSOLUTION file | ||
<code> | <code> | ||
Line 27: | Line 27: | ||
Mtp: -1.190000e+28 | Mtp: -1.190000e+28 | ||
</code> | </code> | ||
- | (for more detailed information about file formats, please refer to [[wphase:documentation#Data formats|the "File formats" section of the W-phase documentation page]]). | ||
- | This file correspond to the Global CMT solution except that we fixed the centroid latitute, longitude and depth to the PDE preliminary hypocenter values and that the time-shift and half duration have been fixed a priori using some initial magnitude estimate. | + | This file correspond to the Global CMT solution except that we fixed the centroid latitute, longitude and depth to the PDE preliminary hypocenter values and that the time-shift and half duration have been fixed a priori using some initial magnitude estimate (for more detailed information about file formats, please refer to [[wphase:documentation#Data formats|the "File formats" section of the W-phase documentation page]]). |
The i_master file is the following : | The i_master file is the following : | ||
Line 107: | Line 106: | ||
</code> | </code> | ||
- | ==== Centroid grid-search ==== | + | ==== Final inversion after centroid grid-search ==== |
The previous solution was computed with a fixed location (PDE). It is also possible to determine an optimum centroid position and timing by running the wp_grid_search.py script: | The previous solution was computed with a fixed location (PDE). It is also possible to determine an optimum centroid position and timing by running the wp_grid_search.py script: | ||
Line 172: | Line 171: | ||
Each spatial and temporal location explored during the grid-search is detailed in grid_search_ts_out (time-shift grid search) grid_search_xyz_out (lat, lon, dep grid-search) and grid_search_xy_out (finer lat, lon grid-search which is performed at the optimum depth found during the lat, lon, dep global exploration). | Each spatial and temporal location explored during the grid-search is detailed in grid_search_ts_out (time-shift grid search) grid_search_xyz_out (lat, lon, dep grid-search) and grid_search_xy_out (finer lat, lon grid-search which is performed at the optimum depth found during the lat, lon, dep global exploration). | ||
- | === Plotting the results === | + | ==== Plotting the results ==== |
The grid-search exploration result can be displayed using the python script make_grids.py : | The grid-search exploration result can be displayed using the python script make_grids.py : | ||
<code> | <code> | ||
Line 178: | Line 177: | ||
</code> | </code> | ||
- | The optional parameter –b activates the use of the basemap toolkit in order to draw the topography and coastlines. The resulting figures shown here in Fig. 1 and Fig. 2 are written in the pdf files grid_search_ts.pdf and grid_search_xy.pdf. | + | The optional parameter –b activates the use of the cartopy toolkit in order to draw the topography and coastlines. The resulting figures shown here in Fig. 1 and Fig. 2 are written in the pdf files grid_search_ts.pdf and grid_search_xy.pdf. |
There are two script which allows to compare predicted and observed waveforms. | There are two script which allows to compare predicted and observed waveforms. | ||
Line 189: | Line 188: | ||
The second script plot the waveform fit: | The second script plot the waveform fit: | ||
<code> | <code> | ||
- | ${WPHASE_HOME}/bin/traces_global.py | + | ${WPHASE_HOME}/bin/traces.py |
</code> | </code> | ||
- | The resulting figure is given in file wp_pages.pdf which is shown here in Fig. 4. If the basemap python module is available, a map is drawn showing the station location with respect to the centroid location. Otherwise, a simple polar representation is used to display the station location. | + | The resulting figure is given in file wp_pages.pdf which is shown here in Fig. 4. If the cartopy python module is available, a map is drawn showing the station location with respect to the centroid location. Otherwise, a simple polar representation is used to display the station location. |