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 [2021/10/22 18:42] wphase [RUNA scripts] |
wphase:documentation [2021/10/22 18:51] wphase |
||
---|---|---|---|
Line 27: | Line 27: | ||
- gcc and gfortran | - gcc and gfortran | ||
- python2.7 (or later) | - python2.7 (or later) | ||
- | - You have to install numpy, matplotlib, basemap and netCDF4 to run some python scripts which make figures. | + | - You have to install numpy, matplotlib, cartopy and netCDF4 to run some python scripts which make figures. |
- rdseed | - rdseed | ||
Line 90: | Line 90: | ||
==== Extracting data from a gzipped tar archive ==== | ==== Extracting data from a gzipped tar archive ==== | ||
- | As SEED volumes are no longer supported, waveforms and instrument response can also be retrived from a gzipped tar archive (hereafter referred to as "tgz") including SAC and SAC_PZ (poles and zeros) files. Instead of using the "extract.csh" script to extract data from SEED, we can then extract data from a tgz archive using: | + | As SEED volumes are no longer supported, waveforms and instrument response can also be extracted from a gzipped tar archive (hereafter referred to as "tgz") including both SAC and SAC_PZ (poles and zeros) files. |
+ | Instead of using the "extract.csh" script to extract data from SEED, we can then extract data from a tgz archive using: | ||
<code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz</code> | <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz</code> | ||
- | where "path_to_tgz_archive.tgz" is the path to the tgz archive including SAC and SAC_PZ files. | + | where "path_to_tgz_archive.tgz" is the path to the tgz archive including SAC and SAC_PZ files. This script will extract SAC and SAC_PZ files a the directory named "DATA_org". |
Two tgz archive formats are supported: WILBER (default) and NIED. | Two tgz archive formats are supported: WILBER (default) and NIED. | ||
Line 187: | Line 188: | ||
==== 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 cartopy is also needed for plotting maps but it is optional (even if we recommend to install it). |
There are 3 script which can be used to plot the W phase inversion results: | There are 3 script which can be used to plot the W phase inversion results: | ||
* 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.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. | + | * 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 cartopy 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. | ||
Line 203: | Line 204: | ||
<http://www.ngdc.noaa.gov/mgg/global/global.html> | <http://www.ngdc.noaa.gov/mgg/global/global.html> | ||
- | If Basemap is installed, you can optionally draw ETOPO01 topography and | + | If Cartopy is installed, you can optionally draw ETOPO01 topography and |
bathymetry. To do so, the path to the ETOPO01 netCDF file must be | bathymetry. To do so, the path to the ETOPO01 netCDF file must be | ||
assigned to the environment variable $ETOPOFILE which can be | assigned to the environment variable $ETOPOFILE which can be | ||
Line 391: | Line 392: | ||
${WPHASE_HOME}/bin/make_grids.py –p -b | ${WPHASE_HOME}/bin/make_grids.py –p -b | ||
</code> | </code> | ||
- | to create the file grid_search_xy.pdf containing Fig. 1 which corresponds to the above example. The option –b is optional: it allows to use the basemap module in order to plot coastlines and topography. | + | to create the file grid_search_xy.pdf containing Fig. 1 which corresponds to the above example. The option –b is optional: it allows to use the cartopy module in order to plot coastlines and topography. |