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/08/04 13:55] wphase |
wphase:documentation [2018/05/03 16:55] wphase [Inversion] |
||
---|---|---|---|
Line 29: | Line 29: | ||
- You have to install numpy, matplotlib, basemap and netCDF4 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, [[wphase:macos|you can refer to this page]]. | To install these dependencies on MacOs, [[wphase:macos|you can refer to this page]]. | ||
Line 36: | Line 35: | ||
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 126: | 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 148: | 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 185: | 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 ==== | + | ===== Using ETOPO01 Global Relief Model ===== |
ETOPO01 is a 1 arc-minute global relief model of the Earth's surface | ETOPO01 is a 1 arc-minute global relief model of the Earth's surface |