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/27 17:39] wphase [Dependencies] |
wphase:documentation [2018/05/03 16:55] wphase [Inversion] |
||
---|---|---|---|
Line 8: | 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 28: | 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. |
- rdseed | - rdseed | ||
- | - GMT4 | ||
- | To install these dependencies on MacOs, you can refer to [[wphase:macos|this page]]. | + | 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 128: | Line 123: | ||
$WPHASE_HOME/bin/wpinversion -ifil o_wpinv.th1.3 -ofil o_wpinv.th0.9 -th 0.9 | $WPHASE_HOME/bin/wpinversion -ifil o_wpinv.th1.3 -ofil o_wpinv.th0.9 -th 0.9 | ||
</code> | </code> | ||
- | - Finally, it is possible to remove the channels showing a large rms ratios (observed/predicted and predicted/observed) using the option –nr (e.g. –nr 2.0): <code> | + | - Finally, it is possible to remove the channels showing a large rms amplitude ratios (observed/predicted and predicted/observed) using the option –nr (e.g. –nr 2.0): <code> |
$WPHASE_HOME/bin/wpinversion -ifil o_wpinv.th0.9 -ofil o_wpinv.th1.3 -nr 2.0 | $WPHASE_HOME/bin/wpinversion -ifil o_wpinv.th0.9 -ofil o_wpinv.th1.3 -nr 2.0 | ||
</code> | </code> | ||
Line 150: | 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 187: | 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> | ||
+ | |||
------- | ------- |