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
wphase:documentation [2015/01/27 18:01]
wphase [Dependencies]
wphase:documentation [2018/05/03 16:57] (current)
wphase [RUNA scripts]
Line 27: Line 27:
   - gcc and gfortran   - gcc and gfortran
   - python2.7 (or later)   - 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, matplotlibbasemap ​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 138: Line 135:
 ----- -----
 ==== RUNA scripts ==== ==== RUNA scripts ====
-Data extraction, screening and inversion described above can be performed by running one of the RUNA3*csh scripts (usually in bin). These scripts perform data extraction/​screening and calculate a moment tensor solution after median data screening and after rejecting ​ traces associated with large misfit using the threshold 5.0 3.0 0.9 (i.e. -th, see 2. in section ​4):+Data extraction, screening and inversion described above can be performed by running one of the RUNA3*csh scripts (usually in bin). These scripts perform data extraction/​screening and calculate a moment tensor solution after median data screening and after rejecting ​ traces associated with large misfit using the threshold 5.0 3.0 0.9 (i.e. -th, see 2. in the "​Inversion" ​section ​above):
 <​code>​ <​code>​
 ${WPHASE_HOME}/​bin/​RUNA3.csh ${WPHASE_HOME}/​bin/​RUNA3.csh
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. in the "​Inversion"​ section ​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 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 =====
 +
 +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>​
 +
  
 ------- -------
wphase/documentation.1422378118.txt.gz · Last modified: 2015/01/27 18:01 by wphase