User Tools

Site Tools


wphase:documentation

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
wphase:documentation [2015/01/23 19:54]
wphase [How to run W-phase]
wphase:documentation [2018/05/03 16:56]
wphase [RUNA scripts]
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. 
-  - GMT4 +  - rdseed
- +
-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 rms 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>​
 +
  
 ------- -------
Line 194: Line 206:
  
 ===== Data formats ===== ===== Data formats =====
 +
 +You can also have a look to [[wphase:​tutorial|the W-phase tutorial page]].
  
 === CMTSOLUTION FILE === === CMTSOLUTION FILE ===
wphase/documentation.txt · Last modified: 2022/01/10 07:39 by wphase