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/08/04 13:55]
wphase [Plot routines]
wphase:documentation [2018/05/03 16:57]
wphase [RUNA scripts]
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 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 178: Line 175:
  
 ----- -----
-====== Plot routines ​======+==== Plot routines ====
  
 All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module basemap is also needed for plotting maps but it is optional (even if we recommend to install it). All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module basemap is also needed for plotting maps but it is optional (even if we recommend to install it).
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
wphase/documentation.txt · Last modified: 2018/05/03 16:57 by wphase