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 [2015/01/23 18:39] wphase |
wphase:documentation [2018/05/03 16:54] wphase [RUNA scripts] |
||
---|---|---|---|
Line 2: | Line 2: | ||
====== W-phase documentation ====== | ====== W-phase documentation ====== | ||
+ | In addition to this page, you should probably have a look to [[wphase:tutorial|this tutorial page]]. | ||
===== Installation ===== | ===== Installation ===== | ||
Line 7: | Line 8: | ||
==== Getting the code ==== | ==== 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. | + | Currently, the W-phase package is hosted as a github repository. To check out the W-phase repository: |
- | + | ||
- | Check out W-phase repository: | + | |
<code> | <code> | ||
git clone https://github.com/eost/wphase.git wphase_package | git clone https://github.com/eost/wphase.git wphase_package | ||
Line 20: | Line 19: | ||
</code> | </code> | ||
- | For more details, you can also checkout [[wphase:repository|this wiki page.]] | + | For more details on using git for W-phase, [[wphase:repository|you can read this page.]] |
==== Dependencies ==== | ==== Dependencies ==== | ||
Line 27: | Line 26: | ||
- csh shell | - csh shell | ||
- gcc and gfortran | - gcc and gfortran | ||
- | - python2.7 | + | - python2.7 (or later) |
- | - You have to install numpy, matplotlib and basemap to run some python scripts which make figures. | + | - You have to install numpy, matplotlib, basemap and netCDF4 to run some python scripts which make figures. |
- | - GMT4 | + | - rdseed |
+ | To install these dependencies on MacOs, [[wphase:macos|you can refer to this page]]. | ||
==== Building the code ==== | ==== 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: | ||
<code> | <code> | ||
- | setenv GMT_BIN /path/to/gmt/bin | ||
setenv RDSEED /path/to/rdseed/executable | setenv RDSEED /path/to/rdseed/executable | ||
setenv GF_PATH /path/to/greens/functions/database | setenv GF_PATH /path/to/greens/functions/database | ||
- | setenv WPHASE_HOME /path/to/wphase/package | + | setenv WPHASE_HOME /path/to/wphase |
</code> | </code> | ||
If you use bash: | If you use bash: | ||
<code> | <code> | ||
- | export GMT_BIN=/path/to/gmt/bin | ||
export RDSEED=/path/to/rdseed/executable | export RDSEED=/path/to/rdseed/executable | ||
export GF_PATH=/path/to/greens/functions/database | export GF_PATH=/path/to/greens/functions/database | ||
Line 69: | Line 67: | ||
===== How to run W-phase ===== | ===== How to run W-phase ===== | ||
- | You can checkout [[wphase:tutorial|this wiki tutorial page]]. | + | You should also have a look to [[wphase:tutorial|the W-phase tutorial page]]. |
==== Preparing directories, i_master and CMTSOLUTION ==== | ==== Preparing directories, i_master and CMTSOLUTION ==== | ||
Line 76: | Line 74: | ||
* i_master: a file containing other parameters such as the band-pass filter parameters, minimum and maximum distances, etc. | * i_master: a file containing other parameters such as the band-pass filter parameters, minimum and maximum distances, etc. | ||
- | The format of these files are described in [[wphase#Notes on file formats|notes on file formats]]. | + | The format of these files are described in [[wphase:documentation#Data formats|the file formats section]]. |
Path to SEED file(s) must be correctly specified after the keyword 'SEED:' in the i_master file. If multiple SEED files are used for the same inversion, each of them must be referenced properly in i_master using one 'SEED:' line per file. | Path to SEED file(s) must be correctly specified after the keyword 'SEED:' in the i_master file. If multiple SEED files are used for the same inversion, each of them must be referenced properly in i_master using one 'SEED:' line per file. | ||
Line 147: | Line 145: | ||
will perform median and rms misfit screening for Z channels only. | will perform median and rms misfit screening for Z channels only. | ||
- | These two scripts don’t perform the “ratio screening” (i.e. 3. above). To run the same routines with an additional screening based on the ratio (observed/predicted and predicted/observed), the following scripts can be used: | + | These two scripts don’t perform the “ratio screening” (i.e. 3. above). To run the same routines with an additional screening based on the amplitude ratio (observed/predicted and predicted/observed), the following script can be used: |
<code> | <code> | ||
${WPHASE_HOME}/bin/RUNA3r.csh | ${WPHASE_HOME}/bin/RUNA3r.csh | ||
Line 184: | Line 182: | ||
* The first script is <code>${WPHASE_HOME}/bin/make_grids.py</code> which be used to display grid-searches results. Use <code>${WPHASE_HOME}/bin/make_grids.py -h</code> to have more details on available options and arguments. The 2 other scripts plot observed and synthetic seismograms after W phase inversion. | * The first script is <code>${WPHASE_HOME}/bin/make_grids.py</code> which be used to display grid-searches results. Use <code>${WPHASE_HOME}/bin/make_grids.py -h</code> to have more details on available options and arguments. The 2 other scripts plot observed and synthetic seismograms after W phase inversion. | ||
- | * In order to plot complete seismograms individually and place station on a map: <code>${WPHASE_HOME}/bin/traces_global.py</code> which draw complete seismograms and show station location on a map (if basemap is available). Please use <code>${WPHASE_HOME}/bin/traces_global.py -h</code> for more details on available options. | + | * In order to plot complete seismograms individually and place station on a map: <code>${WPHASE_HOME}/bin/traces.py</code> which draw complete seismograms and show station location on a map (if basemap is available). Please use <code>${WPHASE_HOME}/bin/traces.py -h</code> for more details on available options. |
* To plot concatenated waveforms:<code>${WPHASE_HOME}/bin/make_cwp.py</code> For more details on available options and arguments see <code>${WPHASE_HOME}/bin/make_cwp.py -h</code> to see the options. | * To plot concatenated waveforms:<code>${WPHASE_HOME}/bin/make_cwp.py</code> For more details on available options and arguments see <code>${WPHASE_HOME}/bin/make_cwp.py -h</code> to see the options. | ||
+ | |||
+ | ------- | ||
+ | |||
+ | ===== Using ETOPO01 Global Relief Model ===== | ||
+ | |||
+ | ETOPO01 is a 1 arc-minute global relief model of the Earth's surface | ||
+ | provided by Amante et al. (2009): | ||
+ | <http://www.ngdc.noaa.gov/mgg/global/global.html> | ||
+ | |||
+ | If Basemap is installed, you can optionally draw ETOPO01 topography and | ||
+ | bathymetry. To do so, the path to the ETOPO01 netCDF file must be | ||
+ | assigned to the environment variable $ETOPOFILE which can be | ||
+ | done in your .tcshrc fil (or .bachrc, .cshrc, etc.). The make_grids.py | ||
+ | script have been tested only for grid-registered netCDF file of the | ||
+ | ETOPO1 bedrock file available at: | ||
+ | <http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/bedrock/grid_registered/netcdf/ETOPO1_Bed_g_gmt4.grd.gz> | ||
+ | |||
------- | ------- | ||
Line 191: | Line 206: | ||
===== Data formats ===== | ===== Data formats ===== | ||
+ | |||
+ | You can also have a look to [[wphase:tutorial|the W-phase tutorial page]]. | ||
=== CMTSOLUTION FILE === | === CMTSOLUTION FILE === |